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VORTRAGSBERICHTE 
A. Praktische Analysis 


Ergebnisse neuer Versuche zur Periodenanalyse, die im Istituto 
Nazionale per le Applicazioni del Calcolo ausgeführt wurden 


Von Mauro Picone und Fabio Conforto in Rom 


Sei f(t) eine empirisch gegebene reell- oder komplexwertige Funktion im Intervalle 
—1<st<s-H1. Die Periodenanalyse der Funktion f(£) bedeutet hier folgendes: Bestim- 
mung derersten natürlichen Zahl nund der Konstanten 71; as --+» Ins Zi ?as -++> Zu» 
so daß die Ungleichung 


ke z,t 
ft) =. DR NE RER a (1) 
gilt. 
Die Zahl e ist dabei eine bekannte positive Konstante, die mit der Genauigkeit der Ver- 
suche zusammenhängt, so daß der wirkliche Wert der Funktion f(t) zwischen f(t) —e und 
fi) + & liegt. 

Im folgenden wird eine Methode zur Periodenanalyse kurz skizziert. Zuerst sei folgendes 
bemerkt: Da f(t) nur empirisch gegeben ist, darf man nie mit Ableitungen von f(f) rechnen. Wohl 
kann man dagegen von bestimmten Integralen von f(t) oder des-Produktes von f(t) mit bekannten 
Funktionen sprechen. Die Zahl » darf nie zu groß sein. Die Methode ist nun die folgende: 

Hätte man exakt 


ge zyt 
f (t) 7 & Yne ’ 
h-1 
so wäre f(t) Lösung einer Differentialgleichung mit konstanten Koeffizienten: 
N 
24/P)—0 Ba N: (2). 
h=0 


Ist nun 9,(t) (s=1, 2, ...) ein im Hilbertschen Sinne vollständiges System von Funktionen 
im Intervalle [—1, +1], so ist die Gl.(2) äquivalent den Gleichungen: 


x n 
[ed af” () di= 0 1,2...) 
— h=0 


die man nach einer gewissen Anzahl von partiellen Integrationen auch folgendermaßen schreiben 


kann: 
1 


Sen). at rd at n 9m (fl) a — (1) [pw (ft) di 5 


a 
Kein.) 
‚mit 
Ms n 
2, Al) Pi, — (1), fr? N)=a,. 
h=i h=i 
_ Ist aber f(t) nur empirisch gegeben, so kann man ohne Weiteres die Integrale in (3) berechnen, 


- aber nicht die Konstanten a, und f;, dain diesen auch Werte der Ableitungen von f{t) fürt= +1 
16 
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vorkommen. Man betrachtet daher einfach die &;, ß; in (3) als neue Unbekannte. Formel (3) 
ist dann ein System von unendlich vielen linearen Gleichungen mit 3%» Unbekannten, nämlich 
BB Ü=L 2% %..,.nprundYo, (Rh Ober, NZ). x 


Die Koeffizientenmatrix des Systems (3) hat genau den Rang 3n, wenn z.B. 9,(t) = est 
gesetzt wird und die Folge der A, die Eigenschaft hat: 
Maxlimi,=-+ o, Min lim A, = —- wolf sa 


so» so» 


vorausgesetzt, daß f(t) nicht schon einer linearen homogenen Differentialgleichung mit kon- 
stanten Koeffizienten einer Ordnung m <n genügt. 

Die Folge der A, kann z.B. die Folge 0, —1, +1, —2, +2, —3, ... sein. Zusammen- 
fassend kann man also für jedes empirisch gegebene f(t) das System (3) Bi mit 9, (t)—=e est, 
wobei für die Folge der A, die Gl.(4) gelten. 

Zahlreiche numerische Versuche, die i in Rom im Istituto Nazionale per le Applicazioni del 
Calcolo ausgeführt wurden, haben nun gezeigt, daß das mit der Methode der kleinsten Quadrate 
gelöste System einer zunehmenden Anzahl p von Gl.(3) tatsächlich die richtigen Werte der 
%(h=0,1,..., n—1) liefert für immer größer werdendes 7, falls solche Werte bekannt sind. 

Sind daher die o,(k=0, 1, ...., n—1) unbekannt und f(t) empirisch gegeben, so liegt 
es nahe, das System von p Gleichungen des Systems (3) mit der Methode der kleinsten Quadrate 
zu lösen fürn=1,2,.... Ist n die erste Zahl, für welche für zunehmendes 9 die co, gegen 
bestimmte Konstanten c, konvergieren, so kann man vermuten, daß f(t) sehr nahe an einem 
Integral der mit diesen c, gebildeten Differentialgleichung (2) ist. Man findet daher die Wurzeln 
z(h=1, 2, ...., n) der charakteristischen Gleichung dieser Differentialgleichung und bestimmt 
dann die y, durch Lösung mit der Methode der kleinsten Quadrate einer zunehmenden Anzahl 9 
von Gleichungen des De 


1 
Je #70 ul. = 7 Jet a ENE 


Konvergiert für zunehmendes p die Lösung dieses Systems gegen gewisse y7, so sind die y, und 2, 
bestimmt, und man kann in numerischer Weise probieren, ob (1) erfüllt ist. Ist (1) erfüllt, so 
ist die Periodenanalyse erledigt; im entgegengesetzten Falle kann man mit n+]1, und dann 
mit n-+2 usw. probieren. Wird aber n zu groß, so muß man eine Periodenanalyse für f(t) als 
unmöglich betrachten. 


Einige neuere Forschungen über numerische Behandlung 
von Differentialgleichungen 


Von L. Oollatz in Hannover 


Hier können nur die Ideen einiger neuerer Methoden angedeutet werden; wegen ausführlicherer Dar- 
stellungen muß auf die Literatur verwiesen werden. Vollständigkeit ist nicht angestrebt. 


I. Ausbau einiger spezieller Methoden 
1. Differenzenverfahren 


a) Iterative Lösung. In letzter Zeit wurden die iterativen Methoden zur Lösung der bei 
kleiner Maschenweite umfangreichen Gleichungssysteme, insbesondere Relaxationsmethoden, viel 
verwendet und ausgebaut?). 


b) Nichtiterative Lösung. Hyman?) hat eine Methode aufgestellt, die am Beispiel der 

1. Randwertaufgabe der Potentialtheorie bei einem Rechtecksbereich O<x<sa, 0 <y<sb 

erläutert werde. Es seien für eine Lösung u(z, y) der Potentialgleichung Au=0 die Randwerte 

auf dem Rechtecksrande /’vorgegeben. Bei einem Gitter ll, den Maschenweiten h,l und den 
Gitterpunkten 

3 a b x (j = 0, .... 

,=jh, Y: — kl, een ee 0.1. od) 


1) Aus der umfangreichen Literatur sei nur genannt: R. v. Southwell: Relaxation methods in Ent. 
neering Science, Oxford 1934, Relaxation methods in Theoretical Physics, Oxford 1946; L. Fox: Some i impro- 
vements in the use of relaxation methods, Proc. roy. Soc. London 190 (1947), 8. 31 — 59; G.M. Dusinberre: 
Numerical Analysis of Heath Flow, New York, Toronto, London 1949. 

2) M. Hyman: On the Non-Iterative Numerical Solution of Boundary Value Problems. Vortrag 
Washington, April 1950. e 


t ce gleichungen. Vortrag Freiburg 30. März 1951. Z. angew. Math. Mech. 31 (1951), 5. 240— 241. , 
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haben die Differenzengleichungen Partikularlösungen der Form 
- k 
A N, re STE ER ra 
wobei A aus einer quadratischen Gleichung zu bestimmen ist; man findet 
1 ]\2 
ui mit p=1-—(}) (cos "®*_1). 
2n h a 
Durch Überlagerung erhält man eine Partikularlösung 


M 
> sin ut (P, An En Q, Man) ’ 
n=1 q 

welche auf den Rechtecksseiten &=0 und 2=a verschwindet und auf y=0 und y=b die 
vorgegebenen Randwerte annimmt, sofern man die P, und Q, passend bestimmt (dazu ist nur 
harmonische Analyse nötig, die nach den bekannten Verfahren von Runge usw. durchgeführt 
werden kann). Durch Überlagerung mit der entsprechenden Partikularlösung, die auf x&—=0 
und x—=a die vorgegebenen Randwerte annimmt, gewinnt man die exakte Lösung derDifferenzen- 
gleichungen. Das Verfahren ist auch bei anderen Differentialgleichungen mit konstanten Koeffi- 
zienten (Potentialgleichung in drei Dimensionen, Plattengleichung usw.), bei der 2. und 3. Rand- 
wertaufgabe und bei Hinzunahme der Kollokation (vgl. hier Nr. 2a) auch bei krummlinigen 
Rändern anwendbar. 

c) Mehrstellenverfahren. Bei der linearen Differentialgleichung m-ter Ordnung 


Z[ ] y 4 0P ta, ? 3 
1 ya \% Sm art, DO RATEN Sr 
ich 2a ( er: ayı ( Y%) ( ) 
macht man bei Zugrundelegung eines etwa rechteckigen Gitters wie in (1) einen Ansatz: 
9=2 USE Dur eh ne un en LE (A 


wobei u, bzw. (L [w]),. den Wert von « bzw. L[w] an der Gitterstelle &,, y, bedeutet, die Summe 
über eine noch zu wählende Anzahl von Gitterpunkten zu erstrecken ist und die C,. und D,. 
so zu bestimmen sind, daß beiTaylorentwicklung von ® nachu und den partiellen Ableitungen 
an einer festen Gitterstelle ı,, x, die Faktoren von % und deren Ableitungen bis zu möglichst 
hoher Ordnung verschwinden; man erhält dann eine Näherungsgleichung ®=0, wenn man 4,, 
durch Näherungswerte U,, und L[w] durch r ersetzt. Das gewöhnliche Differenzenverfahren 
ordnet sich hier ein, wenn man nur eine Größe D,,, verwendet und alle anderen D,„—0 setzt. 

Das Mehrstellenverfahren®) hat bei einfach gebauten Differentialgleichungen gegenüber 
dem gewöhnlicheu Differenzenverfahren den Vorzug, bei nur unwesentlich vergrößerter Rechen- 
arbeit im allgemeinen zu genaueren Ergebnissen zu führen. Das Verfahren wurde auf ver- 
schiedene Differentialgleichungen angewendet, z. B. auf partielle Differentialgleichungen 4. Ord- 
nung der Plattenbiegung®), auf parabolische und hyperbolische Differentialgleichungen; um nur 
einen Ausdruck zu nennen, sei die inhomogene Wärmeleitungsgleichung 


ou, u 
In es =— pr 
oy  9a% 
herausgegriffen. In dem Gitter (1) erhält man für u bei konstantem K und r die Beziehung 
l 
User, EAU/ESR Zu a + I; 2m en K' + Restglied 4. Ordnung . . (5). 


Bei komplizierteren Differentialgleichungen wird die Durchführung des Mehrstellenver- 
fahrens mühsam; z.B. hätte man an sich bei 


U nla)yn@)—r() 


für jede Ableitung y” (x) die Ausdrücke des Mehrstellenverfahrens aufzustellen und erhielte bei m 
Abszissen ein Gleichungssystem für (größenordnungsmäßig) n »m Unbekannte. Durch Um- 
schreibung der Differentialgleichung in eine Form, wo an Stelle der Ableitungen wiederholte 
Integrale über y(x) mit gewissen Faktoren auftreten, gelingt es Sassenfeld°), auch hier ein 
Gleichungssystem für nur m Unbekannte aufzustellen. 


3) Eine allgemeine Darstellung bei L. Collatz: Numerische Behandlung von Differentialgleichungen. 


Berlin 1951, S. 116, 235, 311. 
4) Z. angew. Math. Mech. 30 (1950), S. 385—388. REN 
5) Sassenfeld: Ein Summenverfahren für Rand- und Eigenwertaufgaben linearer Differential- 
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d) Verallgemeinerung des Mehrstellenverfahrens. Der: Ansatz (4) läßt sich erweitern zu 


P=2 {Cru + Dr Ed} +2 (MH). A Et: 


wobei M, gewählte, also gegebene Operatoren sind. In der Näherungsgleichung ®—=0 ist dann 
auch M,„[r] bekannt. Durch Hinzunahme derartiger Glieder läßt sich ein Abgleich in der 
Taylorentwicklung bis zu höheren Ableitungen erzielen. Die Anwendung der so entstehenden 
genaueren Formeln hat natürlich nur Sinn, wenn die N Ableitungen genügend 
Ru 

02° Ba oy? 

angegeben, wobei als Operator M wieder A gewählt ist und ein quadratisches Gitter der Maschen- 
weite A angenommen ist: 


D®— 20 u, — 4 — du, 2 TS - (89 Au, Rn; %+L%Au,) ER ar 704 Awu,=Restglied8.Ordnung, 


hoher Ordnung besitzt. Als Beispiel seien für Lan zwei Ausdrücke ® 


d—=-20u,-48%W-8 ut = (34 Au,— dA ee &0 “(34 AA 2 3A Au,)— Restglied 8. Ordnung. 


Dabei ist a ein beliebiger Gitterpunkt, %, der Wert von u an dieser Stelle, I, bedeutet die 
Summe der #-Werte in den vier Nachbarpunkten, &w, die Summe der u-Werte in den vier 
vom Punkte a um Ah y2 entfernten Gitterpunkten, entsprechend sind & Au, usw. zu verstehen. 

e) Gewöhnliche Differentialgleichungen im Komplexen. Salzer®) gibt für Differential- 
gleichungen 1. und 2. Ordnung für eine gesuchte Funktion w(2) Formeln an, nach denen man 
in einem quadratischen Gitter in der 2-Ebene aus den z-Werten in einigen Gitterpunkten die 
z-Werte in weiteren Gitterpunkten angenähert berechnen kann. 


2. Fehlerabgleichsmethoden 


Für eine Funktion u(&,, ..., %,) von n unabhängigen Veränderlichen seien eine Diffe- 
rentialgleichung 
F(&, ... Ins %; %ıs ne, Uns %1> 22} Unns ” .)=0 in B Er Br (7) 


und Randbedingungen 
Vals Da a U ee Uns el RAU I a DE ee 


vorgelegt (die 7‘, können ganz oder teilweise zusammenfallen). Bei « und in (10) bei » be- 


. N u 
deuten Indizes partielle Ableitungen, z.B. u. = ERSEYSE 
; 9% 
&,*** ©,-Raumes, 7‘, sind (n—1)-dimensionale ‚‚Hyperflächen“, F und V,, sind gegebene Funk- 
tionen ihrer Argumente; ein Näherungsansatz für % 


BEmIGmN . ms O re ee 2a DR REN! (9) 


sei so gewählt, daß er für beliebige a. entweder (Fall 1) die Differentialgleichung oder (Fall 2) 
die Randbedingungen exakt erfüllt. Die a, sind dann so zu bestimmen, daß im Fall 1 die Rand- 
bedingungen und im Fall 2 die Differentialgleichung möglichst gut erfüllt werden. Im Fall 1 
spricht man von einer Randmethode, im Fall 2 von einer Gebietsmethode. Es werde etwa 
Fall 2 genauer beschrieben, in Falll hat man Differentialgleichung und Randbedingungen sinn- 


gemäß zu vertauschen. Beim Einsetzen von w in die Differentialgleichung bleibt eine Fehler- 
funktion e übrig 


; 8 ist ein gegebener Bereich des 


EL, 72 En  ) Flair. 
Die a, sollen so bestimmt werden, daß die Fehlerfunktion g möglichst gut die Funktion Null 
annähert. Von den verschiedenen Möglichkeiten seien einige genannt: 

a) Kollokation. Der Fehler e soll in p Punkten P,, ..., P,, den „Kollokationspunkten“ 
verschwinden. Die p Punkte versucht man einigermaßen gleichmäßig über den Bereich ® zu 
verteilen. Sind 219; » .., Zne die Koordinaten von P,, so lauten die Bestimmungsgleichungen 
für die @,: 

E (dig +. ner A pH )=O (= 1, 1 PP N REN! 
Es wäre wichtig, über die zweckmäßige Wahl der Kollokationspunkte Untersuchungen an- 
zustellen. 

b) Kollokation mit Ableitungen. Die Rechnung wird zuweilen einfacher, wenn man nicht & 
an p Kollokationspunkten gleich Null setzt, sondern nur an g Stellen, an einigen Stellen (den- 


6) H. E. Salzer: Formulas for numerical De of first and second order differential equations 
in the complex plane. J. Math. Physics 29 (1950), S. 207—216. 


EI ale Math, aa \ ; 
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selben oder anderen Stellen) das Verschwinden von Ableitungen von & verlangt, und so insgesamt 
p Bestimmungsgleichungen aufstellt. 
c) Gewöhnliche Fehlerquadratmethode. Man verlangt 


=: [da ...da, = Min. AT ER TAT RE 


AH 
und erhält in PR —= 0 insgesamt p Bestimmungsgleichungen für die «,. 


d) Fehlerquadratmethode mit Gewichtsfunktion?). Es sei p(®,, ...., 2,) eine in B positive 
Funktion. An Stelle von (12) tritt 


S-: „[pedan. ER ih a N) 

e) Fehlerorthogonalität. Der Fehler & soll zu p gewählten, voneinander linear unab- 
hängigen Ortsfunktionen g,, ..., 9, orthogonal sein: 

I: [en da ...da,— 0 En. a a a FAR 


-e*) Galerkinsche Gleichungen. Sind die Randbedingungen linear, so kann man für w 
einen Ansatz verwenden 


p 
N ED er ae DR 
01 
wobei v, die gegebenen (inhomogenen) Randbedingungen und v,, ..., vd, die zugehörigen 
4 homogenen Randbedingungen erfüllen. Für v,—g, ergeben sich aus (14) die Galerkinschen 
Gleichungen. 
f) Teilgebietsmethode. Der Bereich ® wird in p Teilbereiche ®,, ..., ®, aufgeteilt. 
Dann soll gelten 
„Sedan...d—0 BEE EIH DIE ee RELON 
e 


Alle hier als Gebietsmethoden formulierten Verfahren: sind, wie oben betont, auch als 
Randmethoden anwendbar. Die Randmethoden haben gegenüber den Gebietsmethoden den 
Vorzug, zu ihrer Durchführung nur die Auswertung von Randintegralen, aber nicht von Gebiets- 
integralen zu erfordern. (Bei der Kollokation sind keine Integrale auszuwerten, aber auch dort 
ist die Randkollokation angenehmer und weniger unsicher als die Gebietskollokation.) Für alle 
hier genannten Verfahren ist eine weitere praktische Erprobung und theoretische Durchforschung 
sehr erwünscht. 

Für die drei wohl noch wenig verwendeten Methoden a), b), f) sei ein ganz einfaches Beispiel 
ER) werde der 
nl 
die Randbedingungen erfüllende Ansatz w= x + a, (& — x) + a, (@° — x) gemacht; dann 
wirde=w’+w=z-+a,(®— 2 -+2)+a,(2°-+-5 x). Bei Methode a) (Kollokation) setzt 
man &(2,) = &(2,) = 0 und bestimmt daraus a, und as; die Tabelle gibt die Werte für mehrere 


gegeben. Bei y’’+y=0; y(0)=0, y(1)=1 (mit der exakten Lösung y(x) = 


1 
Wahlen von &, und &5; beiMethode b)wird man etwa e und eg’ gleich 0 setzen an der Stelle 2 


2) ’ 
> und bei f) liegt es nahe, 
’ 1/2 1 
edt— f edz=0 
fi) 1/2 
zu verlangen. 
1 1 
4, @y | w (5) Fehler von w (3) 
$ 2 2 


+ 0,9% 
0,56897 — 0,14% 


&) Kollokation 


+ 010% 


?) M. Picone: Analisi quantitativa ed esistenziale nei problemi di propagazione, Atti del 1° Congresso 
> dell’Uniono Mathematica Italiana (1937). 
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3. Ritzsches und Trefftzsches Verfahren 

a) Gemischter Ritz-Ansatz?). 

Es sei J[w] ein Integralausdruck (wie beim Ritzschen Verfahren), der für die Lösung « 
einer vorgelegten, etwa linearen Randwertaufgabe einen stationären Wert, z.B. ein Minimum 
annimmt. Beim Ritzschen Verfahren würde man für w den Ansatz (15) machen. Nun möge 
für eine Veränderliche, etwa x,, in den Randbedingungen keine Differentiation nach %, vor- 
kommen. Dann erfüllt w in (15) auch die Randbedingungen, wenn man die a, als noch zu be- 
stimmende Funktionen von x, ansieht. Die Forderung 


J[Iw]=J.{a, (2,)0==Min.o 2 Er. SiRene re | 


führt dann zu einem System Eulerscher Differentialgleichungen, und zwar zu p gewöhnlichen 
Differentialgleichungen für die a,(%,). 

b) Trefftzsches Verfahren. 

Während das Ritzsche Verfahren eine Gebietsmethode ist, ist das Trefftzsche Verfahren 
eine Randmethode und fand als solche in letzter Zeit wieder mehr Beachtung. Ursprünglich 
von Trefftz nur für die 1. Randwertaufgabe der Potentialtheorie aufgestellt, wurde das Trefftz- 
sche Verfahren auf allgemeinere Aufgaben ausgedehnt und ist auch bei der 2.und 3. Randwert- 
aufgabe bei der elliptischen selbstadjungierten Differentialgleichung 


4 2 ou 
U=— D zu lAıg2)tmu=r: FE rd ad (18) 
7 2 


kl 


u are 


anwendbar. Dabei sind g, r gegebene, stetige und die A;; stetig differenzierbare Funktionen 
von 2%, ».., 2,; für zwei mit stetigen partiellen zweiten Ableitungen versehene Funktionen 9, 
y lautet dann die Greensche Formel 


eye nen (NM. 


Dabei ist I’ der (gewissen Glattheilsforderungen genügende) Rand eines Bereiches ®, df das 
Oberflächenelement auf 7‘, dr das Volumenelement in ® und ZL* bis auf einen Faktor der Diffe- 
rentiationsoperator in Richtung der Konormalen°), bei der Potentialgleichung ist die Konormale“ 
einfach die äußere Normale. 


Nun sei füf u die Randbedingung vorgelegt 
Au ZH ule Alan 1 ein Or eh . (20). 


Für « wird ein Näherungsansatz gemacht 


vi 
VEW—UT DEN N RER TR, 
e= 


> ie u 


EM T, 


mit L[[v]J=r, L[w]=0 für o=1, ..., p. Die Greensche Formel (19) liefert dann für 
yeU—v, y=y R 


u— vo) L*][ v,] — v, L* [u — vo]} df = 0, (el a En. (22). 


—_ 
- 


Bei der 1. Randwertaufgabe ist u am Rande gegeben; ersetzt man Z*[u—v,] angenähert 
P p . R 

durch Z*] 2 4,0] = I a, L*[v,], so sind die Gleichungen (22) p lineare Bestimmungsgleichungen 

. = 1 { 

für die a,. Bei der 2. Randwertaufgabe ist Z*[ «] auf I’ bekannt und man ersetzt in (22) in der 


p 
runden Klammer u—v, durch 3 a,v,; man erhält wieder p lineare Gleichungen für die a,. Bei 

7=T ö i 
der 3. Randwerlaufgabe kann man «—v, und L*[ u—v,] auf 7’ durch einander ausdrücken und 
so ebenfalls p lineare Gleichungen für die a, aufstellen !). | a 


®) L. Kantorovie: Sur une möthode direete de la solution approximative du probleme du minimum — 
d’une intögrable double (Russisch), Leningrad, Bull. Ac. Sc.7 (1933), 8. 647—652. E. Mettler: Allgemeine ae 
Theorie der Stabilität erzwungener Schwingungen elastischer Körper. Ing.-Arch. 17 (1949), S.418-—49,. 
F. Weidenhammer: Der eingespannte, axial pulsierend belastete Stab als Stabilitätsproblem. Z.angew.. 
Math. am), re > ge re 
) Vgl. A. @. Webster u. G. Szegö: Partielle Differentialgleie ı ischen ik. 
Leipzig und Berlin 1930, 8. 311. s 5 ee m gene a Ka 
10) Ähnliche Überlegungen bei M. Picone: Sur le caleul de la d6formation d’un solide 6lastique encastr6. 
tila 


7. Internat. Kongress für angewandte Mechanik. London 1948 und G. Fichera: Risultati concernenti® 
risoluzione delle equazioni funzionali lineari dovuti all’ Istituto Nazionale per le applicazioni del calcolc 
Memorie dell’ Accademia Nazionale dei Lincei. Sei. fis. mat. serie VIII, vol. III, sez. % fasc. 1. Ro, 
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II. Theoretische Durehforsehung einiger bekannter Methoden. 


1. Stabilität des Differenzenverfahrens 
a) Wellengleichung. Der Wellengleichung mit zwei Raumkoordinaten x, y% und der Zeit- 
koordinate 2 
u. u „o2Uu 
022 Ay2 022 ° 
(K gegebene Konstante) entspricht in einem quaderförmigen Gitter mit den Maschenweiten A, I 
und den Gitterpunkten 


. (23) 


Is =utih yzy=y+jih ey =uHtkl... . .'. (24) 
die Differenzengleichung 
Ui; = 0 (Virı,,a + Ui; + Ui ja + Ui,e) + (2—40) Oi, — Dia: - - (25). 


Dabei ist U,;, ;,; ein Näherungswert für u(z,, Y%,2,) und ?=Koh? gesetzt. Würde man nun 
o=1, also I=Kh? wählen, so würde (25) für die numerische Rechnung unbrauchbar werden; 
ein kleiner Fehler e in den U-Werten, wie er z. B. allein schon durch die Abrundung entsteht, 
würde sich so fortpflanzen, daß die Rechnung durch die anwachsenden Fehler völlig überdeckt 
werden würde (ein Fehler e an der Stelle :, 7, k würde an der Stelle ö, 7, k-H1 den Fehler — 2e, 
an der Stelle 2,7, +2 den Fehler 7& und an der Stelle &, 5, k-+3 bereits den Fehler — 28 & 
hervorrufen). Das Verfahren wird daher „instabil“ genannt. Wählt man jedoch e=1}, so 
lautet (25) 


1 
Unyen = 5 Wa Ve Vie FÜ 214) Vene 0 >. 4026). 


Ein Fehler & bei 2=2,+ kl verteilt sich jetzt bei größeren 2-Werten und kann nicht anwachsen; 
das Verfahren ist jetzt ‚‚stabil‘“ gegenüber kleinen Störungen. 

b) Einzelfehler. Den Einfluß eines Einzelfehlers e auf die weitere Rechnung untersucht 
John v. Neumann!) für ziemlich allgemeine lineare Differentialgleichungen unter Benutzung 
weitreichender mathematischer Hilfsmittel. 

c) Fehlerverteilung. Viel einfacher in der Durchführung (allerdings auch nicht so weit- 
reichend) ist eine andere Methode ?), nicht von einem Einzelfehler, sondern von einer Fehler- 
verteilung auszugehen. Die Methode ist in einer Arbeit!?) für die Wärmeleitungsgleichung und 
die Wellengleichung mit einer Raumdimension durchgeführt und werde hier für die Wellen- 
gleichung (23) mit zwei Raumdimensionen beschrieben. Die zugehörige Differenzengleichung (25) 
hat die Partikularlösung 


Ua: sinifsinjn - 2... 20... (27), 
wenn zwischen £&, n, & der Zusammenhang 
SE —1=0(os5sE+c08sNn—2). . .. 2.2.2.2. (28) 


besteht. Liegt etwa die Anfangs-Randwertaufgabe vor: 
| - ou 
u(0, y,2)= u(a, y,2)=u(2,0,2)=u(z,a,2)=0, u(z,y,0)=fı (8, Y), 5) Nena (2, Y), 
h U; 


mit gegebenem f, f» (quadratische, ringsum eingespannte Membran mit gegebenem Anfangs- 
zustand), und wählt man ein quadratisches Gitter in der &-y-Ebene mit der Maschenweite 
h=a/N, so kann manmit $=&, = = le Tr die Fehlerverteilung in der Ebene 2=1 


in der Form 


N—1 : ; 
._ UN . IV 3 
&,j1 > auesin 7 sin on a, ne (297 
uv—l N N 
ansetzen; man erhält als allgemeine Fehlerverteilung 
a : 
du UT a 
EEE >" sin "7" sin 7 (be wo, e Be 10 4-0 (30) 
u a \ 


11) Eine Darstellung gibt R. P. Eddy: Stability in the Numerical Solution of Initial Value Problems 

in Partial Differential Equations, Naval Ordnance Laboratory Memorandum, Bd. 10 (1949), S. 232. 

. 12 0. G. Brien, M. Hyman u. S. Kaplan: A Study of the Numerical Solutions of Partial Differential 
_ Equations; Journal of Mathematics and Physics 29 (1951), 8. 223—251, eine ausführlichere Darstellung der 
Methode soll am Beispiel der Stabschwingungsgleichung demnächst in einer kleinen Mitteilung in dieser Zeit- 
schrift gegeben werden. 
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wobei b,, und @,, sich in einfacher Weise aus den @,, berechnen. Stabilität besteht, wenn &;,,,; 
beschränkt bleibt für k > ©, d.h. wenn |e*!ur| = 1 ist für alle £,,, die sich mit &=£,, aus (28) 


für &=&,, nn, ergeben. Die Diskussion zeigt, daß el erfüllt ist für O<e s— y> 
aus diesem Intervall darf also oe gewählt werden. 


2. Fehlerabschätzungen 


a) Milderung der Lipschitzbedingung. Bei einer gewöhnlichen Differentialgleichung 


1. Ordnung y’=f(z, y) und bei Systemen 
s—F(® Yo». Yu) Gel... a.n nn ine) 


pflegte man bei Fehlerabschätzungen bisher stets die Lipschitzbedingung in der Form mit 
absoluten Beträgen zu verwenden. Eltermann?®) trifft eine mildere Voraussetzung: für festes & 
und zwei verschiedene Wertesysteme y; und y, genüge die Funktion 


P 15 @ ALIEN nF; (® 1» als Yn)] (4 — 5) 
L(& Du 2: In ya 90) = ; a) 


innerhalb eines betrachteten Bereiches der Bedingung 


Le 9. I Null SL ee 


insbesondere für negalives Z gelangt er zu wesentlich besseren Abschätzungen, als sie bisher 
vorlagen. 


b) Ritzsches Verfahren bei Eigenwerlaufgaben. Abschätzungen von N. Kryloff ver- 
bessert Bertram!#), zugleich für selbstadjungierte volldefinite Eigenwertaufgaben. 


c) Fehlerabschätzung nach v. Mises. Während man sonst bei einer Aufgabe mit der 
Differentialgleichung (7) und Anfangs- oder Randbedingungen (8) nach dem Fehler «&, d.h. der 
Abweichung w— u einer Näherungslösung w von der exakten Lösung « fragt, hält v. Mises"°) 
es vom physikalischen Standpunkt aus für vernünftiger, bei gegebener Toleranz ö nach einer 
Näherung w mit 


Fila, 8.2) V la en a (34) 


zu fragen. Diese Fragestellung verlangt zugleich die Aufstellung weiterer Näherungsverfahren *®). 


d) Fehlerabschätzung für Ableitungen. Bisher hat man bei einer Aufgabe (7), (8) ge- 
wöhnlich nach einer Abschätzung für den Fehlerbetrag |a| = |w»— w| der Näherungslösung w 
gefragt. Bei den Anwendungen interessiert man sich aber häufig mehr für den Fehler bei den 
Ableitungen, z. B. bei Plattenaufgaben ist der Fehler bei der Durchbiegung meist nicht so wichtig, 
als der Fehler bei den Spannungen, Momenten usw. In dieser Richtung ist erst wenig Er 
- worden. 


e) Satz vom Randmaximum. Der klassische Satz, daß eine in einem abgeschlossenen 
beschränkten Bereich reguläre nichtkonstante Potentialfunktion ihr Maximum und ihr Minimum 
nur auf dem Rande annimmt, ist auf allgemeinere elliptische Differentialgleichungen übertragbar 


und kann zu einer Fehlerabschätzung für die 1. Randwertaufgabe benutzt werden. Es fehlen 


noch die entsprechenden Sätze für die 2. Irand werkautzabezn: ” 


13) H. Eltermann: Verbesserung der Fehlerabschätzung bei Differentialgleichungssystemen, Vortrag 
Freiburg 30. März 1951. 

14) G. Bertram: Zur Fehlerabschätzung für das Ritzsche Verfahren bei Eigenwertaufgaben, Vortrag 
Freiburg 30. März 1951, Z. angew. Math. Mech. 31 (1951), S. 141. 


15) R. v. Mises: On Network Methods in Conformal Mapping and in Related Problems (erscheint 


demnächst). 
1°) Vortrag von L. Oollatz 1951 in Hannover, soll as veröffentlicht werden. 


17) Eine Feblerabschätzung von J. Grünsch für die 3. Randwertaufgabe soll demnächst ver- a 


öffentlicht werden, 
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Über die numerische Auflösung des dritten Randwertproblems 
bei elliptischen Differentialgleichungen 


Von Eduard Batschelet in Basel 
Vorgelegt sei die Differentialgleichung 


u u au ou 
ee) ey, reaylamıey). 


die in einem Gebiet @ der zy-Ebene mit dem Rand ( so aufgelöst werden soll, daß auf C eine 
Randbedingung dritter Art 


kl) 2 Als) u—/f(e) 


erfüllt ist. Da %(s) und h(s) auf Teilbögen des Randes auch verschwinden können, so ist der 
Fall von gemischten Randbedingungen eingeschlossen. Über den Rand (©, sowie über die Koeffi- 
zienten der Differentialgleichung und der Randbedingung werden Annahmen gemacht, welche 
die Existenz einer eindeutigen Lösung sichern, deren partielle Ableitungen genügend hoher 
Ordnung stetig sind. Solche Probleme ergeben sich in der Hydrodynamik, der Elastizität, sowie 
bei der Untersuchung der stationären Wärmeleitung und der Diffusion. 

Zur numerischen Bearbeitung solcher Aufgaben empfiehlt sich das Differenzenver- 
fahren, mit dessen Hilfe das kontinuierliche Problem durch ein diskontinuierliches ersetzt wird. 
Zu diesem Zweck überdeckt man @ mit einem (z. B. quadratischen) Punktgitter mit der Gitter- 
seite h. Die Randkurven werden in geeigneter Weise durch einen Gitterrand ersetzt, wobei der 
Approximation der Normalableitung von % durch einen Differenzenquotienten besondere Auf- 
merksamkeit geschenkt werden muß. Schwierigkeiten können auftreten, wenn der Rand Ecken 
besitzt. In.inneren Punkten werden die partiellen Differentialquotienten in bekannter Weise 
durch entsprechende Differenzenquotienten angenähert. Auf diese Weise erhält man an Stelle 
der gegebenen Differentialgleichung und der Randbedingung ein einfach gebautes System von 
linearen, algebraischen Gleichungen. Ihre Anzahl stimmt mit der Punktzahl des approximieren- 
den Gitters überein. 

Dieses Gleichungssystem wird am besten mit Hilfe eines Iterationsverfahrens näherungs- 
weise aufgelöst, das sich auf Blockrelaxationen stützt, wie dies vor allem von R. V. South- 
well vorgeschlagen wurde. Man kann zeigen, daß bei Randbedingungen dritter Art die Kon- 
vergenz des Iterationsverfahrens ebenso gesichert ist wie bei Randbedingungen erster Art. 

Von großer Bedeutung ist die Frage, wie gut die Gitterlösung die exakte Lösung der 
Differentialgleichung approximiert. Durch Modifikation eines Verfahrens von S. Gerschgorin 
lassen sich für den Fehler obere Schranken finden. Aus ihnen kann unmittelbar ersehen werden, 
daß bei nach null strebender Gitterseite h die Gitterlösung gegen die exakte Lösung des Problems 
konvergiert. Selbst bei einer Approximation der Differentialquotienten durch die einfachsten 
finiten Ausdrücke ist lineare Konvergenz gesichert. Ein besserer Konvergenzgrad wird erreicht, 
wenn die Approximation der Normalableitung am Rande durch eine Dreipunkteformel der 
numerischen Differentiation vorgenommen wird. 

Die gefundenen Schranken gestatten auch eine Diskussion der Abrundungsfehler, sowie der 
Fehler, die bei näherungsweiser Auflösung des vorliegenden Systems von linearen, algebraischen 
Gleichungen entstehen. 


Numerische Integration von Differentialgleichungen bei 
Approximation durch trisonometrische Ausdrücke 


Von W. Quade in Hannover 


Ist y=u(«) ein durch den Punkt a, b gehendes Integral der expliziten Differentialgleichung 
- erster Ordnung %y’=f(z,y), so besteht bei stetigem f(x, y) die Identität 


ua) —b+ [FE ul dE. 


"Bei den verschiedenen nach Adams benannten Verfahren zur numerischen Integration solcher 
- Differentialgleichungen wird der Integrand durch ein Polynom P (x) angenähert, das so bestimmt 
_ wird, daß die Differentialgleichung in isolierten Punkten erfüllt ist. 
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Das Problem der Annäherung des Integranden durch ein Polynom gehört der Approxi- 
mationstheorie!) an. Dort werden Näherungen durch Polynome und trigonometrische Ausdrücke 
untersucht, letztere bei periodischen Funktionen. Da jede Approximation durch Polynome auf 
eine solche durch trigonometrische Ausdrücke zurückgeführt werden kann, und der letztere Vor- 
gang sich in der Theorie als besonders fruchtbar erwiesen hat, war zu erwarten, daß er auch bei 
der numerischen Integration von Differentialgleichungen von Nutzen sein werde. 

Über die in dem Gebietsstreifen @: , 


aszsato, — oO <y<o, RS) 


erklärte rechte Seite von y'’=f(x,y) mögen die üblichen Voraussetzungen gemacht werden, 
welche die Existenz und Eindeutigkeit des durch den Punkt a, b gehenden Integrales sichern. 
Um Näherungslösungen mit Hilfe trigonometrischer Ausdrücke zu konstruieren, substituiere man 


RN) ’ 
ı=4-+0sin?’z; (StR. wi ar el 


t 
Dabei geht das Integral v(x) über in « (a + sin? - 3) =v(t), wobei v(t) dasjenige Integral von 


dy 1 


ROT ssinuf(a + asin 4 ») EEE 2) 


ist, welches der Anfangsbedingung v(0)=b genügt. Durch Fortsetzung in das Intervall (—,0) 
wird v(t) eine gerade Funktion mit der Periode 27. Die ebenfalls gerade periodische Funktion 


fa +asin? =; ve möge durch den trigonometrischen Ba NE, T,„(t) von der 
Ordnung n approximier t werden, der mit fin den äquidistanten Punkten t,— — ke 0,122.) 
übereinstimmt, M 
n—1 1 
RAR = ae 2 A,cosvt 5 A,cosnt; 


seine Konstanten sind 
= 1 1 n—1 } 
A,=— | -ht=(-1P fan t% fr cos 4), = N, 1 IE Me Pas 
n\2 2 k=1 


Über die Annäherung einer periodischen Funktion durch diesen Interpolationsausdruck 
gelten die folgenden Sätze: 


1FIstr (t) stetig und von beschränkter Schwankung, dann strebt 7‘, (£) mit zunehmendem n 
gleichmäßig im Periodenintervall gegen f{t). 

2. Genügt f(t) einer Lipschitzbedingung von der Ordnung e(0<o<1), dann strebt 7', (t) 
mit zunehmendem n gleichmäßig im Periodenintervall gegen FM: 

Bei der Approximation durch Polynome, wie sie bei den eingangs erwähnten Verfahren 
benutzt wird, stehen solche weittragenden Aussagen nicht zur Verfügung, da die Folge der 
Lagrangeschen Interpolationspolynome beiäquidistantenOrdinaten im allgemeinen divergentist. 


Wird 7',(t) so bestimmt, daß die Differentialgleichung (2) in den äquidistanten Punkten z, 
erfüllt ist, so folgen die unbekannten Ordinaten 9, ..., Y, der sich auf diese Weise ergebenden 
Näherungslösung aus dem System 


Y%—%Y>- (oft art us + Cm fa)» k=1, 2, A N... (3). 


Hierin bedeuten die c;, Konstanten, die sich bei festem n ein für allemal im voraus berechnen 
lassen, und es ist f;—=f(x; ,y,), wobei die x; sich gemäß (1) aus den ?, ergeben. 


Aus (3) können auf Grund der über f gemachten Voraussetzungen die n Unbekannten %, 


durch Iteration bestimmt werden, sofern nur & hinreichend klein ist; das System (3) besitzt genau 
eine Lösung. Mit zunehmendem n’ geht die Folge der die Lösung v(t) annähernden trigono- 
metrischen Ausdrücke gleichmäßig gegen die Lösung v(t). 

An einigen Beispielen wurde gezeigt, daß dieses Verfahren den bisher Henn: in bezug 
auf Genauigkeit überlegen ist; eine vorbereitende Rechnung ist bei ihm nicht erforderlich, 
sondern lediglich Iteration. Das Verfahren kann sowohl auf Systeme von expliziten Ditferential- 


gleichungen erster Ordnung als auch auf Differentialgleichungen höherer Ordnung ausgedehnt 


werden. Eine ausführliche Darstellung befindet sich in Vorbereitung. 


1) Vgl. C. de la Vallöe-Poussin: Lecons sur l’approximation des fonctions Fun variable röelle, 


Paris 1919; D. Jackson: The Theory of Approximation, New York 1930. 
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Fehlerabschätzungen bei Differenzenmethoden zur numerischen 
Integration Sewöhnlicher Differentialgleichungen 


Von @ermund Dahlguist in Stockholm 


Bei einer Arbeit mit der,schwedischen Relaismaschine BARK haben wir Fehlerabschätzungen 
für eine Reihe von Integrationsmethoden gemacht. Die Ergebnisse und Methoden werden hier 
der Kürze halber nur an einem besonders einfachen Beispiel angedeutet. Wesentlich ist, daß 
das Integrationsverfahren eine Ditfferenzengleichung liefert, die von höherer Ordnung ist als die 
zu integrierende Differentialgleichung. 

Wir wollen das Integral der Gleichung 


Y 
179% 
berechnen und gebrauchen dazu die Formel 


. Yn+ı = In-ı + 2h fo bye 

y, ist die exakte Lösung, %, ist die Näherungslösung. Man hat 

Ynt = Yn-ı + 2h Fly d+ 6; 
wo der Quadraturfehler ö, durch die Formel 

BES Schr 

6,= ER Y(r) (7 ist ein Mittelwert) 
gegeben ist. Setzen wir 
Yn Y len 


so erhalten wir nach Vernachlässigung von Gliedern höherer Ordnung in e,, 


of 
En+i = &n—1 -H 2 h n ir) ” Ey -H Ö, . . . . . . . . . . (ii 

/n { 
(Die vernachlässigten Glieder können übrigens auch durch eine Korrektur von 6, repräsentiert 
werden.) Die Differenzengleichung (1) kann etwa in folgender Weise angenähert gelöst werden. 
Erstens hat sie eine partikuläre Lösung &,, die sich nur wenig von der Lösung &(t) der linearen 
Differentialgleichung 


BE Aa 2 , 
Arie: DR e-+öf(t); N FE ES Ä 
unterscheidet, wo öft) = 6,/2h = h?y (1/6. iR Folge u, = e,— e, genügt der homogenen „2 
Gleichung zweiter Ordnung di 
B) 2 
Un+1ı = Umn-1 nn 2h Bi ERRENHEE VÜRRO I RS ER RE (3), E 
die zwei linear unabhängige Lösungen hat. Eine dieser Lösungen unterscheidet sich nur wenig % 
von einer Lösung v der Differentialgleichung 
dv of 
#2 — 2y mM NET TE REN ee . . . . . (4), 


und man findet leicht nach der Substitution v,—(— 1)" w,, daß die aus (3) entstandene Gleichung 
für w, eine Lösung hat, die sich nur wenig von einer Lösung der Gleichung 
dw of 
unterscheidet. Wir erhalten somit 
&, zEl(t,) + A v(t,)+B- Dr 1)" w(t,).- Pr 
Die Koeffizienten A und B können aus den Anfangsbedingungen bestimmt werden. Bezüglich 
des ersten Gliedes ist es bemerkenswert, daß ö(t) im wesentlichen von der dritten Ableitung 


von y bestimmt ist, und daß ö(£) darum gewissermaßen inResonanz mit derLösung der homogenen 
Gleichung steht. Der Fehler hat auch einen alternierenden Teil, d.h. (—1)*w,, der von Be- 


deutung sein kann, wenn der Realteil von 9f/dy negativ ist. Für die Gleichung dy/dt = — y 
erhält man z.B. F: 
= h? 
Mil FEN FYo HT + ds are R 


» 
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was man auch durch explizite Lösung von der Differentialgleichung und der Differenzengleichung 
finden kann. Für h=0,1 ist das letzte Glied das größte, wenn £>2. Sein Betrag wächst expo- 
nentiell. 

Die Zerlegung des Fehlers auf der rechten Seite von (6) ist für mehrere andere Integrations- 
methoden möglich. In (3) und (4) wird nur der Ausdruck für ö(f) verändert. Gl. (5) lautet dann 


dw Al 
— W, 
dr 0" ay 
wo ceine für das Verfahren charakteristische Konstante ist. Z.B. ist c = — 1/3 für die Verfahren, 


wo eine Extrapolation durch mehrfache Iterationen mit Hilfe der Simpsonschen Regel ver- 
bessert wird. Wenn man die Simpsonsche Regel in jedem Schritt nur einmal zur Verbesserung 
benutzt, so ist e von der benutzten Extrapolationsformel abhängig. 

Ich beabsichtige, die Herleitungen und Resultate später vollständig zu veröffentlichen. 


Ein Summenverfahren für Rand- und Eigenwertaufgaben 
linearer Differentialgleichungen!) 


Von H. Sassenfeld in Darmstadt 


Das Summenverfahren dient zur Lösung von Rand- und Eigenwerlaufgaben und stellt 
eine Verbesserung des bekannten Differenzenverfahrens dar. Es macht von Integrationen Ge- 
brauch und legt Quadratur- statt Differenzenformeln zugrunde. In einfachen Sonderfällen führt 
das Verfahren auf das von L. Collatz!) eingeführte Mehrstellenverfahren. Die Anwendung des 
Summenverfahrens ist genau so einfach wie das Differenzenverfahren und führt auch bei kompli- 
zierten Differentialgleichungen zu übersichtlichen Ansätzen. Im folgenden wird das Verfahren 
nur für gewöhnliche Differentialgleichungen hergeleitet, es ist jedoch as weiteres auf partielle 
Differentialgleichungen anwendbar. 


Beim Differenzenverfahren werden in dem linearen Differentialausdruck 


N \ 
SR Ne ER A N AN 
y= r 
die Differentialquolienlten 4) approximiert durch 


ye Zn +R=By HR VI RR OT (2) 
mit den Abkürzungen %; = y(x, + th) A R = Restglied. 


Die Ausdrücke (2) werden (ohne Restglied) für alle im Integrationsbereich liegenden 
Punkte 2, x,—+thin (1) eingeführt und die Lösung des entstehenden Gleichungssystems liefert 
die gesuchten Y;- 


Beim Summenverfahren integriert man zunächst die Differentialgleichung (1) n-mal. Setzt 
man die a,(x) v-mal differenzierbar voraus (bei selbstadjungierter Form immer erfüllt), s so wird 
durch wiederholte Anwendung von Teilintegrationen aus (1): 

Mo) m]... [Ldn— If. [b,(a)ydartmf-.- [gdr + Pr ı@—m)=0 (8). 
& To v=0 & % %o . 


Hier werden jetzt die Integrale durch Quadraturformeln ersetzt, von der Art: 


oft [fasr =eßul +R=8i+® nr N (4). 


* Das in (3) noch auftretende Ba le Grades P„_1(2—x,), das die bei Randwert- 


aufgaben z. T. unbekannten Anfangswerte enthält, eliminiert man durch Differenzenbildung: 


N 


= () (IH M@)=d,N= > 9, 95,, #3, 5,9 =0; Time by Zee 
v0 


i= 


wobei die ;s, für die nenne zur Matrix ©, zusammengefaßt sind. Die d„©,, die die Kompo- 


nenten von d,6,f= - nyW) J. sind, können ein für allemal unabhängig von der speziellen Form 


der a,(z) in (1) ausgerechnet werden. . 


466 8. 


> 


Fa 


1) L. Collatz: Eigenwertaufgaben Br technischen Anwendungen. Akad. De, "Leipzig 1949, NA 2 


E42 


TER an 


N 
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Die Übersetzung von Randbedingungen mit y® am Rande geschieht durch Linearkom- 
. bination der Gleichung (3) in Form von ‚‚links“- oder „rechtsseitigen‘ Differenzenformeln dr: 


k n 
(kV) - 
BEL ER 2 RO re A 
v—=0 v0 
Hiermit sind die Ableitungen von y bis zur k-ten Ordnung an der Stelle &— x, durch Linear- 
kombinationen der %,(z), d.h. wegen u,—b,% durch Linearkombination der y selbst dargestellt. 
Es lassen sich zahlreiche Verschärfungen des Summenverfahrens angeben. So kann man 
Hermitesche Interpolation anwenden. Auch kann die Differentialgleichung L[y] mehr als 
n-mal integriert werden, womit man ebenfalls zu genaueren Ergebnissen kommt. Abschließend 
sei noch erwähnt, daß für a, —= const und «a, bis «,_ı identisch null das besprochene Summen- 
verfahren mit dem Mehrstellenverfahren identisch wird ?). 


Zur Fehlerabschätzung für das Ritzsche Verfahren 
bei Eigenwertaufgaben 
Von @. Bertram in Hannover 


Das Ritzsche Verfahren liefert bei der Anwendung auf Eigenwertprobleme unter ge- 
eigneten Voraussetzungen obere Schranken für die Eigenwerte. Es wird eine Methode zur 
’ gleichzeitigen Gewinnung unterer Schranken aus dem Ritzschen Verfahren angegeben. Be- 

trachtet werden gewöhnliche, selbstadjungierte, volldefinite Eigenwertaufgaben der Form 


M[y(@)] = AN [y(e)] 
mit den Randbedingungen 


U,[y]=0 Mazalıen 2m), 
wobei M und N lineare homogene Differentialoperatoren der Gestalt 
m NR 
Miyl= 2 1) [na y]” ; Ny=2 1?) y”]” (m>n) 
= y= 


sind. Die f, und g, unterliegen geeigneten Definitheits- und Differenzierbarkeitsvoraus- 
setzungen. Die Randbedingungen sollen linear homogen in den Randwerten und Randablei- 
tungen bis zur (2m —1)-teu Ordnung, voneinander linear unabhängig und frei von A sein. Als 
Ansatzfunktionen werden die Eigenfunktionen eines ähnlich gebauten, definiten Problems 


Pliyl=Ay:; U.yJ=0 
H 
mit denselben Randbedingungen und P[y] = 3% (—1)’ [p,(z) y"]® verwendet. Es werden Be- 
v0 


‚ziehungen des Ritzschen Verfahrens in der Galerkinschen Fassung zu Integrodifferential- 
> gleichungen aufgezeigt. Weiterhin wird der Spezialfall der Eingliedklasse, d.h. 


Ny] = — 1’) y]” 


diskutiert. Hier gelingt unter geeigneten einschränkenden Voraussetzungen die Zurück- 
führung der Integrodifferentialgleichungen auf Integralgleichungen. Da sich sowohl für die 
Eigenfunktionen als auch für die Ritzschen Näherungen ähnlich gebaute Integralgleichungen 
angeben lassen, Kann man durch Einführung eines Störungsparameters & eine Gleichungsschar 
angeben, die speziell für «=1 die exakte, für «—=0 die Näherungslösung ergibt. Die Betrachtung 
dieser verallgemeinerten Integralgleichung erlaubt im Anschluß an eine Idee von N.Kryloff die 
Aufstellung einer Alternativgleichung, aus der sich Abschätzungen für den Fehler A,—), der 
genäherten und der exakten »v-ten Eigenwerte gewinnen lassen. Nebenbei ergeben sich asympto- 
tische Formeln für das Konvergenzverhalten der Eigenwertfehler in Abhängigkeit von der 
Ansatzgliederzahl, die einfache Konvergenzaussagen gestatten. Für eine Reihe technisch wich- 
tiger Problemklassen lassen sich von Fall zu Fall die Abschätzungsformeln noch wesentlich ver- 
bessern. Die aus dem Verfahren gewonnenen Ergebnisse haben sich als praktisch brauchbar 
erwiesen. 
Einzelheiten soll eine in Kürze erscheinende Arbeit enthalten. 


2) H. Sassenfeld: Über lineare Gleichungssysteme und ihre Anwendung bei Summenverfahren für 
Rand- und Eigenwertaufgaben. Dissertation Technische Hochschule Darmstadt 1951, 99 S. 


u 


DRG ETN, 
. A 


” 


var 
Er 


> 


er 


er N u 


TE DI 


949 A. Praktis che Analysis 


e Z.ang 
Bd.31 Nr. 8/9 Aug./Sept. 1951 


Ein Analogon zu Grammels Methode der graphischen Integration 
gewöhnlicher Differentialgleichungen 


Von Walter Gautscht in Basel 


Bedient man sich an Stelle rechtwinkliger Koordinaten der Polarkoordinaten, so gelingt 
es, Ditferentialgleichungen beliebig hoher Ordnung in einfacher Weise graphisch-rechnerisch zu 
integrieren, falls Anfangswerte gegeben sind. Eine Möglichkeit in dieser Richtung hat R. Gram- 
mel!) durch sein Orthopolarenverfahren aufgezeigt. Es soll hier ein entsprechendes Verfah- 
ren mitgeteilt werden, das vermeidet, mit reziproken Funktionswerten zu arbeiten. 

Sei u(z) eine genügend oft differenzierbare Funktion. x und % (für negative u: +7 
und |«]) als Polarkoordinaten aufgetragen, ergeben in einem gewissen Variationsbereich von & 
eine Kurve (, das sog. Polarbild der Funktion v(z). In entsprechender Weise mögen die Ortho- 
polarbilder 0’, C",.... aus w (x), u" (x), ... hervorgehen, wobei jeweils die Nullage der Winkel- 
messung um 7/2 zu drehen ist (Bild 1). Sind P, P', P", ... entsprechende Punkte der Kurven 
0,0',C",..., so steht die Tangentialrichtung des Polarbildes C im Punkte P senkrecht auf 
P' P, diejenige von 0” in P’ senkrecht auf P" P’ usw. 


Bild 1 Bild 2 


Diese Eigenschaft wird zur Lösung von Differentialgleichungen ausgenützt. Beispielsweise 
sei die Gleichung zweiter Ordnung 


Ye EU) RS Le BER Ed 
mit den Anfangsbedingungen 
‚ yo) (&,) — yo») (v = 0, 1) 1 RE TE (2) j 


vorgelegt. Man wird zunächst die Anfangspunkte P® vermöge den Anfangswerten (2) markieren 


‚Bild.2). Zur Bestimmung von P5 dient (1). In den Punkten P® hat man senkrecht zu den 


Verbindungslinien PP+t P®) (die nicht notwendig zu zeichnen sind) ‘die Fortschreitungs- 
richtungen. Diesen wird man folgen etwa bis zum Schnitt mit dem um einen genügend kleinen 
Winkel Ax gedrehten Polarenkreuz. Sind P®) die Schnittpunkte, so gilt näherungsweise 
OPW— yN(xz,+4Ax). Pi berechnet man, unter Verwendung dieser Näherungswerte, wie- 
derum aus (1). Die nächsten Punkte P, gewinnt man aus den P, so, wie man die P, aus den 
P, konstruiert hat. Nach einigen Schritten kommen genäherte Polar- und Orthopolarbilder 
der gesuchten Lösungsfunktion y(x) heraus. 

Es besteht keine Schwierigkeit, die Verfeinerung der Konstruktion, wie sie Grammel 
seinem Verfahren beigelegt hat, den neuen Gegebenheiten anzupassen: Nachdem man, ausgehend 
etwa von den Punkten P, die Punkte P, konstruiert hat, ermittle man die Mittelpunkte M,, 
Mo, Mo der Strecken P,P,, PoPı, Po Pi. Schreitet man jetzt senkrecht zu den Verbindungs- 
linien MP+D MW weiter, so gelangt man zu genaueren Punkten P®), von denen man zeigen 
kann, daß sie die Funktionswerte 4) bis zur zweiten Ordnung genau wiedergeben. Pj kann 
nötigenfalls aus den P®) genauer berechnet werden. Die zuerst konstruierten Fortschreitungs- 
richtungen PP) P®) bleiben als Tangenten in den P®) weiterhin nützlich. | 


1) R. Grammel: Ing.-Arch. 10 (1939), 8. 395. 


a vs Br 
ew. Math. Mech. 


el ra 1951 - A. Praktische Analysis r - 243 


Sollte eine Nullstelle der Lösungsfunktion oder ihrer Ableitung (radiale Fortschreitungs- 
richtung!) die Konstruktion aufhalten, so berechne man sich die nächsten Punkte durch Taylor- 
entwicklung. Die nötigen Ableitungen stehen, teils in der Konstruktion, teils in der Differential- 
gleichung zur Verfügung. Unter Umständen bringt auch die Einführung einer neuen geeigneten 
Funktion Vorteile?). 


Bemerkungen zu den Iterationsverfahren für lineare Gleichungsysteme 
Von U. Wegner in Heidelberg 


Ist Ar=;3 ein lineares Gleichungssystem mit quadratischer Matrix A, deren Elemente 
Q;;(t, k=1,2,...n) heißen mögen, so ist das allgemeinste lineare Iterationsverfahren bestimmt 
durch den folgenden Ansatz für die (m-+1)-te Approximation tn+1 der Lösung r: 


Im+i RR u TER P„(4 Em uT, 3) TER the Pre Pi (I), 
Em+1 = Om+1 Em+ı + Am+ı m—=0,l,... :.. (b) 


wobei t, willkürlich gewählt werden darf, und P,, eine Matrix bezeichnet, die von %,, und a;; 
abhängig sein kann. Qn+ı und ay+ı Können ebenfalls von den a;, abhängig sein und mit 
dem Index m variieren. (I,) können wir schreiben: Imıı= Cyim + Pnz mit („= E— PA. 
In der verschiedenen Wahl der P,, und Q,, unterscheiden sich die einzelnen Iterationsverfahren. 
Ist z.B. A=R-+S mit existentem R-!, so erhält man in P„=R, Qm+sı=E, a„=0 das 
Hertwigsche Verfahren, was beispielsweise konvergiert (wichtiger Fall für die Baustatik), wenn 
A=R-+S und R—S positiv definit sind. — Ist P„=4'-U mit U=(uw;,), uwr=0; i+k, 
2M, ® 
Um m . 3; 9, = i»und willkürlichem M;>0, so erhält man bei — u 
TS UHML+. IMs 5; Ge un illkürlichem M,>0, so erh n bei Qm+ı1 
Am =0, ein stets konvergentes — jedoch nur für lineare Integralgleichungen praktisch ver- 
wendbares — Verfahren (Cimmino). — Ist P,=P,, falls s=v-9-+r mit 0<r<g ist mit 
ganzzahligem v® und g, so erhält man wieder mit Q,=E und a„=0 ein konvergentes Ver- 
fahren, falls Adi .< ,i=1,2,...9 OG=BE—P;-A, APi= charakteristische Wurzel von 
0,0, (= Norm von (;) ist. Speziell P,=t;- E mit komplexer Zahl t, liefert die Konvergenz- 
bedingung: |1—1t;-44|<1; 44 = charakteristische Wurzel von A (Bückner). — Setzt man 
P„= als n-reihige, quadratische Matrix an, so erhält man mit Q,=E, u„—=0— also T„—= Im — 
ein sehr brauchbares Verfahren zur Approximiation von A’! bei der Konvergenzbedingung: 
A#-&4| <1 (Günther-Schulz). — Setzt man P„=D„= (air), am —_0 für i+%k und 
akt f. mit d;,im> Wobei fm folgendermaßen bestimmt ist: 1. Entweder: 


n 
v,., (Em) Max lo; (x„)| mit v%;(ru) = — | e3 a) — x). 2.. oder: 9%, (£m)=%,,(Em) mit 
=1,23..'.n = 
Y„=m (mod n) undO<r„<n (1. und 2. Gauß-Seidelsches Verfahren). Ist A positiv definit, 
oder [44] <1 und a,;-#0, so tritt Konvergenz ein. Auch hier ist Q,—=E und a„—=0. — Eine 
Erweiterung liefert das folgende Verfahren: Ist A nicht-negativ definit, @;;-#0, so bestimmt die 
obige Folge r,„, den Grenzwert: lim Ar„ = 3. Bei hochgradig statisch unbestimmten Systemen 
(z.B. Tragflügel) findet das Verfahren in der folgenden Form Verwendung: A = positiv 
definite Matrix a, b, c= willkürliche Vektoren. Bildung des Gleichungssystems: 
' Ay+y9Aa +:04-5 +MAcr =;, 
aAy+ yPa’Aa+zMa’Ab +tWa’Ac=u';z, 
0’Ay-+ yOb’Aa + 20h’Ab +9’ Ac—b}5, 
eAy-+ yMec’Aa +2 d Ab + Ar = cz f 
mit den n-H3 Unbekannten 9, yO, 20, {0 Anwendung des Gauß-Seidelschen Verfahrens 
liefert: lim (y,+ at b+in c)—= A}. (Die Erweiterungen der Matrix A werden bei- 
Mm © p 


spielsweise durch die statischen Gleichgewichtsbedingungen gegeben, die zwar in Ar—=3 stecken, 
aber explizite nicht erkennbar sind. Der obige Ansatz zwingt zu schnellerer Konvergenz.) — 


2) Durch eine freundliche Zuschrift von Prof. K. Federhofer (Graz) bin ich auf die Arbeiten von 
R. Neuendorft in der Z. angew. Math. Mech. 2 und 3 aufmerksam geworden. Darin ist die Idee einer solchen 
Konstruktion bereits ausgesprochen. Als Fortschreitungselemente werden Krümmungskreisbögen und zur 
Verbesserung der Näherungslösung eine graphische Auswertung der Methode der sukzessiven Approximation 
angewendet. 
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Nach Gauß-Seidel wäre z.B. bei existentem A-! das Verfahren tn+1 = Im — Dy (4’A A) 
' stets konvergent, praktisch aber unbrauchbar. ‘Setzt man dagegen an: 
Em+1 Ta Ir TB Di (A’A Im sr 4'3) D 
Im+ı = Om+ı Imtı + Am+ı 
me On, nr dr 3, SO erhält man durch diese Transformation ein stets kon- 


vergentes Verfahren mit den zwei Variabilitäten, bedingt durch die Formen für D,,, das ein- 
BR fach und schnell arbeitet und zusammengezogen werden kann in die einfache Form: 


Lm+1 = Em AD,„«4A' In 3) . 
(Die geometrische Interpretation ist sehr einfach.) 


Praktisch rechnet man so: 
! a; n € 
=" +0, 5 Um(Em)s vahy..m mer; Sir. 
% (Man beachte besonders, daß v,,, (%,.) vom Index » unabhängig ist, was das Verfahren so brauchbar 
Rz macht) (Kaczmarz). Ist ,= (E—4’D,4) (E—4’D.14)--- (E— A'D,A) (k=1,2,...,n); 
ds $,—kte Spalte von A, so kann für jede der Variablen x; eine independente Darstellung gegeben 
werden in der Form: 


Er IB I,)23; (Eee 


1 R (E—II,)* wird dabei durch den Iterationsprozeß geliefert. Der Formalismus der Gauß- Ü) 
Seidelschen Verfahren und des Kaczmarzschen Verfahrens eignet sich besonders — in 

$ der Auffassung von A als linearer Operator — zur Auflösung der Fredholmschen Integral- 

ra gleichungen 1. und 2. Art. — Durch Transformationen im Gauß-Seidelschen Verfahren 

Be der Gestalt: 


Em+1 = Im — Dm (A PA’ Gi (4), mn Am — Pit 41; 
oder zusammengezogen in 
Im+1 = im — A'GD„@’APAM(Arm — 3) 


mit beliebiger nicht-singulärer Matrix @ und positivem definiten P, erhält man viele stets kon- 


vergente — und darunter alle bis jetzt bekannten — Verfahren (z.B. P=E; @beliebig liefert 


[6 
As in 


das Hardy-Cross’sche Verfahren). — Die Konvergenzgüte hängt vom Quotienten im 


Bine 
allgemeinen ab. Je Br dieser Quotient ist, um so besser ist die Konvergenz. Daher sind Ab- 


schätzungen von Ale von unten und für: I von oben her notwendig. f 
Zu den bekannten Verfahren seien noch einige a die für die Praxis von großer Be- 


deutung sind, angeführt. Es sei s= Spur der Norm von (;; d = Determinante der Norm von 0, 
Dann ist: ar 
1. Au>s— Ye nam ii 2d. Bez 3 IK \ 
A N n—1 d r 
9, Es sei mp =s— (n—1) Sr Dann ist: ee 2 uZ Ad, wobei 
m - y 
u—E"t ist, und &= größte positive Wurzel von = s: 2 — (n— 1) ya ist. 


61101 


Qı C, 


Gr = Gı + 0 (e- E— 0) As. @=1 Be) 
I = Ct 01 N C)"a, | 


mit g, beginnend, wobei g, aus Ayraz > gV) > AQjur gewählt wird. Ist 208 einfach, ‚so gilt: » 
AHIR<p<t; G, ART, ah > BE, 
mit lim g, = lim 6, — Afyaz und x RE Er 
v© vo» AN 
%: = (9: E —Oy’a>%; 9 (GE — 003 u 0 En 
so daß (5 ) Eigenvektor von O zu Ayyax wird. (94. E08 wird natürlich, Br S 


3. C sei symmetrisch; (= ( ) dann bilde man: 


u 
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4. Ist C=(6,,); %,>0; so gilt nach Perron-Frobenius: 


Om+1 . Omtl 
Om N Is Hr >0 3 Nm — Min 2  Alfar 


M,„ = Max 
Cm vs 


mit beliebigem x, > 0. Aber es giltsogar: 
M,=zM,2M;,>--- ee >. ZMSMSM 


was für die praktische Berechnung von größter Wichtigkeit ist. 


n ’ 
1% Gir — Am ik 


616; 


; 2, ==: "Min 


Brust 0, 
i=2, 8,40. ‚n el CGıı Am 


so konvergiert x„—A,, gegen die Wurzel, die in der Nähe von c,, liegt, falls x,+c,, in der Nähe 
von A gewählt wird, und wenn diese charakteristische Wurzel von C’ einfach ist. 
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Über das Iterationsverfahren 
Von Johannes Weissinger in Hamburg 


Die Anwendungen des Iterationsverfahrens auf die verschiedensten. linearen und nicht- 
linearen Probleme lassen sich als Spezialfälle eines abstrakten Satzes auffassen, der zunächst 
formuliert werden soll. 

Gegeben seien ein vollständiger, metrischer Raum!) 2 mit Elementen &, n, ..., deren 
Abstand |£&. | also die Dreiecksungleichung |&, n|<]&,&]-+I&, n| befriedigen soll, und eine Halb- 
gruppe?) HZ von Abbildungen 7 des Raumes Q in sich. Wir sagen, in H sei ein Betrag erklärt, 
wenn jedem T eine reelle Zahl |7|> 0 zugeordnet ist, so daß gilt 


BEE Ernie für beliebige En) 44. 0 ene nee 20h 

Gewöhnlich ist Z die von einem festen 7’ erzeugte Halbgruppe H = {T"}. Dann ist durch 
|7*r|=|7|® ein Betrag in H definiert, falls (1) für die Erzeugende T gilt. 

Hauptsatz: Die Folge &, = T*r&,—= T&„_ı Konvergiert bei beliebigem &, sicher dann 


und zwar gegen die eindeutig bestimmte Lösung & der Gleichung &= TE£, wenn & |T”| konver- 
n=1 
giert, also sicher wenn |7|<1 ist. Im letzteren Fall gilt die Fehlerabschätzung 


IT| 
[En &| = 1—|7] | En nl. 
Der sehr einfache Beweis soll hier unterdrückt werden, um statt dessen einige (lineare, 
reelle) Beispiele zu betrachten. 


Satz 13): Zerlegt man die symmetrische, positiv definite Matrix %: 
Y=-D-G.-% B=8+5F (8 Transponierte von 8), 


1) Für die meisten Anwendungen genügt es, 2 als additive Gruppe vorauszusetzen. Die Verallgemeine- 
rung zu einem metrischen Raum verdanke ich Herrn W. Maak. 

2) D.h. mit 7,, 7, soll auch 7,7, in H liegen. 

3) Dieser Satz enthält als Spezialfall ein bekanntes Kriterium für das Einzelschrittverfahren [von 


Mises und Pollaczek-Geiringer: Z. angew. Math. Mech. 9 (1929). Hinsichtlich weiterer Literaturhin- 


weise muß auf eine demnächst erscheinende, ausführlichere Darstellung verwiesen werden. 
17 


rn 
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so daß (D—B,)! existiert und D-+-% symmetrisch und positiv definit ist, so Konvergiert die 
Vektorfolge 
In+1 = Cr. +t, = (D 2x2 3) B, 
bei beliebigem x, gegen die Lösung von r—= Cr —+t (und damit von ı =r= (8 — B)t). 
Beweis: 2 sei der Raum der Vektoren x mit der durch X bestimmten Metrik 


ee  IeHl=le-9l. 
T ist die Transformation 6x -+t, H= {T"}. Dann ist in Z durch |7|=|E|, |7*| = zn 
ein Betrag definiert, wenn |&x| < |E| |x} für alle x gilt, und der Satz ist bewiesen, wenn |E] < 
ist. — Nun folgt aus der leicht verifizierbaren Identität 
EAIE=-A—%, P-(E-E)(D-H+FT)lE— OB), & —= Einheitsmatrix 
durch Multiplikation mit x’ bzw. r 
Sr=k—rPe abo JEr<[le mitt JE=YI—p<I1, 

da ® positiv definit und daher !PBr>ple, 1>p>0 ist. 

Der Satz läßt sich mühelos auf semidefinite Matrizen verallgemeinern und auf Fredholm- 
sche Integralgleichungen zweiter Art übertragen, so daß er als Grundlage einer allgemeinen 
Lösungstheorie für solche Integralgleichungen dienen kann. 

Ist jeder m-reihigen quadratischen Matrix V = (a;,) ein „Betrag“ A|>0 zugeordnet 
mit folgenden Eigenschaften e 

a) |Y]=0 ist äquivalent mit Y=0, Q 

b) aus |YW,|>0 folgt ,>0, 

JH SIH HUN .. ..r° 

- für beliebige W,, Xs, 

LAS 
so gilt 

Satz2: Ist |8|<1, so konvergiert die Matrizenfolge &„4ı = BX, + © bei beliebigem %, 
gegen die Lösung & der Matrizengleichung £&=DBXE- ©, und es ist 


f B 
| & —H| = “ Er | &, &n—ı] 
Be -. Der Beweis ergibt sich sofort aus nn Hauptsatz, wenn man für 2 den Matrizenring 
F nimmt und H mit 2 identifiziert. Man bestätigt nun leicht, daß z.B. durch |Y]| = nz = le; :| 


A oder |Y| = uk 2| a; oder |X| = Y21a | a,,|?ein „‚Betrag“ im Matrizenring definiert ist. Jede 


Ey. dieser Beizesdefisitioneh liefert, in Satz PD eingesetzt, ein bekanntes Konvergenzkriterium nebst 
| Fehlerabschätzung für das Gesamtschrittverfahren (wenn ©, &, und damit alle &, als Vektoren 
deutbare ‚‚Spaltenmatrizen‘“ sind). Natürlich kann man Satz 2 auch ohne Berufung auf den 
abstrakten Hauptsatz unmittelbar aus den Eigenschaften a) bis d) herleiten, die ja besagen, 


; daß man mit Matrizen und ihren Beträgen im wesentlichen genau so rechnen darf wie mit Zahlen. 
T Man hat damit eine sehr einfache und elementare Methode zur einheitlichen Formulierung und 
a Herleitung der üblichen Kriterien. 

3 Lagrange-Hermitesche Interpolation im Komplexen 


Von H. Unger in Darmstadt 


Die Interpolation im Komplexen ist eine immer wichtiger werdende Teilaufgabe der Ver- 
tafelung spezieller Funktionen für komplexes Argument. Solche Tafeln erfordern doppelten 
Eingang und können im allgemeinen nicht so engmaschig tabelliert werden, daß eine Inter- 
polation mit einem Polynom niederen Grades durchführbar ist. Um mit möglichst wenig Stütz- 
werten auszukommen, muß eine hochgradige Interpolation herangezogen werden. Für diese 
Zwecke erscheint die von Hermite gebrachte Erweiterung der gewöhnlichen Eagraußsschen 
Interpolationsformel unter Heranziehung der ER 


P@)= > > L.2(@) gie | j 


% besonders geeignet. Legt man die Stützetellen 1 in die Eckpunkte eines regelmäßigen Polyebae? 
> dann kann man die Koeffizienten Z,,,(z) auf ua der Beziehung 


Fe pn 
L,,»(2) = e "In „lei] 


a0 235 Pi » rs rn 
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auf die Koeffizienten mit dem Index 1 an den abgeänderten Argumentstellen 
27 

Er 3 m tm) 

zurückführen. 

Für die praktische Durchführung bedeutet dies im Falle eines Quadrates, dem man im 
allgemeinen den Vorzug geben wird, daß man für eine Interpolation mit einem Polynom vom 
beispielsweise 11. Grade an Stelle der 12 Koeffizienten nur drei zu vertafeln braucht, nämlich Lo» 
L,,, und L,.. Dabei müssen in den Eckpunkten des Quadrates die Funktionswerte, die ersten 
und zweiten Ableitungen der zu interpolierenden Funktion bekannt sein. Bei vielen speziellen 
Funktionen können die Ableitungen einfach berechnet werden (z. B. durch Rekursionsformeln, 
Differentialgleichungen). 

Eine ausführliche Darstellung erscheint in der Zeitschrift für angewandte Mathematik 
und Physik, Zürich. 
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Eine Bemerkung zum Ausgleichen fehlerhafter Meßreihen 
Von K. Stange in Karlsruhe 


Sind die durch Messung bestimmten Koordinaten (x, %, ...) einer Punktreihe, welche der 
Gleichung F(z, y,...; a, b,...)= 0 genügen soll, mit Meßfehlern behaftet, so bestimmt man 
die Parameter a, b,... nach der ‚‚Methode der kleinsten Quadrate‘ so, daß gewisse ‚‚verbesserte 
Meßwerte“ x, 9, ... der vorgegebenen Gleichung #—=0 genügen. Dabei schiebt man die Meß- 
fehler oft mit gewisser Willkür einseitig einer der beiden Veränderlichen zu. Da das Ergebnis 
des Ausgleichens jedoch von den getroffenen Annahmen über die Verteilung der Meßfehler 
abhängt, so fälscht man mit der eben genannten vereinfachenden Annahme unter Umständen 
das Ergebnis von vornherein in nicht übersehbarer Weise. Ein sinnvolles Ausgleichen ist nur 
möglich, wenn man die dem Meßverfahren zugeordnete Streumatrix kennt. 

Im allgemeinen Falle werden die Parameter a, b,... unter Verzicht auf einschränkende 
Voraussetzungen so bestimmt, daß die Wahrscheinlichkeitsdichte an den Meßpunkten ? in 
bezug auf die ausgeglichenen Punkte P möglichst groß wird. Der Weg führt auf ein System 
nichtlinearer ‚‚Normalgleichungen“. Die Lösung dieses verwickelten Systems läßt sich jedoch 
iterativ schrittweise gewinnen, wobei man bei jedem Einzelschritt nur ein System linearer 
Gleichungen zu lösen hat. Die unter vereinfachenden Annahmen geltenden bekannten Lösungen 
der Ausgleichaufgabe lassen sich als Sonderfälle in den allgemeineren Rahmen einordnen. 


Fehlerabschätzungen zu den Verfahren von Graeffe 
und Brodetsky-Smeal 


Von E. Weinel in Jena 


Durch einen unmittelbar an das Graeffesche Verfahren anschließenden rekursiven 
Rechenprozeß lassen sich nach dem Vorgang von Brodetsky und Smealt) die komplexen 
Nullstellen eines reellen oder komplexen Polynoms vollständig bestimmen. Für den Fehler, der 
sich beim Abbrechen des Verfahrens nach dem r-ten Schritt einstellt, lassen sich einfache Ab- 
schätzungen angeben, die unmittelbar aus den im Zuge der Rechnung entstehenden Zahlenwerten 


gebildet werden können. 
Ein ausführlicher Bericht folgt in einem späteren Heft dieser Zeitschrift. 


Bestimmung der Sonnenfleckenperioden mit Hilfe des Periodogramms 
= Von H. Wundt in Freiburg i. Br. 


Eine empirische, in einem Intervall der Länge 7 tabellarisch gegebene Funktion f(t), deren 
Zusammensetzung aus periodischen Funktionen unbekannter Periodenlänge uns interessiert, 


können wir nach Fourier analysieren. Die Gesamtamplituden h(A) — Ja?(A)-+b?(A) bilden, 
über der Wellenlänge A als unabhängiger Veränderlicher aufgetragen, ein diskretes Linien- 


spektrum. Führen wir statt A die Kreisfrequenz = 7 ein, so werden die „Fourierabszissen‘‘ 


1) Brodetsky und Smeal: Proc. Cambridge philos. Soc. 22 (1924). 
17* 
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= un äquidistant. Die zugehörigen ‚„‚Fourierordinaten‘“ sinnvoll zu interpolieren ist 


Aufgabe der Periodogrammrechnung. 
Eine spezielle Methode, eine solche interpolierende Funktion IJ/(w) zu finden, ist die folgende 
(vorgeschlagen von Ostrowski): 


(t) wird mit der Einzelschwingung Ah - cos (wo? — 9) verglichen und das mittlere Abwei- 
chungsquadrat im Intervall 7 durch Variation der drei verfügbaren Parameter h, o und @ 
minimal gemacht: 


Tre —R- cos (wt— p)]2dt= Min. aa: Een el, 
J 3 
Differenziert man zunächst nach A, so findet man, daß der Ausdruck 
. i ° 
ro: cos (wi — p) dt 
0 - a 


7 
[ cos? (ot — p) dt 
0 


f maximal gemacht werden muß. Sein Nenner N ist vom speziellen Problem unabhängig. Es 
X bedeutet daher eine Rechenvereinfachung, wenn wir ihn als konstant, also von und von @ 
er unabhängig, betrachten könnten. Dies ist der Fall für o7T>1, nämlich 
7 
Iim.N Sa ee a 3). 
und | | ) 


Es bedeutet, daß wir uns auf Wellenlängen beschränken, die viel kürzer als Tsind. Unter dieser 
Voraussetzung bleibt schließlich zur Erfüllung von (1) eine Funktion von & allein, die zum 
Maximum zu machen ist: 


To)= 7 ([10 ara) +(Frosimara)‘. BRRIR A U. (A). 


Diese Funktion ist das von Schuster 1905 definierte Periodogramm und löst offenbar die 
eingangs gestellte Interpolationsaufgabe. Ihre einzelnen Maxima geben diejenigen Schwingungen 
an, die (jede für sich) f(t) optimal annähern. Man bestimmt weiterhin die Maxima durch Auf- 
suchen der Nullstellen der Ableitung mittels geeigneter Näherungsmethoden für jedes spezielle 
Problem. Die Rechnungen sind sehr umfangreich. 

Verschiedene Mängel der Methode treten bei unserer Herleitung sofort in Erscheinung: 

1. Voraussetzung T>1. 

Läßt man diese fallen, so findet man nach einiger Rechnung das sog. „vollständige“ 
Periodogramm: 


2. Gegenseitige Störung der Maxima. 


ra en (w) bzw. h(w) zeigen u. U. eine ganze Anzahl m von Perioden an. Jede einzelne Schwin- 
e gung befriedigt (1), doch ist dies für ihre Summe nicht mehr der Fall. Im allein durchgerechneten 


Periodogramm bleiben zu können: 


Differenzschwingung von der Größenordnung T ist. 


Cr — Sg)? +(8p — Cq) N ' 
h(o)— , >)- DELL IR ey 
GR Ppr—P) | 
wobei C und 8 die schon oben gebrauchten verallgemeinerten Fourierkoeffizienten 
T 7 | 
C(o) — [ fit) cos wi dt; S(w) — [ f(t) Bin bat nu RE (6) 
{ f) 0 
sind, und 7, 9, r die folgenden Korrekturgrößen bedeuten: 
sin2wT T sin? 7 T sin2wT 
alt PYAYı Tor); Sa DER alien; m. 


ballm=2 findet man, daß man nunmehr folgende Forderungen stellen muß, um beim Schuster- 


TF>1;: aT>1; 4, lo To ea Da | 
Letzteres besagt, daß die Perioden nicht zu benachbart sein N da sonst die Länge der 
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3. Periodogramm prinzipiell nicht invariant gegen Niveauverschiebung von f(?t). 


Ändert man alle Ordinaten der Ausgangsfunktion f(t) um einen gewissen konstanten 
Betrag c, so ändert das nichts an den diskreten Fourier Amplituden des Spektrums wegen der 
Orthogonalitätsrelationen der trigonometrischen Funktionen in jedem Intervall, das ihre Periode 
ganzzahlig enthält (natürlich mit Ausnahme des uninteressanten Absolutgliedes). Alle nicht- 


i ; } 2 
harmonischen Zwischenordinaten (A = nn 1, 2,....) dagegen erfahren Veränderungen, 
A . 


wie man aus der Gegenüberstellung 


S(w) = [ro sin wi di | 
% ö N ' (9) 
S(w) = J [f(t) + e] sin ot dt = S(w) -| „A — cosoT) | 


ersieht, und entsprechend für C(w) und C(w). Der Zusatzterm kann durch Wahl von c beliebig 


groß gemacht werden. Aber selbst bei der naheliegenden Wahl e = — f(t) (Reduktion auf Mittel- 
niveau 0) kann er für kleine & durchaus die Größenordnung von S(w) erreichen. 


Die dargelegten Verhältnisse wurden am Beispiel der Reihe der Sonnenfleckenrelativzahlen 
illustriert. Das Ausgangsmaterial sind die 2400 (ungeglätteten) Monatsmittel von 1749 bis 1948 
einschließlich. Der gesamte Verlauf der Periodogramme zwischen den Periodenlängen A=5 Jahre 
und A=100 Jahre (,‚Durchmusterung‘“) kann hier nicht wiedergegeben werden, so daß auch die 
Unterschiede zwischen Schusterschem, ‚vollständigem‘ und ‚‚auf Nullniveau reduziertem‘ 
Periodogramm nicht zum Ausdruck kommen. Die Abweichungen sind nur für Wellen länger 
als etwa 30 Jahre beträchtlich. Stattdessen gebe ich hier die aus 7(w) mittels Näherungsmethoden 
bestimmten Maxima, d.h. alle so erkannten Sonnenfleckenperioden an, die, von mathematischen 
Einwänden abgesehen, die z. Z. genauesten sein dürften. Unter « (Jahr) ist hier ein Zeitraum 
verstanden, der entsteht, wenn alle Monatsmittel ohne Rücksicht auf die Monatslänge mit 
gleichem Gewicht in die Rechnung eingehen (,,‚mittleres Jahr“). 


Die Periodizitäten der Reihe der Sonnenfleckenrelativzahlen 


: Ampli- ee: Zeitpunkt des : Ampli- 2 Zeitpunkt des 
Perioden ade Phasen I Mshnrnd Perioden a Phasen ee Maximums 
P7 ERTE tg p vor 30. 12. 1948 [7 TTı948 tg p vor 30. 12. 1948 
« 5,498 13 1,012 93.10.1947 13,204. 9,0 1,044 19. 1. 1936 
5,773 6,0 0,863 27. «1. 1946 14,347 8m: 5,728 7. 8.1948 
6,997 5,4 2,544 3. 10..1947 19,565 6,1 0,846 15. 4. 1947 
7,663 7,2 5,579 20. 5.1947 26,643 9,3 2327 2922. 5.1945 
8,413 12,5 3,349 24. 12. 1946 31,670 5,3 0,690 Ar 47: 1943 
8,847 12,6 2,181 17. 9:1946 37,747 8,7 0,980 16. 8.1943 
9,311 16,8 1,694 16. 1.1947 45,414 12,1 3,411 l. 11. 1941 
9,375 21,8 0,851 |. 4. 11. 1947 82,527 27,6 1,764 29. 3.1937 
11.211 27,6 5,920 19. 12. 1938 


Referate der Vorträge: 
H. Bilharz, Freiburg i. Br.: Über die Gaußsche Methode zur angenäherten Be- 
rechnung bestimmter Integrale, 
H. Eltermann, Braunschweig: Verbesserung der Fehlerabschätzung bei Dilferen- 
tialgleichungssystemen, 


S. T. A. Ödman, Djursholm: Ein Extremalprinzip und dessen Anwendung als 
Entwicklungstheorem für die Lösung linearer Differentialgleichungen, 


sind nicht eingegangen. 
Ein ausführlicher Bericht über den Vortrag W. Richter, Frankenthal: Graphische 


Lösung von gewöhnlichen Differentialgleichungen mit nomographischen Hilfsmitteln, wird 
in einem späteren Heft der Zeitschrift erscheinen. 
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Geometrisches über Potenzbetragflächen 
Von Egon Ullrich in Gießen 


Die allgemeinen Potenzen 2*+%P erwiesen sich, innerhalb einer sehr weiten Klasse analy- 
tischer Funktionen w(z), als ausgezeichnet durch eine Krümmungseigenschaft ihrer Betrag- 
flächen im Großen: Sie sind nämlich die einzigen Funktionen der Klasse, deren Betragfläche im 
Gesamtverlauf einsinnige Gau ß krümmung K zeigt,d.h. daß überall im Existenzgebiet K>0 
oder überall K<O gelten soll; das tritt hier ein, je nachdem «21 gilt, während X=0 überhaupt 
nur für die Potenzen 2!%iP und die Exponentialfunktion vorkommen kann. Eine ausführliche 
Darstellung erscheint im 54. Bande der Math. Zeitschrift; dort finden sich auch Figuren über 
die Betragflächen zu 2*+iß. Jede solche Potenzbetragfläche (kurz Pbfl. &, ß) ist in Zylinder- 
koordinaten r, @, h dargestellt durch» 

= fe Pr, at I EI Der Eee Br 


Sie schneidet die Koordinatenflächen A = const in logarıthmischen Spiralen, = const in 
Parabeln h = cr* (Parabel bezieht sich hier ausnahmsweise auf alle festen a0), r = const 
endlich längs Raumkurven, die entstehen, wenn man eine Exponentialkurve auf diesen Zylinder 
aufrollt ;wir haben sie Exponentialwendeln genannt 
(kurz: Wendel). Nur wenn’ ß verschwindet, sind 
die Pbfl. schlichte Rotationsflächen, die Wendeln 
arten in Kreise aus; für jedes $=+0 aber haben 
wir unendlich vielblättrige Flächen, die man er- 
zeugen Kann, indem man die Parabel cr*in jedem 
Fächerschnitt mit e = e”?? so überhöht, daß sie 
für laufendes längs der Wendel über r=1, 
h= ef? geführt wird. Für «—=1 entsteht der 
Kegel dieser Wendel, mit der Spitze in z=0, der 
aber auch als Kegel einer log. Spirale verstanden 
werden kann; für &=0 ergibt sich das gerade 
Konoid der Wendel h= ef. 


Ich habe die Pbfl. geometrisch näher unter- 
sucht. Als Ausschnitt aus dieser Arbeit berichte 
ich hier über eine Klasse von Raumkurven, die 
mit diesen Flächen eng verbunden ist, und bisher 
nicht bemerkt worden zu sein scheint. Man 
stößt darauf schon, wenn man nach den Kurven 
festen Arguments fragt; es ist argz’tP—y 
längs der Durchdringung der Pbfl. mit dem 
Zylinder Z der log. Spirale 


& y a 
r=exp|— — %_ann2). 
get ganz 

Auf denselben Kurventypus kommt man 

DR RR RER CRS N. aber auch, wenn man zwei beliebige Pbfl. (a,, ß,) 

abirelia ud Erponen äh und (%, ß,) zum Schnitt bringt. Die Durchdrin- 

gung geschieht im allgemeinen Falle längs einer 

Folge von Kurven S,n—=0, +1, +2, ...), die aus einer, etwa dem Hauptschnitt S, durch 

(r, 9, h) = (1, 0, 1)in einfacher Art zu gewinnen sind: Man gehe etwa von den Durchstoßpunkten 

des Fächerschnittes 90 an der Pbfl. x, ß, mit allen Blättern g=2nr von &,, ß, aus. Jede 

dieser Kurven kann auch erzeugt werden, wenn man den Zylinder L einer log. Spirale mit dem 
Konoid M einer Exponentialwendel schneidet. Für S, gibt das: 


L.r= er) PR P—Pı B2%ı = Pı 
M:.h= en] mit Im Te Be ae 2 

Daraus ergeben sich alle Schnittkurven $, der beiden Pbfl. in der Form r„—=r?'e*?, h,—rtıren® 

mit 1, = exp {2 nßz: (& — a,)}. Zwei Pbfl. mit %,—=ß, schneiden sich nur in Wendeln <A,u> 


nn — 
9 . 


Da unsere Kurven mit den log. Spiralen als deren Überlagerungen so eng verbunden sind, 
sei der Name Spirella angeregt. 
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Durch jede Spirella (2) mit den Hochzahlkoeffizienten <A, u> geht eine einparametrige Schar 
von Pbfl. Wir wählen den Realteil des Potenzexponenten, x, als Parameter und finden dazu 
die Schar der Potenzen z*+?(4x=@, bzw. ihre Betragflächen. Darunter gibt es für A+0 genau 
eine schlichte Fläche für «= u:4. Jede Spirella kann also auch erzeugt werden, indem man 
den Spiralzylinder Z..r = e*? mit der Rotations,,paraboloid‘fläche P..h = r# schneidet. 
Diese Erzeugung bietet den Vorteil guter Übersicht; nur die Exponentialwendeln <0, u) ent- 
ziehen sich natürlich der Einlagerung in eine schlichte Fläche über der z-Ebene. 

Für u=4 wird Pin einen Kreiskegel ausarten; die Spirella (A, A) ist dessen Schnittkurve 
mit dem Zylinder Z der log. Spirale, also die aus der elementaren Differentialgeometrie durch 
ihre natürlichen Gleichungen bekannte Konische Schraubenlinie: Krümmung und Windung sind 
dann der Bogenlänge proportional. Für #=+-/ aber haben wir in den Spirellen echte Verall- 
gemeinerungen der konischen Schraubenlinien vor uns. Sie liegen nicht auf einem Kreiskegel. 
Wohl aber kann jede von ihnen auf-dem Kegel einer Exponentialwendel untergebracht werden; 
wir müssen dann aus der Schar der Pbfl. durch die Spirella «—=1 herausgreifen, also P=A—u 
nehmen. Jede Spirella gehört so drei Regelflächen an: Außer dem in (2) genannten Spiral- 
zylinder Z und dem Wendelkonoid M noch dem Wendelkegel K..h= rem Das bietet 
Vorteile im Rahmen der Darstellenden Geometrie. 

Unser Bild zeigt die Durchdringung von K und M und deutet im Grundriß (zur Über- 
sichtlichkeit tiefer gezeichnet) den Kreis und die logarithmische Spirale an, welche aus dem 
Konoid die Exponentialwendel und die X, L, M gemeinsame Spirella herausheben. Der Flächen- 
streifen auf dem Konoid, zwischen Spirella und Wendel ist hervorgehoben. Um mehrere Umläufe 
in der Zeichnung unterzubringen wurden die Koeffizienten A, a passend, klein gewählt. 

Dem Charakter dieses Vortragsauszuges entsprechend sind entbehrliche Konstanten bei- 
seite gelassen, Ausartungen nur angedeutet worden. Eine ausführlichere Darstellung ist für 
die Monatshefte der Mathematik vorgesehen. 


Der Zusammenhang allgemeiner Eigenwertaufgaben bei 
Differentialgleichungen mit der Integralgleichungstheorie 


Von N. Joachim Lehmann in Dresden 


Viele technische Eigenwertaufgaben führen zu Differentialgleichungen 
m Nn\ı< m) 
a3 (IP (py=i 23 (-1rHy®)” kurz My)=AN(y) ....() 
v=0 v=0 


mit in y und A linearen Randbedingungen: U,(y)=A8,(y), u=1, ..., 2m. 

Spezielle Aufgaben hiervon — z.B. wenn N (y) eingliedrig ist — lassen sich mit einer 
Greenschen Funktion der Theorie symmetrischer Integralgleichungen einordnen. Im allge- 
meinen Fall schien aber dieser einfache Weg versperrt. Es soll gezeigt werden, wie hier die 
Einführung einer geeigneten Maßform weitergeführt; explizit werden dabei die Integral- durch 
Integrodifferentialgleichungen ersetzt. Das Vorgehen sei unter den Kamkeschen Selbstadjun- 
giertheitsforderungen zu (1) skizziert. Dabei ist im Grundgebiet a...b $S,=0, f„ +0. Zum 
Operator M soll die (symmetrische) Greensche Funktion @(&, s) existieren. (1) ist damit äqui- 
valent zu 


y(a) =A[G@ (x, s) N(y(s))ds. (Integrationen immer über @...b verstanden.) : . . (2). 
Für Lösungen gewinnt man durch partielle Integration mit Died Rd 
[ ng» nl 2 ik or r | 
va—a8 55ER IT I A we 
vo 08? v,0=0 i,k-1 os i 


Dabei konnten nach Kamke mittels der Randbedingungen in den Randausdrücken alle Ablei- 
tungen n-ter und höherer Ordnung hinausgeworfen werden. 

Gl. (2)’ dient als Ausgangspunkt für eine Verallgemeinerung der üblichen Integralglei- 
chungen. Diese werden durch Integralformen Jw(®) v(2)dz bzw. [u?(z)da beherrscht, die 
jetzt durch allgemeinere Bilinearformen 


n f y 
(w)= 3 I [umvWdon. +3 u) 2222.) 
v‚e=0 ü,k J 
“ersetzt werden. Dabei seien o,, gegebene Belastungsfunktionen beschränkter Schwankung, 
a’ sind feste Konstanten, &;, x; endlich viele Punkte aus dem Grundgebiet. Die Form (u, v) 


4) Es dürfen auch Glieder z. B. der Form ff Fye(z,5)u®)(x)v(e(s) dx ds dazukommen. 


ne > = Er eo 
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kann z.B. für alle in @...b stückweise stetigen Funktionen, die ebensolche Ableitungen bis 
zur n-ten Ordnung besitzen, gebildet werden — diese Gesamtheit soll als Grundbereich (oder 
Raum) N bezeichnet werden. Mit jedem Kern K(z,s) €”. können nun ‚„‚Integralgleichungen“ 
hingeschrieben werden, die (2)’ umfassen: 


y(o)=AlKl,s),yle))+fle) » >> N (4). 
EN ist eine Störfunktion, (K(z, s), y(s)) wird bzgl. der Variablen s gebildet. 


Hier interessiert besonders der Fall, daß K(x,s) symmetrisch ist und die Ableitungen 
or+e 
0x’ 05° 

Man erhält das wichtige Ergebnis: 

Gehört solch ein Kern für jede Variable zu einem Funktionsbereich ® (d.h. einer linearen 
Funktionengesamtheit) aus N, für welchen die Form (u, v) u, vED symmetrisch und positiv- 
definit ist?), so läßt sich auf (4) sofort die bekannte Integralgleichungstheorie für symmetrische 
Kerne ausdehnen. Wie früher lassen sich Existenz und Abzählbarkeit der Eigenwerte, Ent- 
wicklungssätze für iterierte Kerne und quellenmäßig dargestellte Funktionen beweisen. (Bei 
den Entwicklungssätzen ist die Konvergenz sogar der n-ten Ableitung noch gleichmäßig) ?). 

Gibt es keinen Bereich ®, der K(z,s) enthält, so liegt eine Verallgemeinerung der polaren 
Integralgleichung vor. Auch sie läßt sich bei symmetrischem (u, v) und definitem Kern?) sofort 
der gleichen Beweisführung unterordnen — die störende Diskrepanz bei der Behandlung ge- 
wöhnlicher und polarer Integralgleichung verschwindet. Dazu wird die Ausgangsgleichung 
iteriert. Bei z.B. der homogenen Be erhält man 


y(2)=R(K(a,s), (K(s, 8), yd)))=u<K(,s),y(s)) - SR Sata )® 


K(z, s) für v,0—=(, ..., n im Grundgebiet in beiden Variablen gleichzeitig stetig sind. 


und es kann 
us), vo) = (u), (E68), 0 E)) = (RC, ©), ul), vO)) 


als neue in ganz X symmetrische und definite Bilinearform verwendet werden. In dieser Maßform 
gehört das Problem wieder zu (4) und es it X = D und K(z,s)€®. Die hierfür gültigen Er- 
gebnisse können sofort auf die Ausgangsgleichung zurückübertragen werden. 

Die allgemeine Theorie zur Gl. (4) umfaßt über Gl. (2)’ alle Kamkeschen selbstadjun- 
gierten links (M)- oder rechts (N)-definiten Randwertprobleme der Gl.(1)°). An Hand der 
Extremaldarstellung der Eigenwerte erkennt man sofort, daß sogar alle (M)- oder (N)-definiten 
„natürlichen Eigenwertaufgaben‘“®) erfaßt werden. Von den speziellen Ergebnissen sei nur ein 
Entwicklungssatz erwähnt. Er umfaßt bei absolut gleichmäßiger Konvergenz bis zu den 
(m —1)-ten Ableitungen alle Funktionen, die ebensoviele stetige und eine stückweise stetige 
m-te Ableitung besitzen und die wesentlichen Randbedingungen (nach Kamke) erfüllen. 

Weitere Anwendungsgebiete der Gl. (4) lassen sich nur andeuten. Es werden z.B. die 
Aufgaben ’”) 


M(y) + [F(«, s)y(s)ds =A {IN (y)-+ [E(z, 3) Ya) BEN ae ee (5) 
ebenso eingeordnet wie die Gleichung für die Schubbeulung einer Platte: 
> £ 0°2 ni A ER 
Ad en, mit BF er 0 am Rand (Einspannungsfall) . . (6). 


Bei (6) (einem Beispiel für zweidimensionale Aufgaben) kann wieder die bei vernünftigen Ge- 

bieten existierende positiv-definite Greensche Funktion J/'(z, s; &,n) ausgenutzt werden. Man 

kommt nach Iteration auf 

2 I'(2,y; &n) ou Ov 
PRArE ay yon 


(dr yon) mit ||? 


und. das läßt sich wieder der Theorie zu Gl. (4) einordnen. 


— dadyd£&dn (6) 


2) Bei z.B. (u, v) =/gu"v"dae-+fhuvdz, gh>O0 ist einfach DEN. 

,) Bei unsymmetrischen Kernen bleibt die Fredholmsche Theorie gültig. 

4) Es genügt sogar, wenn erst die Form <u,») nach Gl. (4)’ in einem Bereich ®,, dem wieder K(z, s) 
angehören muß, symmetrisch und positiv-definit ist. 

5) Es genügte allerdings die Behandlung (M)-definiter Aufgaben, (N)-definite Probleme NN zufolge 
fm==0 nur für ein Vorzeichen oo viele Eigenwerte. Eine Spektralverschiebung führt daher immer zu (M)- 
definiten Formen. 

6) E. Stiefel und H. Ziegler: Natürliche Eigenwertprobleme. ZAMPT (1950), S. 111—138. 

?) L. Collatz: Z. angew. Math. Mech. 25—27 (1947), S. 129—130. 


u ET Ve 


NN 1951 6. Mathematische Maschinen und Instrumente 253 


Zusammenfassend läßt sich feststellen, daß die Auffassung der Integralgleichung als 
‚„‚Bilinearform mit Kern“ Erweiterungen ermöglicht, welche die Mehrzahl aller praktischen 
Eigenwertaufgaben umfassen®). Dabei bleiben die eleganten und elementaren Beweisführungen 
der klassischen Integralgleichungstheorie erhalten. (Im Hinblick auf numerische Fragen wurde 
vom Vf. die konstruktive Beweisanordnung von E. Schmidt bevorzugt.) 


Referate der Vorträge: 


H. Friedman, Oxford/Heidelberg: Zur Mathematik der Gestalt, 

H. Berger, Wien: Eine Erweiterung des Abbildungssatzes von H. A. Schwarz auf 
zweifach zusammenhängende Bereiche, 

J. Dörr, Darmstadt: Über zwei Integralgleichungen 1. Art, die sich mit Hilfe 
Mathieuscher Funktionen lösen lassen 


sind nicht eingegangen. 


C. Mathematische Maschinen und Instrumente 


Einige neue mechanische Integrierseräte 
Von Th. Erismann in Schaffhausen 
(A. J. Amsler & Co. in Schaffhausen, Schweiz) 


Bei der Verarbeitung unregelmäßiger Funktionen und insbesondere bei der Auswertung 
registrierter Meßwerte sind die Analogiemethoden den heute sehr beliebten digitalen 
Methoden grundsätzlich überlegen, da die letzteren in diesen Fällen einer bedeutenden Vor- 
bereitungsarbeit bedürfen. Gerade bei derartigen Aufgaben spielt die beschränkte Genauigkeit 
der Analogiemethoden meist keine Rolle. 

Obwohl die Entwicklung der mechanischen Integriergeräte (wohl der wichtigsten 
Analogiegeräte) heute schon über hundert Jahre alt ist, so ist sie dennoch durchaus nicht ab- 
geschlossen. 

So konnte das Planimeter in den letzten Jahren beträchtlich verbessert werden: Das 
übliche Schneckengetriebe mit der Zählscheibe zum Zählen der vollen Rollenumdrehungen wurde 
durch ein Zahnradgetriebe ersetzt, das die 
koaxiale Anordnung zweier Ablesetrommeln 
gestattet (Bild 1). Dieses sog. Integrier- 
zählwerk kann also fast wie ein Zähler ab- 
gelesen werden, was die Ermüdung und die 
Wahrscheinlichkeit von Ablesefehlern stark 
verringert. Der Preis liegt dank rationeller 
Herstellung nicht höher als bei älteren 
Planimetern. 


Naturgemäß wurde das neue Integrier- 
zählwerk für eine Anzahl weiterer Geräte ver- 
wendet, so z.B. für ein Scheibenpolar- 
planimeter, eine Reihe von Potenzinte - 
gratoren und ein Integrimeter. Beidie- 
sem letzteren ist es dank der geschickten 
Anordnung möglich, das Integral fv ds im 
gesamten Bereich —250 mm < y< --250 mm 
laufend zu vermessen. 

Eine andere Reihe von Geräten beruht 
auf dem 1886 durch Hele Shaw erfundenen Kugelintegrator, z.T. auch auf dem ent- 
sprechenden Kugeldifferenziator von J. Amsler. WER 

Das einfachste dieser Geräte ist ein Integraph zur laufenden Registrierung des Integrals 
[ydas (Bild 2). Bei diesem Instrument erhält die Kugel ihren Antrieb nach dx durch unmittel- 
bares Abrollen auf dem Lineal. Ihre Drehachse wird von zwei Steuerrollen bestimmt, die vom 
Fahrstift F über einen Pantographen und eine Tangentenmechanik geschwenkt werden. Das 
Resultat wird von der bei J gelagerten Integrierrolle über eine Reibungskupplung und ein 


Bild1i 


8) Der Bereich der vollstetigen Transformationen wird noch nicht verlassen. 
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Zahnrad auf die Zahnstange Zübertragen, die die Schreibfeder Strägt. Wird der Hub der Schreib- 
feder erschöpft, so kann diese unter Überwindung der erwähnten Kupplung in die entgegengesetzte 
Endlage zurückgeschoben werden. Damit ist eine Auswertung praktisch beliebig langer Dia- 
gramme möglich. Das Erreichen der Endlagen wird durch ein akustisches Signal angezeigt. 

Auf Grund seiner besonderen Eigenschaften (gute“Genauigkeit, Übertragung eines be- 
trächtlichen Drehmomentes ohne Nachlaufsteuerungen) wird das Kugelgetriebe seit langem mit 
Erfolg bei Integrieranlagen verwendet. Die erste davon war der Registriertisch für 
Dynamometerwagen mit vier Kugelgetrieben, der seit 1913 in vielen Exemplaren fast 
unverändert hergestellt wurde. Es handelt sich hierbei um eine ausgesprochene Einzweck- 
Anlage, d.h. um eine Anlage mit einer einzigen unveränderlichen Schaltung. 


Bild 2 Bild 3 


Die guten Erfahrungen, die mit diesen Registrierlischen gesammelt werden konnten, 
wurden in den letzten Jahren an einer Reihe weiterer Einzweck-Anlagen ausgewertet. Als ver- 
kleinerte und modernisierte Ausgabe des Registriertisches entstand ein Prüfaggregat für 
Triebwagen und Traktoren. 

Ferner ein Fahrdiagraph, der nach der bekannten Knorr-Schaltung arbeitet. In Bild 3 
sind bei J, und J, die beiden Integratoren (Beschleunigung-Geschwindigkeit und Geschwindig- 
keit-Zeit) sichtbar, bei © die Zugs-Charakteristik, bei R das Zifferblatt zum Einstellen der 
Rampen, bei Z, und Z, die beiden Zähler für Zeit und Strecke und bei P das Papierband, auf 
dem vor Beginn der Arbeit die Rampen eingetragen wurden und das zur laufenden Registrierung 
der Geschwindigkeit und der Zeit in Funktion der Strecke dient. Eine Reihe von Hilfseinrich- 
tungen ermöglicht ein außerordentlich einfaches Arbeiten, so vor allem eine automatische 
Kupplung K, die dasÜberschreiten einer beliebig einstellbaren Geschwindigkeitsgrenze verhindert. 

Während des zweiten Weltkrieges wurde für die Schweizer Armee.eine Anlage zur Be- 
rechnung von Geschoßflugbahnen gebaut, welche sowohl die Aufzeichnung der Flugbahnen 
als auch die Aufnahme der Schießtabellen erlaubt. Durch besondere Maßnahmen konnte hierbei 
eine hervorragende Gesamtgenauigkeit erreicht werden. 

Die erwähnten Einzweck-Integrieranlagen gestatten bei verhältnismäßig bescheidenem 
Aufwand eine sehr gute Anpassung an die gestellten Aufgaben. Ihre Wirtschaftlichkeit ist da- 
durch bedingt, daß die betreffenden Aufgaben immer wieder neu gelöst werden müssen. Wo dies 
nicht der Fall ist, kann eine universelle Integrieranlage interessant sein. Eine solche steht 
gegenwärtig in Entwicklung, z.T. schon in Erprobung und dürfte noch im Laufe dieses Jahres 
die ersten praktischen Aufgaben bearbeiten. 


Erste Erfahrungen mit dem programmsgesteuerten 
Rechengerät Z4 von K. Zuse 


Von Heinz Rutishauser in Zürich 


Der Vortragende bespricht einige numerische Probleme, welche am Institut für angewandte 
Mathematik der ETH (Zürich) mit Hilfe der Z4 gelöst wurden und hebt einige Fälle hervor, 
wo die Verwendung einer programmgesteuerten Rechenmaschine eine Modifikation der Lösungs- 
methoden erfordert (Inversion von Matrizen und Lösung der Bipoientnigi Ze durch Relaxa- 
tionsrechnung). 

Der Inhalt des Vortrags wird im Rahmen einer Mitteilung des Instituts für angewandte 
Mathematik der ETH ausführlich publiziert werden, 


7. angew. Math. Mech. 2 h 
Bd. 31 Nr. 8/9 Aug./Sept. 1951 D. Stochastik 298 


Über automatische Rechenplanfertisung bei programmsgesteuerten 
Rechenmaschinen 


Von Heinz Rutishauser in Zürich 


Kurzer Bericht des Vortragenden über seine Untersuchungen. Es handelt sich darum, die 
Aufstellung eines Rechenplans für eine komplizierte Formel (oder Formelgruppe), welche oft 
eine erhebliche Arbeit ist und für welche die Konstruktion komplizierter Hilfsgeräte in Erwägung 
gezogen wird, durch die Rechenmaschine selbst ausführen zu lassen. Die beschriebenen Methoden 
ermöglichen auch das Rechnen im Körper der rationalen Funktionen mit einer programm- 
gesteuerten Rechenmaschine. 

Ein ausführlicher Bericht über diese Untersuchungen wird als Mitteilung des Instituts 
für angewandte Mathematik der ETH erscheinen. 


Ein Referat des Vortrages: 


A. Kuhlenkamp, Braunschweig: Neue Wege der maschinellen Lösung elemen- 
tarer Rechenaufgaben bei hoher Stellenzahl 


ist nicht eingegangen. 


D. Stochastik 


Die Zuverlässigkeit üblicher Vaterschaftsteste 
Von W. Ludwig in Heidelberg 


Der Vaterschaftstest zielt auf den Entscheid der Fragen ab, ob ein vorgegebener Mann (der 
„mögliche Vater“ V) als ‚‚wahrer Vater“ (V*) eines vorgegebenen Kindes (K), bzw. welcher von 
mehreren V als plausibelster V* erklärt werden kann. Hierzu gibt es zunächst denegativen 
(Ausschluß-)\ Verfahren, z.B. daß, wenn K und Mutter die Blutgruppe 0 aufweisen, der V+ 
nicht AB besitzen kann. Indes würden diese Methoden, auch wenn man alle bekannten streng 
erblichen Merkmale einbezöge, nur in einem praktisch noch ganz unzureichenden Ausmaße ver- 
antwortbare gerichtliche Entscheide gestatten. Ihnen gegenüber beruhen die positiven 
Verfahren auf überdurchschnittlicher Übereinstimmung zwischen K und V*+ gegenüber K 
und beliebigem V hinsichtlich irgendwelcher Merkmale, — geschaffen von Essen-Möller 
(1937f.). Möge ein K mit einem beliebigen Mann der Bevölkerung durchschnittlich im Bruch- 
teil 4, der Fälle bez. Merkmal M; übereinstimmen, hingegen K mit seinem V* im Bruchteil «,, 
wobei meist 2; > %; sein wird, so läßt sich bei Betrachten von » Merkmalen und für den Fall nur 
eines nominierten V die Formel 


w=1:[14+0-2(%) NE N 


aufstellen, in der bedeuten: W die (der Wahrscheinlichkeit proportionale) Plausibilität (likeli- 


hood), daß der nominierte V der V+ sei; = en), wobei w = Häufigkeit, in der eine Mutter 


als V den V+* nominiert; und J/ das Produkt der » Übereinstimmungsquotienten der als un korre- 
liert vorausgesetzten Merkmale. (Bei Nichtübereinstimmung bez. M, ist x; durch 1 — x; und 
y, durch 1 — y; zu ersetzen.) Gegen diese Formel sind Bedenken erhoben worden, jedoch zu 
einem Großteile unerheblicher, unrichtiger oder unscharf begründeter Art, weitere noch nie. Zu- 
nächst fehlt inEssen-Möllers Originalformel unser @ — d.h. es ist = 1 gesetzt —, weil 
für einen „‚unbefangenen Gutachter‘ ein nominierter V ebensogut der V* sein könne oder nicht. 
Indes ist, mindestens in städtischen Bezirken, »>> 0,5, also Q<1. Betrachten wir nur I nicht- 
erbliches Merkmal (also & = %) und setzen w = 75%, so folgte nach Essen-Möller für den 
Beklagten W = 0,5, während die richtige Formel (1) notwendig W = w = 0,75 liefern muß. 
_ Über das Ausmaß der durch Ignorierung irgendwelcher Korrelationen der M, untereinander 
hatte sich bisher noch niemand geäußert. Berechnungen!) ergeben: Stimmen z.B. V und Kin 
2 Merkmalen mit je & = 17%, y = 3% überein, so variiert das zugehörige W, sinkend von absolut 
positiver über fehlende zu absolut negativer Korrelation, zwischen 85% — 97% — 100%. Schließ- 


2) Vgl. W.Ludwig, Homo (Mainz), i. Dr.; R.Wartmann,i.Dr.; W. Ludwig u.R. Wart- 
mann,i.Dr. 
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lich ordnete man bisher diese Plausibilität in die Intervalle 1 — 0,9973 — 0,955 — 0,682 
— 0,318 — 0,045 — 0,0027 — 0 ein und schrieb ihnen die Bezeichnungen,, Vaterschaft praktisch 
erwiesen“, „sehr wahrscheinlich‘, ‚‚wahrscheinlich‘“ usw. bis ‚„‚praktisch ausgeschlossen‘ zu, die 
dem Richter als Unterlage dienten, — ohne zu bedenken, daß diese für T o t al wahrscheinlich- 
keiten geschaffenen Spannen nur sehr bedingt auf Einzelwahrscheinlichkeiten übertragen werden 
dürfen. Schließlich wurden verschiedentlich, wenn sich zunächst ein niedriges W ergab, solange 
weitere Merkmale beizuziehen versucht, bis W eine entsprechende Höhe erreicht hatte, nicht 
aber, wenn sich schon anfänglich ein hohes W errechnete: ein weder statistisch noch juristisch 
einwandfreies Vorgehen. Andere Bedenken seien hier unterdrückt. Da sich die negativen Ver- 
fahren vorerst zu wenig effizient erweisen, verbleibt wohl nur die Möglichkeit sorgfältiger Einzel- 
gutachten, ohne einen falsche Exaktheit vortäuschenden scharfen Zahlwert, oder die 1936 von 
R. A. Fisher inaugurierte Diskriminanzanalyse?), die unserem Problem voll ge- 
recht würde und je Distrikt nur eine einmalige und nicht sehr zeitraubende Vorarbeit erforderte, 
überdies den Test sehr verfeinern könnte. Daß sie für den Vaterschaftsnachweis bisher noch nicht 
angewandt wurde, beruht darauf, daß es diesen in jenen Ländern (England, Indien, USA), die 
dieses Verfahren kennen, im europäischen Sinne nicht gibt. 


© 


Die statistische Trennung sich in mehreren Merkmalen 
überlappender Individuengsruppen (Diskriminanzanalyse) 


Von R. Wartmann in Darmstadt 


Häufig ergibt sich die Aufgabe, einzelne Elemente einzuordnen in eine von zwei oder 
mehreren Individuengruppen. Da diese Klassifikation oft mühsam oder mit der Zerstörung des 
betreffenden Elementes verbunden ist, ist man bestrebt, diese Einordnung zu ermöglichen mittels 
einfach zu erhaltender Merkmale. 

Von einer Anzahl Elemente mit bekannter Gruppenzugehörigkeit (etwa durch Anwendung 
des mühsamen Verfahrens) werden mehrere Maße abgenommen und deren Verteilung in den 
einzelnen Gruppen (genauer: Stichproben) aufgestellt. Ein Einzelmaß gibt ein Kriterium für die 
Gruppenzugehörigkeit nur dann, wenn die Verteilungen dieses Maßes in den einzelnen Gruppen 
sich sehr wenig überlappen. Man bildet daher aus all den Einzelmaßen ein neues ‚„‚Rechenmaß‘“, 
für das dies möglichst zutrifft. Letzteres entsteht durch Addition der mit Faktoren multipli- 
zierten Einzelmaße, wobei die Faktoren derart gewählt sind, daß die Überlappungen der Rechen- 
maß-Verteilungen der Gruppen ein Minimum werden. Für ein zur Entscheidung stehendes 
Element braucht man dann lediglich auf Grund seiner Einzelmaße seinen Rechenmaßwert zu 
. ermitteln. — Zur Berechnung der Faktoren ist erforderlich die Kenntnis der Streuungen und 
Korrelationen der (Stichproben-)Einzelmaße, die Bestimmung der größten Wurzel eines in 
Determinantenform gegebenen Polynoms und die Auflösung eines linearen Gleichungssystems 
(Grad des Polynoms und Zahl der Unbekannten gleich Zahl der verwendeten Einzelmaße). Bei 
nur zwei Gruppen entfällt das Polynom. — Es ist möglich, die Wahrscheinlichkeit einer Miß- 
klassifikation anzugeben und ferner die Signifikanz der Stichprobenunterschiede (bzw. Gruppen- 
unterschiede) zu prüfen, falls diese nicht a priori festliegen. Schließlich kann man aus der Lage 
der Rechenmaß-Verteilungskurven der einzelnen Gruppen zueinander auf deren gegenseitige 
Verwandtschaft usw. schließen. Dieses ist meist auch beim Vorliegen qualitativer Merkmale 
möglich. 

Neben dem linearen Rechenmaß ist (bei etwa doppeltem Arbeitsaufwand) eine quadratische 
Diskriminanzfunktion entwickelt. Letztere fällt auf die lineare zurück bei gleichen Streuungen 
in den Gruppen. Diese quadratische Diskriminanzfunktion geht von der Annahme aus, daß die 
in den Einzelmaßen mehrdimensionalen Simultanverteilungen der Gruppen Normalverteilungen 
sind und ordnet ein zur Entscheidung stehendes Element derjenigen Gruppe zu, für die die 
betreffende Merkmalskombination die größte Wahrscheinlichkeit hat. 

Da die Güte der Klassifikation, aber auch der Arbeitsaufwand mit der Zahl der verwendeten 
Einzelmaße wächst, wurde eine Approximation aufgestellt, bei der der Fehler, den man durch 
ihre Anwendung begeht, gering ist gegenüber dem Gewinn, der erzielt wird durch die nun mögliche 
größere Zahl der Einzelmaße. Sind alle Korrelationskoeffizienten gleich, ist die Approximation 
mit der exakten Methode gleichwertig, andernfalls verwendet sie (automatisch) einen optimalen 
Durchschnittswert. Der vereinfachte Weg ist der folgende : Normierung der Einzelmaße, Addition 
zum Zwischenmaß ‚‚size‘ und Addition nach Multiplikation mit Gewichten zum Zwischenmaß 
„shape“, Aufstellung der linearen oder quadratischen Diskriminanzfunktion. Die Gewichte 
ergeben sich auf einfache Weise aus den Mittelwert-Differenzen der Einzelmaße. Für ein frag- 


?) Hierzu vgl. den folgenden VortragR.Wartmanns, 
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liches Element erhält man z. B. somit als lineares Rechenmaß: y=b-a--p (a: size, p: shape, 
b: aus der Extremalbedingung errechneter Faktor — der zweite Faktor ist gleich eins gesetzt). 
Die Klassifikation wird getroffen, je nachdem y größer oder kleiner als ein y, ist (bei nur zwei 
Gruppen). In einem a/p-Koordinatensystem stellt aber p= — b-a-+- y, eine Gerade dar. Liegt 
diese einmal gezeichnet vor, so klassifiziert man einfach dadurch, daß man feststellt, ob der 
durch den a- und p-Wert des Elementes gelieferte Bildpunkt auf dieser oder jener Seite der 
Geraden liegt. Die quadratische Diskriminanzfunktion liefert entsprechend als Grenzkurve eine 
solche zweiten Grades. 


Poissonsche Verteilungskurven bei Versuchen mit Drillmaschinen 
Von Hermann Blenk in Braunschweig 
(Aus dem Inst. f. Landtechn. Grundlagenforschung in der Forschungsanstalt für Landwirtschaft) 


Drillmaschinen sollen das Saatgut in gleichmäßiger Tiefe und in regelmäßiger Folge im 
Ackerboden ablegen. Um die Regelmäßigkeit der Körnerablage zu prüfen, werden mit den Drill- 
maschinen im Laboratorium sog. Leimstreifenversuche angestellt. Die Körner fallen dabei auf 
einen Streifen Papier, der mit Leim bestrichen ist und die Körner an der Stelle festhält, wo sie 
auf das Papier auftreffen. 

Die Gleichmäßigkeit der Tiefenablage kann man mit dem Leimstreifenversuch nicht prüfen. 

) Man verzichtet meist auch auf eine Prüfung der horizontalen Komponente quer zur Bewegungs- 
richtung; d.h. man projiziert alle Körner auf eine Gerade in Richtung der Furche und prüft 
nur die Regelmäßigkeit der Körner in dieser Reihe. 

Für die Beurteilung der Regelmäßigkeit einer Körnerreihe sind wiederholt Vorschläge 
entwickelt worden, so von R. Hege [1], A. König [2] und H. Zödler [3]. Man untersucht die 
Verteilung der Körnerabstände statistisch und leitet daraus eine mehr oder weniger empirische 
Formel für die Regelmäßigkeit ab. Statt die Körnerabstände zu bestimmen, kann man auch den 
ganzen Streifen durch Felder gleicher Breite b überdecken und nun auszählen, wie viele Felder 
gerade n Körner enthalten. Bezeichnet man mit h die relative Häufigkeit der Felder mit n 
Körnern, so müssen die Kurven Äk(n) sich mit den bekannten Verteilungskurven von Poisson 
(Wahrscheinlichkeit seltener Ereignisse) 


| __ Zahl der Körner 
“Zahl der Felder 


vergleichen lassen. Dieser — m. W. erstmals angestellte — Vergleich brachte nun bei der Aus- 
wertung von Leimstreifenversuchen mit einer normalen Drillmaschine ein verblüffendes Ergebnis, 
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z Normale Drillmaschine , Saatgut: Zuckerrübe 2 Gleichstandsdhillmaschine, Saatgut : Weizen 
i ‚Rh = relative Häufigkeit der | Gesamtzahl der Körner «281 * h=relative Häufigkeit der | Gesamtzahl der Körner:652 
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Bild 1 Bild 2 


nämlich eine auffallend gute Übereinstimmung der experimentellen mit den Poissonschen 

Kurven, wie Bild 1 als Beispiel zeigt. Da die Wahl der Breite 5 eine Willkür bedeutet, ist die 

Auswertung für drei verschiedene Werte von b (2, 4 und 6cm) durchgeführt worden. In allen 

drei Fällen stimmen die experimentelle und die Poissonsche Kurve gut überein. Daß sich 

_ für b=6 cm etwas größere Unterschiede zeigen als für 2 und 4 cm, liegt sicherlich daran, daß 
mit wachsendem b die Zahl der Felder je Streifen geringer wird. b ® 

Bild 2 zeigt das Ergebnis eines Versuchs mit einer sog. Gleichstandsdrillmaschine. Durch 

eine besondere Vorrichtung ist für eine Bade Regelmäßigkeit beim Abgang der Körner gesorgt. 


* 
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Der Vergleich mit den Poissonschen Kürven bestätigt den Erfolg der Konstruktion: die Werte n 
in der Nähe des Durchschnitts » sind häufiger geworden, während die größeren Abweichungen 
seltener auftreten. 

In Bild 3 sind eine Reihe von Versuchen dargestellt, bei denen nebeneinander eine normale 
und eine Gleichstandsdrillschar benutzt wurden. Über der Felderbreite b ist das Verhältnis 
Streuungsquadrat 0? 

Durchschnitt 2 
Für die normale Drillschar ergeben sich Werte bei 1 und sogar merklich darüber, während die 
Werte für die Gleichstandsdrillschar um 0,6 herum streuen. An der Überlegenheit der Gleich- 
standsdrillschar kann also kein Zweifel bestehen. 

Auf Grund dieser Untersuchungen könnte man zur Beurteilung der Güte einer Drill- 

13 \ In 
maschine die Größe Z=1—a-% verwenden, wobei & = BrunBeverlat u 
n Ertrag bei 0 —=0 


und durch biologische Versuche zu er- 


aufgetragen. (Für die Poissonsche Verteilung sind o? und % beide —=4.) 


zu setzen 


Vergleich der Drillschare mitteln wäre. 

normale Schar j,= mittlere Körnerzahl im Feld Die Streuung 0° wird, wie Bild 3 

———— bleichstandsschar 02= Streuungsquadrat se R 
3 3 5 Rüben zeigt, von der Wahl der Felderbreite b 

Weizen Weizen Weizen Gerste Rüben Rüben Erbsen A Ben E ne £ : 
15 - 15 nicht unwesentlich beeinflußt. Um diese 
Gas IT ee 8 
o2l ||| Willkür!) auszuschalten, wäre es er- 
HT Im! wünscht, die Verteilung der Körnerab- 
„1 stände mit einer entsprechenden theore- 
Lischen Verteilung vergleichen zu können. 
\ “tt Diese Vergleichsverteilung könnte man 
05 IR sich auf folgende Weise verschaffen: Man 
ordnet die Körner zunächst in gleichen 
L Abständen « voneinander an und ver- 


schiebt dann jedes Korn nach vorn oder 

hinten um einen Betrag «. Für dieses & 

wird eine Gaußsche Verteilung mit be- 

Bu stimmter Streuung s? (bezogen auf a) 

angenommen. Bei dieser Verschiebung 

nehmen die Körner eine neue Reihenfolge an. Es käme nun darauf an, dieser Vorschrift 

entsprechend die Verteilungsfunktion der Abstände aufeinanderfolgender Körner, abhängig 
von s/a, zu berechnen. 


RENTE, ZIHD II 02 
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E. Mechanik 


Die Trägheitspolkurve, ein Beitrag zur Dynamik des ebenen Systems 
Von Karl Federhofer in Graz 


Die Trägheitspolkurve entsteht als Ort der Trägheitspole eines kontinuierlich bewegten 
ebenen Systems. Mit der Kenntnis der Lage des Drehpoles P, des Wendepoles J und mit den 
bekannten dynamischen Eigenschaften der der Scheibe (Schwerpunkt $S und Trägheitshalb- 
messer %s) ist die Lage des Trägheitspoles 7 bestimmt, und zwar gilt hierfür die einfache 
Vektorgleichung: 

5 
1p*15 


IT=—Yj 


!) Die weitere Willkür, die in der Wahl des Anfangspunktes der Felder über den Leimstreifen steckt, 
dürfte kaum einen nennenswerten Einfluß haben können. 


» 
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wo die den Punkten T', P, J entsprechenden Ortsvektoren t auf den Schwerpunkt S als Aufpunkt 
bezogen sind. 

Man kann die Trägheitspolkurve auf ein festes oder auf das bewegte System beziehen, 
analog wie die Rast- und Gangpolbahn des Drehpoles, bzw. wie die Kurven für den Beschleu- 
nigungspol. 

Bei Anwendung auf das zentrische Schubkurbelgetriebe zeigt sich, daß die Trägheitspol- 
kurve der Koppel, bezogen auf die bewegte Koppel, mit großer "Näherung durch eine Gerade 
senkrecht zur Koppelachse ersetzt werden kann, auf der innerhalb eines ganz bestimmten Be- 
reiches die Trägheitspole für den vollen Umlauf der Kurbel liegen. 

Ausführliche Darstellung erscheint im Bd. V (1951) des Österr. Ing.-Archivs. 


Die stabilen Lagen eines n-fachen ebenen Pendels mit vertikal 
periodisch erschüttertem Aufhängepunkt!) 
Von W. Haacke in Braunschweig 


Stabilisierungsphänomene eines erschütterten einfachen Pendels und seines elektromagne- 
tischen Analogons sind in den letzten Jahrzehnten mehrfach behandelt worden. Gerade in 
Hinblick auf die elektrotechnischen Abbildungen liegt es nahe, ein mehrfaches ebenes Pendel mit 


erschüttertem Aufhängepunkt zu betrachten. Als erster Schritt auf diesem Wege soll hier der 


Fall der vertikalen Erschütterung eines n-fachen ebenen physikalischen Pendels untersucht 
werden. Wir erhalten ein Differentialgleichungssystem, das sich durch eine lineare Transfor- 
mation entkoppeln läßt. Auf die einzelnen entkoppelten Differentialgleichungen läßt sich dann 
die bekannte Theorie der Mathieuschen Differentialgleichung, besonders die der Struttschen 
Karte anwenden. 

Das ebene Pendel schwinge in einer &,y-Ebene. ®, sei der Winkel, den der o-te Pendelarm 
mit der positiven y-Achse bildet, die Erdbeschleunigung g wirke in negativer y-Richtung. Es 
ergibt sich ein System von n gekoppelten, nichtlinearen a EnE n 2. Ordnung. Wir 


wollen nach der Existenz von Mittellagen &, = lim * 2 D,(r)dr fragen, um die die Pendelarme 
tab. 


Schwingungen ausführen sollen. Setzen wir 9, = D,—o, und beschränken uns auf die Betrach- 


tung solcher kleinen Schwingungen, so läßt sich das Differentialgleichungssystem für die 9° 


linearisieren und erhält folgende Gestalt: 


HuoR ti B, SE 008 (0. — 0x) Pate 2, B, cos (&, — 0%) 95 — B, (25 —p eos u) cos a, Zn 


em, (Gele n). 


Darin sind die A, und B, positive Konstante, die sich aus den Massen, Längen und Trägheits- 
momenten der Pendelarme zusammensetzen, l; der Abstand der Drehpunkte des s-ten Pendel- 
arms und y=p cos wi die Erschütterung des Aufhängepunktes. Die unabhängige Veränder- 
liche ist u=wt. Dieses lineare System hat die konstante Partikularlösung 9, = — tg &,, wie 
man durch Einsetzen sofort sieht. Daher brauchen wir uns nur noch mit dem homogenen System 
zu befassen. . Dieses schreiben wir in der Gestalt &r’’—=h(w)Dr. Dabei ist & eine konstante, 
symmetrische und positiv definite Matrix. Letzteres zeigt man durch eine Umformung der 
Era: quadratischen Form. ® ist eine konstante indefinite Diagonalmatrix. Weiter ist 


h(u) = -; — pcos u. Wir transformieren dieses System mit Hilfe des Ansatzes = %9, wobei 3 


eine en Matrix mit nichtverschwindender Determinante sein soll. Dann ergibt sich 
C5Yy’=h(u)D%HYy. Da det 5-0, können wir diese Matrizengleichung von links mit der 
Transponierten 3° multiplizieren: $° CFy’=hu)FDTHY. Da C und D symmetrisch sind, 
außerdem € positiv definit, gibt eseine reguläre Matrix %, so daß gleichzeitig F’ C% die Einheits- 
matrix und $’D9OF—D* eine Diagonalmatrix wird, deren Elemente die Wurzeln der Gleichung 


det (e& —D)=0sind. Dann lautetalso das System: y/ (4 —pcos «) &%=0 (e=1...n). Nach 


Flo quet haben die Lösungen dieser Gleichungen die Form: Y,—= dı, EN fe(u) + de Ser l. (—u) 
mit f,(u) = f, (u-+ 2). Nur wenn alle w, reell und +5 (m —=0, 1, 2, ....) sind, können die 


1) Die ausführliche Arbeit wird im Journal für reine und angewandte Mathematik erscheinen. 
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Lösungen ) und damit auch x stabil sein. — Die oben angeführte Bedingung für die Mittellagen 


U>oo 


N - 4 AR 
ist äquivalent mit lim — 9. (7) dr=0. Da diese Grenzwerte über alle homogenen Lösungen 
0 


verschwinden, muß notwendig für die Existenz kleiner Schwingungen tg =0 (e=1...n) 
gelten. Jeder Pendelarm muß also um eine vertikale Lage schwingen. Daher gibt es 2% mögliche 
Mittellagen. Um nun für eine spezielle dieserMittellagen entscheiden zu können, ob sietatsächlich 
vorhanden ist, setzen wir in det (eC— D)—=0 die gewünschten &,-Werte ein und können nun in 


der Struttschen Karte feststellen, ob alle Wertepaare A, = 5% 1 Yu DENKEN n) 


in stabilen Bereichen liegen. Dann und nur dann, wenn dies der Fall ist, existieren kleine Schwin- 
gungen um die fragliche Mittellage. 

Zum Schluß soll ein zahlenmäßiges ‚Beispiel betrachtet werden. Wir wählen ein Doppel- 
pendel mit zwei zylindrischen Pendelarmen, 25cm lang, Radius lcm. Die Erschütterungs- 
amplitude sei 0,5 cm. 

Egli ee. (Das Pendel ist nach oben gestreckt.) 
Diese Lage ist stabil für alle Frequenzen 5 =»y>»,—=83 Hz. 

Lage 2: u =0, 8 =nr. Stabil für alle Frequenzen ' 9->v,— 56 Hz. 

Lage 3:9, =n, %&=0. Stabil für alle Frequenzen »>»,—=45Hz. 

Lage 4: sv = —=n. (Die natürliche Gleichgewichtslage.) Das Pendel ist in dieser Lage 
im allgemeinen stabil. Es gibt aber zwei Folgen von Instabilitäts- 
intervallen, deren Anfänge lauten: 


a NEL, 
0,85 2,28 
0,57 1652 


Für p=5cm sind die Lagen 1 und 2 nie stabil, die Lage 3 nur in einen schmalen Frequenz- 
intervall zwischen 5 und 10 Hz. Bei der Lage 4 treten ähnliche Verhältnisse wie oben auf. 


Verzögerungsfunktion 
Von E. Altenkirch in Neunhagen bei Berlin 
Die analytische Verfolgung eines mit Verzögerung auf den Fühler rückwirkenden Regel- 


vorganges führt zu der Funktion 
vr (Dry ya 


Vans 1) o+NI 1) ae 


die eine Folge endlicher Potenzreihen darstellt. a Untersuchung zeigt, daß bei al und 


bei a=> singuläre Punkte auftreten. Bei a<- ist die Funktion monoton konvergent, bei 


1 
3 er S 5 konvergiert sie unter Pendelungen, bei « >5 divergiert sie. 
(Die Variablen sind dimensionslos. x enthält die Zeit, & die Verzögerung.) 


a=5 liefert eine Sinusschwingung mit gedehnter Amplitude und Phasenverschiebung. 


Bei «=1 sind die Abweichungen vom Beharrungszustand um diesen gleichmäßig verteilt. 
Die Annäherung an ihn ist gegenüber dem verzögerungsfreien Zustand verlangsamt. 
Die «-Grenze gleicher Annäherungsgeschwindigkeit folgt aus 


y= = 1 — 077% 
für größere x zu &—0,66. Ist & kleiner, so ist die Annäherungsgeschwindigkeit stets Fgrößkr, 


1 
BTER gewährt die schnellste Annäherung, die um ein Vielfaches größer ist als bei 
verzögerungsfreiem Vorgang. 


Referate der Vorträge: 


H. Pailloux, Saarbrücken: Quelques asservissements pouvant se traiter par la 
methode de Lagrange, 


H. Bückner, Minden i. W.: Regelstabilität und Übergangsfunktion 
sind nicht eingegangen. 
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F. Elastomechanik 


Kinematographische Auflösung des 
Zerreißvorganges 


Von Hubert Schardin in Weil a. Rh. 


Der Standardversuch zur Bestimmung der Material- 
festigkeit ist das Zerreißen von Probestäben. 

Im allgemeinen Falle schnürt sich ein Zugstab zu- 
nächst infolge plastischer Deformation (d.h. durch Gleiten 
und Zwillingsbildung) ein, wird schließlich spröde und zer- 
reißt dann mit hoher Geschwindigkeit. Der Anteil des 
spröden Bruchs wird um so höher, je schneller die Be- 
lastung erfolgt. Amorphe Körper, bei denen keine neuen 
Gleichgewichtslagen existieren, sind von vornherein spröde. 

‘Die im folgenden kurz wiedergegebenen Untersuchun- 
gen beziehen sich nur auf den spröden Bruch. 

Ausgangspunkt waren Versuche an Glasplatten, die 
in den Glastechnischen Berichten!) wiedergegeben sind. 
Als Resultat derselben konnte festgestellt werden, daß nach 
kurzer Anlaufstrecke der Bruch mit der konstanten Aus- 
breitungsgeschwindigkeit von 1500 m/s sich durch das Glas 
bewegt. Mitunter findet eine Gabelung im Verlauf statt, 
der zunächst immer sehr genau senkrecht zur vorhandenen 
äußeren Zugspannung erfolgt. 

Über die Dynamik des Abbaus der Zugspannung durch 
die dem Bruch vorauslaufenden elastischen Wellen konnten 
eine Reihe von Aus- 
sagen gemacht werden, 
so u.a. daß in größerem 
Abstand vor der Bruch- 
front in dessen Ebene 
Blehfrische zunächst die ursprüng- 
Auslösung | liche Zugspannung auf- 
recht erhalten bleiben 
muß. 

Unmittelbar vor 
der Bruchfront war auf 
Belastungs- den Aufnahmen ein 
‚gewicht „Spannungshof‘ zu er- 
kennen, der eine Zu- 
nahme der Zugspan- 
nung bedeutet. Dieser 

Spannungshof ver- 
schwindet, wenn der 
Bruch stehen bleibt. 

Die dem Bruch 
vorauslaufenden elasti- 
schen Wellen konnten jedoch mit Glas als Versuchsmaterial 
nicht sichtbar gemacht werden. Daher wurden anschließend 
die im folgenden beschriebenen Versuche mit Plexiglas- 
platten durchgeführt. 


Die Anordnung ist in Bild 1 wiedergegeben. Die Ver- 
suchsplatte ist mit Hilfe von Araldit zwischen Metallbacken 
geklebt. Ein Belastungsgewicht fällt auf einen Teller, 
wodurch die Zugspannung in der Plexiglasplatte hervor- 
gerufen wird. 

Die Auslösung der funkenkinematographischen Appa- 
ratur erfolgt zu Beginn des Zerreißvorganges, indem un- 


Bild 1. Anordnung für den Zerreißversuch von 
Plexiglasplatten 


1) Ergebnisse der kinematographischen Untersuchungen des 
Glasbruchvorganges, Heft 1, 3 u. 12, Jahrgang 23 (1950). 


18 


Bild 2. 5 Teilbilder aus einer kinematographischen Serie. Der schwarze Streifen bildet den Kontakt zur Auslösung. Schliereneinstellung mit horizontaler Schlierenblende. 


Bildfrequenz im Original 120 000/s. Plexiglasstreifen 50 mm breit, 2 mm stark 
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mittelbar neben einer Einkerbung — zur Festlegung des Ortes des Bruchbeginns — ein auf- 
gedampfter Metallstreifen zerrissen wird. Die Aufnahmefrequenz der kinematographischen 
Apparatur beträgt im allgemeinen 100 000 Bilder in der Sekunde. Die Sichtbarmachung der 
elastischen Spannungen erfolgt teils schlieren- teils spannungsoptisch. 

Bild2 gibt eine Auswahl von fünf schlierenoptisch aufgenommenen Bildern einer kinemato- 
graphischen Serie wieder. Man erkennt auf den ersten drei Bildern den laufenden Bruch. Er 
bildet die Grenze zwischen einem oben hellen und unten dunklen Gebiet. Diese Helligkeits- 
verteilung beruht darauf, daß die freien Bruchflanken auseinandergehen und die vorher plan- 
parallele Platte an diesen Stellen keilförmig wird. Das Licht wird infolgedessen wie in zwei 


Prismen (Basis = Bruchfläche) abgelenkt. 


Vor der Bruchfront ist die Helligkeitsverteilung gerade umgekehrt, d.h. hier stoßen die 
„Prismen“ mit der Spitze in der Mitte zusammen: ein Beweis dafür, daß die Mitte vor der Bruch- 


front die höchsten Zugspannungen hat. 


Unmittelbar vor der Bruchfront ist ein schwarzes Gebiet, d.h. ein Gebiet mit hohen Span- 


nungen zu sehen, das sich auch in den spannungsoptischen Aufnahmen (Bild 3) wiederfindet. 


= ‚ 
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A Bild4. Nach vollständigem Durchriß der Plexiglas- 
platte entsteht, vom rechten Rand ausgehend, eine 

; Luftschallwelle. (Der Deutlichkeit in der Reproduk- 

R tion halber ist die ‚‚Tangente‘ an den ‚‚Kreis‘ durch 

er» \ Striche gekennzeichnet) 


form in die Platte hinein ausbreiten. 


beim Einlaufen des Bruches in den Rand plötzlich frei wird, so muß sich die Energie in Wellen- 


Besonders deutlich ist das auch auf den spannungsoptischen Aufnahmen zu sehen Bi) 


2 


Bild3. Spannungsoptische Bildserie, Plexiglasplatte 46 mm breit, 8mm stark. Bildfrequenz im Original 100 000/s 


Ferner sind nun auch auf diesen Aufnahmen die elastischen Wellen zu sehen, die dem 
‚Bruch vorauslaufen. Sie liegen gegeneinander exzentrisch. Nachdem der Bruch das Ende der 
Platte erreicht hat, starten von diesem Punkt am Rande neue sehr intensive Wellen. Der Span- 
nungshof vor der Bruchfront bedeutet eine Konzentration von Spannungsenergie und wenn diese 


sichtbare 
Wellen { 
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Bild5. Auswertung einer kinematographischen Auf- 
nahme des Zerreißversuches einer Plexiglasplatte 
(50 mm breit, 4 mm stark). Bildfrequenz rund 
100 000/s. 
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Die ersten vier Bilder dieser Serie zeigen den laufenden Bruch. Auf dem fünften Bild ist 
der Bruch gerade in den Rand eingelaufen. Und nun erkennt man auf den folgenden drei Bildern 
wieder die Ausbreitung der elastischen Spannung in etwas anderer Form. 

Diese Expansion ist sogar mit der Bildung einer Schallwelle in Luft verknüpft, die man 
auf den Schlierenaufnahmen erkennen kann (Bild 4). Die Luftwelle zeichnet sich ab als ein 
Kreis mit zwei Tangenten, die in den Rand der Platte dort einlaufen, wo sich die am Rande 
entlanglaufende Transversalwelle befindet. Dieser Punkt ist auf den spannungsoptischen Auf- 
nahmen besonders deutlich, aber auch auf den schlierenoptischen Bildern zu erkennen. 

Die vorliegende Form der Luftschallwelle ergibt sich, wenn man annimmt, daß nach dem 
vollständigen Durchriß der rechte Rand sich keilförmig aufweitet. Für diesen Fall ist die Bildung 
der Stoßwelle von M. J. Lighthill!) theoretisch ermittelt worden. 

Bild 5 bringt schließlich ein Beispiel für die Auswertung einer Aufnahme (Bild 4 hierzu 
gehörig). 

Der Bruch selbst breitet sich mit 560 m/s aus. In diesem Falle war jedoch der Bruchweg 
vorher mit kleinen Anbohrungen versehen worden, um stärkere Inhomogenitäten zu haben. 
Ist das nicht der Fall, so ist die Bruchgeschwindigkeit kleiner (Werte zwischen 305 und 460 m/s). 

Mit dem Bruch starten auch die elastischen Wellen, von denen die Lage der jeweils ersten 
sichtbaren eingetragen ist; die Weg-Zeitkurven der primären Longitudinal- und Transversal- 
welle ergeben sich auf Grund anderweitiger Messungen. 

Der Aufnahme unmittelbar sehr genau zu entnehmen ist jedoch die Transversalwelle, die 
am rechten Rande nach vollständigem Durchriß entlangläuft. 

Für die Ausbreitung der Schallwelle in Luft gilt in guter Übereinstimmung mit dem Ver- 
such die normale Schallgeschwindigkeit von 340 m/s. 

Die obigen Untersuchungen sind als experimentelle Grundlage für die Erforschung des 
Bruchvorganges gedacht. Theoretische Betrachtungen müssen immer von mehr oder weniger 
idealen Voraussetzungen ausgehen. Wie diese zu wählen sind, läßt sich leichter an Hand bereits 
vorliegenden Versuchsmaterials entscheiden. 


Zum Problem der Stabilität erzwungener Schwingungen 
elastischer Körper 


Von E. Mettler in Clausthal 


Die bisherige Theorie der Stabilität erzwungener Schwingungen elastischer Körper?) geht 
im allgemeinen von einer Grundbewegung aus, die mittels der klassischen Elastizitätstheorie 
berechnet wird, und untersucht das Verhalten der Nachbarbewegungen mit Hilfe der Methode 
der kleinen Schwingungen unter Berücksichtigung der höheren Glieder, die von der Elastizitäts- 
theorie endlicher Verschiebungen geliefert werden. An einem speziellen Beispiel der erzwungenen 
Querschwingungen des geraden, beiderseits in festen Gelenken gelagerten Stabes soll nun gezeigt 
werden, daß man zu weiteren Instabilitätsfällen kommt, wenn man schon die Grundbewegung 
mittels genauerer, höhere Glieder in den Verschiebungen enthaltender Gleichungen berechnet. 


- Es seien u und w die Verschiebungen der Se 
Stabachse in Längs- und Querrichtung. Die c EN 
x-Achse liege in der Stabachse. Tiefstehende EN 
Indizes x und £ bedeuten wie üblich partielle Be RE 
Ableitungen. Die Querbelastung wird speziell a 2 SS 

zu g= sin T cos wi(g, und w= const, NaRz, 
1— Stablänge) angenommen. Dann lauten die 
Differentialgleichungen für die Längs- und 
Querbewegung des Stabes mit den gewöhn- 
lichen Bezeichnungen 


2 3 
Br [ut tum=0 DEREN) 


“ 2 


1) Proc. Roy. Soc. 200, 8. 564 (1950). 
RN Val. Merle: Ing.-Arch. 16 (1947), S. 136 und 17 (1949), S. 418 und die dort angegebene Literatur. 


F : @. F. Carrier: Proc. Symposia appl. Math. 1 (1949), S. 175. 
; Zur Ableitung dieser Gleichungen im einzelnen vgl. K. Marguerre: Jb. dtsch. Luftfahrtforschg. 


(1938), 8.1433 und E. Mettler: Ing.-Arch. 13 (1942), 8. 97. 
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Dieses nichtlineare Gleichungssystem läßt sich näherungsweise lösen, wenn man die plausible 


2 
(und nachträglich kontrollierbare) Annahme einführt, daß die Dehnung etz der Stab- 


achse von x unabhängig ist. Dann erhält man nämlich durch Integration nach x über die Stab- 


länge und Berücksichtigung der Randbedingungen (3) die Längskraft 
l 
‚2 
S- Er (+) - Eau a NE 
{) 


und durch Einsetzen von (4) in (2) wird aus (2) eine Integrodifferentialgleichung für w allein?), 


die sich durch den Ansatz w= y(t) sin en integrieren läßt. Man findet für y die Gleichung 


l 
7 otyi1 + By?) En cos ot) we 
mit 
n\tEJ En 91 UM 
2 — | _—_ — = — u — ee 
(7) ey" P 4.,°.4308 p u 


Von Gl. (5), der Gleichung des Einmassenschwingers mit nichtlinearer Federkennlinie, ist die 
periodische Lösung 


y=Y,coswt-+Y3cos3wi+ --- . ONE A TH 
bekannt, wobei Y, aus Bild 1 abgelesen werden kann’). Der dort ram Parameter ist yo/2%, 


wo %, die klassische Durchbiegung unter der statischen Last q,sin T “ bedeutet. Ferner ist in 


guter Näherung®) 
I. 2% | =)| 
67 He 2) 7) - e IE 


Bleibt man in Bild 1 rechts von der punktierten Linie, so wird Y,/Y, <0,024. Es genügt 
also für unsre Zwecke, in (6) das erste Glied zu berücksichtigen. Zur Kontrolle der Lösung kann 
man nun wz/2 als bekannte Funktion in (1) einsetzen, diese Gleichung für wintegrieren und dann 
prüfen, ob die Annahme, daß e unabhängig von & sein soll, ausreichend erfüllt ist. Die Rechnung 
liefert 

s 
:—4[5) Y2[ 1-+cos2wt TER or ARE - 2 » 
s\ı j ( [3 ) ma 
Be EF 


‚und dieser Ausdruck ist in der Tat um so weniger mit x veränderlich, je tiefer &® unter dem 
2 
Längseigenfrequenzquadrat 7, = (7) = liegt. Für Stäbe mittlerer Schlankeit ist die Nähe- 


rung befriedigend, wenn & in die Nähe von o, fällt. 

Das Stabilitätsverhalten der Lösung von (5) ist bekannt. Die in dem schräg schraffierten 
Gebiet von Bild 1 liegenden Kurvenäste sind instabil. Die Bewegung zeigt die sog. Kipp- 
erscheinungen. Das läßt sich im einzelnen ableiten, indem man die berechnete Biegeschwingung 


als Grundbewegung betrachtet und daneben eine spezielle gestörte Bewegung annimmt, für die 


gleichfalls die Längskraft örtlich unveränderlich ist. Man kommt in einfacher Weise auf die 


übliche Mathieugleichung, welche die Nachbarbewegung beschreibt und deren erster Instabili- 


tätsbereich sich auf das in Bild 1 schraffierte Gebiet abbildet. 


4) Das gleiche Integral für S und die ER: Integrod:tferentialelasehneig für die freien Stab- Ti 
schwingungen wird (allerdings mit anderer Herleitung) benutzt von 8. FRE Aula J. Appl. Macz ar 


New York 17 (1950), 8. 35. 
5) Nach K. Marguerre: Neuere Festigkeitsprobleme des Ingenieurs. 
1950, S. 136. 


°) N. W. McLachlan: Ordinary non-linear differential equations in engineering an Physical sciences. I Sr RR 


Oxford 1950, 8. 54. 


2. angew. Math. Mech. 


es en SE A SS Sn a ae dl BL; 2 al REEL 2 
vr EEE, ”ıY N 


Bd.31 Nr. 8/9 Aug./Sept. 1951 F. Elastomechanik 265 


Affine oder projektive Kinematik, eine prinzipielle Untersuchung 
zu den Grundlagen der Rheologie 


Von Heinrich Hencky } in Mainz-Gustavsburg 


Bei der Anwendung des Erlanger Programms auf die Mechanik handelt es sich um drei 
Transformationsgruppen, die metrische Gruppe der starren Bewegungen, die alline Gruppe, 
welche die unendlich ferne Ebene in sich transformiert und die projektive Gruppe, welche eine 
Ebene im Endlichen, die Fluchtebene, mit der unendlich fernen Ebene zusammenfallen läßt. 

Vom mathematischen Gesichtspunkt aus gesehen, ist die Mechanik eine Differential- 
geometrie für ein noch von dem Zeitparameter abhängiges dreidimensionales Euklidisches 
Kontinuum. Dabei kommt es, wie S. Lie bereits gezeigt hat,nur auf die infinitesimalen Gruppen 
an, durch welche dann die Raumstruktur festgelegt wird. 

Für die Massenpunktmechanik genügt die affine Gruppe, welche zwar auch Verformungen 
"bewirkt, aber mathematisch am zweckmäßigsten ist, weil die Drehungen ja schon durch eine 
Spezialisierung der affinen Matrix bedingt sind. Man kann also auch starre Körperchen 
in ihrem kinematischen Verhalten mittels der affinen Technik transformieren, wobei man die 
notwendigen sechs Freiheitsgrade zur Verfügung hat. Wendet man aber die affine Technik auf 
das Element eines verformbaren Kontinuums an, so hat man in der infinitesimalen Trans- 
formationsmatrix nur die Ableitungen des Verschiebungsveklors zur Verfügung, was einen 
Verlust von drei Freiheitsgraden im Infinitesimalen bedeutet. 

Es ist aber absurd, eine Mechanik aufzubauen, bei der durch eine Erweichung der Starrheit 
ein Verlust von Freiheitsgraden erfolgt. Für den Physiker hat die Vorstellung eines Punkt- 
transports ohnehin etwas Abwegiges. Transportieren können wir nur ein kleines Teilchen, und 
dabei wird dieses Teilchen auch unwillkürlich gedreht und verdrückt. Als Grundlage der Ver- 
formung wird man daher eine Transformation wählen, welche gleich unsere Grundlfigur in der 
allgemeinst möglichen Weise verwandelt, und das kann nur die projektive Hauptgruppe. Unsere 
Grundfigur besteht aber aus dem infinitesimalen Koordinatenparallelepiped von acht Punkten. 
Nur durch die projektive Verformung können die acht Punkte des Parallelepipeds die allgemeinst 
mögliche Lage bekommen, die affine Gruppe würde nur ein schiefes Parallelepiped daraus 
machen. Die projektive infinitesimale Verformung ist nun gegeben durch den infinitesimalen 
Verschiebungsvektor und die Lage der Fluchtebene zum jeweiligen Aufpunkt. 

Der Unterschied einer affinen und einer projektiven Verformung des Kontinuums läßt sich 
in folgender Weise charakterisieren: 

Bei der affinen Verformung ist die Drehung des Teilchens durch die Ableitungen des 
Verschiebungsvektors schon mitbestimmt, die Feldrotation also mit der individuellen Drehung 
des Elements identisch. Bei der projektiven Verformung dagegen besteht im allgemeinen ein 
Unterschied zwischen der Feldrotation und der Teilchendrehung. Das gleiche gilt auch für die 
reinen Verformungen. 

Der Zusammenhang zwischen Feld und Teilchen ist also im Kontinuum des Euklidischen 
Raumes loser als im affinen Raum, in dem das Teilchen sofort drei seiner Freiheitsgrade verliert. 

Zur Gewinnung einer wirklichkeitsgemäßen Dynamik der Kontinua muß daher die bisher 
übliche affine Kinematik durch die Kinematik der Hauptgruppe ersetzt werden. 


Ein nichtkonserviertes Stabilitätsproblem 
Von Hans Ziegler in Zürich 


Ein fester elastischer Körper ist dann und nur dann stabil, wenn er sich unter dem Einfluß 
genügend kleiner Störungen beliebig wenig aus seiner Gleichgewichtslage entfernt, also eine 
beliebig schwache Schwingung um diese ausführt. 

Sind die inneren und äußeren Kräfte konservativ, so existiert eine potentielle Energie. 
Der Körper ist in diesem Falle so lange stabil, als diese positiv definit ist, und eine allfällig vor- 
handene kritische Belastung ist die kleinste Last, unter der dies nicht mehr zutrifft. ' 

j Da die potentielle Energie mit dem Auftreten nichttrivialer Gleichgewichtslagen aufhört, 
positiv definit zu sein, kann die kritische Belastung auch als diejenige ermittelt werden, unter 
der erstmals eine nichttriviale Gleichgewichtslage existiert. 

Von diesen drei Stabilitätskriterien — dem kinetischen, dem energetischen und dem 
‚statischen — sind die beiden letzten an die Voraussetzung gebunden, daß neben den inneren 
auch die äußeren Kräfte konservativ seien. Trifft dies nicht zu, so muß auf das kinetische 


_ Kriterium zurückgegriffen, d.h. eine Schwingungsaufgabe gelöst werden. 
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Es ist leicht einzusehen, daß das Kräftepaar mit konstantem Momentvektor nichtkonser- 
vativist. Das Stabilitätsproblem der dünnen Welle, welche durch axiale Endmomente belastet 
und damit auf Torsion beansprucht ist, ist daher — gleichgültig, ob man sich für das kritische 
Torsionsmoment oder für die kritische Drehzahl interessiert — im allgemeinen nichtkonservativ. 
In der Tat gelingt es nur in Ausnahmefällen, die Arbeit dieser Momente bei einer virtuellen 
Deformation als Variation eines Integrals darzustellen. 

Ist die Welle beidseitig eingespannt, so existiert ein solches Integral, und in der Tat ist 
hier das Moment konservativ. Die Anwendung des statischen Stabilitätskriteriums ist daher 


4 
noch gestattet, und dieses führt beim Knickproblem auf die Greenhillsche!) Formel W—= Be 


I 
in der W das kritische Torsionsmoment, «& die Biegesteifigkeit und } die Länge der Welle darstellt. 
Bei der einseitig eingespannten Welle versagen die beiden ersten Kriterien, da das Moment 
nichtkonservativ ist; insbesondere gibt es hier keine nichttriviale Gleichgewichtslage. Die 
Schwingungsrechnung liefert für das kritische Moment den unerwarteten Wert W=0°). 
Auch bei der beidseitig gelenkig gelagerten Welle ist das Moment nichtkonservativ, und 
die Schwingungsrechnung zeigt, daß ein kleines Torsionsmoment’auf die Querschwingungen ohne 


Einfluß ist. Es kann daher vermutet werden, daß hier die Greenhillsche Formel W= euer 


I 
samt ihren Erweiterungen, die R.Grammel?) zu verdanken sind, zu Recht bestehen. 

Läßt man die Welle gleichförmig um ihre Achse rotieren, so erhält man durch Zufügen der 
Trägheitskräfte aus der Knickaufgabe ein Drehzahlproblem, das im allgemeinen wieder als 
Schwingungsaufgabe behandelt werden muß. Dabei stellt sich beispielsweise im Falle der 
fliegenden Welle — wie auf Grund des Ergebnisses beim Knickproblem zu erwarten ist — heraus, 
daß für W=+0 jede Drehzahl kritisch ist. 

Die Resultate dieser Untersuchungen, die in der ZAMP ausführlich dargestellt werden 
sollen, decken sich — wenigstens was die Labilität der fliegenden Welle anbelangt — schlecht 
mit der Erfahrung. Sie beruhen aber auf der üblichen Voraussetzung einer ideal schlanken Welle 
und zudem auf der Vernachlässigung der Dämpfung. Da der Knickvorgang hier nicht wie im 
Eulerschen Falle in einem ebenen, monotonen Ausweichen der elastischen Linie besteht, 
sondern den Charakter eines Schwingungsvorgangs hat, bei dem sich die Stabelemente auf 
spiralförmigen Bahnen allmählich von ihrer Ausgangslage entfernen, ist nämlich auch mit einer 
Beeinflussung der Stabilität durch die Dämpfung zu rechnen. 

Jedenfalls zeigen die Ergebnisse, daß selbst ganz einfache Aufgaben nichtkonservativ sein 
und daß sie, falls sie als Stabilitätsprobleme auftreten, nicht mit den üblichen Stabilitätskriterien 
behandelt werden können. 


Zwei Bemerkungen zur allgemeinen Frage der klassischen 
Elastizitätstheorie 
Von @eorg Hamel in Landshut (Bayern) 


I. Z,, Z,, Z,von z unabhängig beiverschwindender Raumkraft verlangt nach den stalischen 
Gleichungen: 
Er 1. eine Airysche Spannungsfunktion A (x, y,2), so daß 


02A 02 A 02 A 
Kr 
2. eine Funktion B(x,y), so daß 
ee) BE 
5 °y ’ SE WRIN ’ 
5% BRAUN 


Indem man damit in die sechs Kompatibilitätsgleichungen hineingeht, bekommt man 
sechs Gleichungen für A, B, En die sich vollständig lösen lassen. Es werden 


1 1—o 65 
A=— rt gel et MH zr+zAt An 

1) A. G. Greenhill: On the strength of shafting when exposed ER to torsion and to end thrust. 
Proe. Inst. Mech. Eng. (1883), S. 182. 

2) Herr Grammel hat dieses Ergebnis vor einer Reihe von Jahren auf anderem Wege erhalten, aber 
nicht publiziert. ß 

®) Vgl. C.B. Biezeno und R. Grammel: Techn. Dynamik, Berlin 1939, S. 538if. 


‚nahme auf eine Festigkeitshypothese mit den Festigkeitskennwerten des Werkstoffes vergleicht. 


= ‚stoffen zu tun hat und auf den Vergleich mit dem Zugversuch angewiesen ist, bei dem nur eine 


u noch nicht das Verhalten des Werkstoffes in der Nähe der Fließgrenze berücksichtigt, das ja 
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wobei I, und % harmonisch in & und ysind, A,=«,, AA, —J, AA,= I, Z,=—k+ Ei Bi. 

und (1-+-0) B-- il, eine analytische Funktion ist. : 
Die zugehörigen Verschiebungen lassen sich ausrechnen, es ist dazu nützlich, eine allgemeine Be; 

biharinonische Funktion B= GE &F,-+-yF, einzuführen, wobei F,+sF, und G,-+?@, analy- 


tische Funktionen von SE sind, worauf dann 


F, RL! 
—=4(1-70)— +const wird und —|AB-+ —— (AA, +-const)|: F\- 
( ei FH Fk De ( 1 cons Ile Ta eg 
II. Die statischen Gleichungen lassen sich bekanntlich im dreidimensionalen Fall bei 
d2 Br 
p = Burn durch ar 
[old 7) 0° x ep 0y Rp. Op | 4 
X,=9+-—+- , Y,=8+-—— +, Z,=2+-H Be 
; fir re 92° R - 2 ö "ag da?’ 2 
0° 0° x Op 14 
A,=-- , S,=--, Me ag 
e 020 ä 9202 ® 9%y9% 
lösen. «+ 


Setzt man 9= 9,49, Yy=P%+93, X=9s3+-9,, So lassen sich aus der zweiten Zeile die 
Verschiebungen «, v, w berechnen und zwar zu 


(7) ya 2 
uu= — St +), m=— ln, U+h), oz Mn. is 9)» 


wo BE T2l. 0 d ne 2), A=h(z, Y) aaa En 
Dies sind die Integrale von 


ee 2) - ; (9-49) 
°Y 020Y 


Setzt man das in die drei ersten Gleichungen ein, so bekommt man drei lineare Ditierential- 
gleichungen zweiter Ordnung für 9,, 3, 9 z.B. 


USW. 


et 
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Wo 
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30 £ ray ; Team 


Über eine integrable, nicht-lineare Elastizitätstheorie 
von H. Neuber in Dresden 


Die technische Festigkeitslehre ist — von wenigen Ausnahmen abgesehen — bei der 
Berechnung der Spannungsverteilung in Bau- und Maschinenteilen seither aufdie lineare Elastizi- 
tätstheorie angewiesen, welche dem Hookeschen Gesetz entspricht. Die so erzielbaren Er- 
gebnisse haben zur Aufstellung von Formeln geführt, aus denen der Konstrukteur die an der 
Stelle der höchsten Beanspruchung auftretenden Spannungswerte errechnet und unter Bezug- 


Hierbei ergeben sich zunächst insofern Schwierigkeiten, als die genaue Berücksichtigung der 
Oberflächenform mathematisch nur in wenigen Fällen gelungen ist 1), Jedoch konnte durch 
Einschaltung geeigneter Interpolationsformeln Abhilfe geschaffen werden?). Bei schwacher Krüm- 
mung der Oberfläche schien ein solches Errechnungsverfahren bereits weitgehend brauchbar zu 
sein, bei stärkerer Oberflächenkrümmung jedoch machte sich ein zu starkes Anwachsen der 
theoretischen Spannungswerte bemerkbar, da die dem Hookeschen Gesetz entsprechende 
Elastizitätstheorie für Spitzkerben unendlich hohe Spannungswerte liefert. Wenn auch vom 
physikalischen Standpunkte aus im atomaren Bereich außerordentlich hohe Spannungen zu 
erwarten sind, ist ihre Ermittlung für den Ingenieur indiskutabel, da er es mit technischen Werk- 


mittlere Spannung des Querschnittes bestimmt wird. Durch Einführung eines endlichen Gefüge- 
teilchens konnte der Vf. diese Diskrepanz in erster Näherung überbrücken’). Damit war jedoch 


ERBE 
1) 2) a H. Neuber: Kerbspannungslehre 1937. 
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gerade bei Kerbwirkung in Betracht zu ziehen ist. Bei zähem Werkstoff ist in Wirklichkeit die 
Sicherheit gegenüber der Zerstörung weitaus größer, als die in der statischen Berechnung auf- 
tretende Sicherheitszahl erwarten läßt, welche auf Vergleichswerten beruht, die aus dem Hooke- 
schen Gesetz errechnet sind. Durch den Werkstoffmechanismus wird von selbst ein teilweiser 
Abbau der Spannungsspitzen und zugleich eine stärkere Belastung der ursprünglich geringer 
beanspruchten Fasern herbeigeführt (‚‚Schlauheit des Materials‘). Die Forderung nach stärkerer 
Ausnützung der Werkstoffe und die Verantwortung für die Sicherheit der Bauteile verlangt daher 
eine eingehende Erforschung des Beanspruchungsvorganges bei höheren Laststufen, insbesondere 
in der Nähe der Fließgrenze. Zur Klärung dieser Frage wurden zunächst umfangreiche mathe- 
matische Untersuchungen durchgeführt, welche zeigten, daß einer erfolgreichen Lösung dieses 
Problems erhebliche Schwierigkeiten gegenüberstehen. Für einen besonderen Fall zweidimen- 
sionaler Spannungsprobleme gelang dem Vf. jedoch die Aufstellung einer integrablen, nicht- 
linearen Elastizitätstheorie, welche zur Berücksichtigung des Einflusses der Oberflächenform 
einerseits und der Spannungsdehnungskurve zäher Werkstoffe andererseits geeignet ist und 
bereits zu fruchtbaren Ergebnissen geführt hat. Die Grundgedanken dieser Theorie lassen eine 
Analogie zu Problemen der Gasdynamik und der Flächentheorie erkennen. Es ergibt sich eine 
Differentialgleichung, deren Lösungen durch Überlagerung von zwei konformen Abbildungen 
darstellbar sind. 


Über eine Möglichkeit zur Kontrolle der physikalischen Widerspruchs- 
freiheit der Gleichungen der mathematischen Elastizitätstheorie!) 


Von @. Fichera in Rom 


Wir betrachten in einem r-dimensionalen euklidischen Raum S, ein berschänktes Gebiet D 
(offene Punktmenge). Der Einfachheit halber setzen wir voraus, daß die Grenze B von D aus 
einer einzigen geschlossenen Oberfläche der Klasse 2 besteht, d. h. lokal abbildbar ist mit Hilfe 
stetiger und mit stetigen partiellen Ableitungen erster und zweiter Ordnung versehener Funktionen. 
Sei u(z) ein Vektor mit’ r Komponenten, die wir mit n,(z), ug(%), ... , u,(&) bezeichnen. Mit 
© = (%, %3, . . ., %,) bezeichnen wir den in $, variablen Punkt. 


Wir setzen 
0? RM; ou; . m; 
I An» li — >; > e re ne). 
A: 2 AP div de, grad u, = 3 Ber 
Wir betrachten in D die folgende Vektorgleichung 
Asu-Kxrgrad IN US$ a ae 


x ist eine bestimmte Konstante und fein gegebener in D definierter Vektor. Wir werden auch 
den auf B definierten Differentialoperator betrachten, der den Vektor « in den Vektor L(w) 
mit den folgenden r Komponenten Bupzahr 


a (x — A) div um; + BE Re, 


Mi 


= 


L.[lu] = 


= (v1, 99, * , v,)ist der zu B normale Versor, der nach dem Innern von D gerichtet ist. d/dv be- 
deutet die Differentiierung in Richtung der inneren Normalen von D. A ist eine bestimmte 
Konstante. 

Wir zerlegen B in zwei offene Hyperflächen B, und B,. Das Problem, mit dem wir uns 
beschäftigen werden, betrifft die Integration der Gleichung (1), wobei der Vektor v auf B, und 
der Vektor L[u] auf B, vorgeschrieben ist. Im Spezialfall «= =0 erhält man das wohl- 
bekannte gemischte Pichler; von Dirichlet-Neumann für die Gleichung von La- 
place. Dieses Problem ist, wie wohlbekannt, für r = 2 äquivalent jenem, das darin besteht, 
eine in D holomorphe Funktion zu bestimmen, deren Realteil auf B, und deren Imaginärteil 
auf B, vorgeschrieben ist. Wenn man annimmt, daß r—=3, x >1,A= 1 befindet man sich im 


- Fall des fundamentalen Problems der mathematischen Theorie der Elastostatik, d.h. des Pro- 


blems des Gleichgewichts eines isotropen und homogenen elastischen Körpers, der auf einem 
Teil seines Randes vorgeschriebenen Bindungen unterworfen ist sowie der Wirkung eines 
bekannten Systems von Massenkräften in D und von Oberflächenkräften auf B, unterliegt. 
Ohne Beeinträchtigung der Allgemeinheit können wir v=0 auf B,, f=0 in D voraussetzen. 
Wir werden übrigens mit & den auf B, vorgeschriebenen Vektor L[w] bezeichnen. 


E) Vorliegende Arbeit wurde im Int ikmt Nazionale perle ABPi CH delCal- 


colo in Rom ausgeführt. 
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Um ein Eindeutigkeitstheorem zu sichern, werden wir voraussetzen, daß x > —1, 


— 1<A<y(x, r), wobei y(x, r) die kleinere der beiden Zahlen 1 und sa ist, oder auch 
en = beiden Fällen ist es möglich, in D einen Vektor /'[w] mit #, Komponenten 
"lu], Fall, **, 7/,.[w] zu definieren, die Linearkombinationen der ersten SE ALEEN der 


Komponenten von uw sind, in der Weise, daß die folgende Formel besteht 
fu . Liv]} do + tel . T'[o]} de=0. 
B 


u ist ein Vektor, der der Klasse U aller Vektoren angehört, die stetige und mit in D integrier- 
‘ barer Norm versehene erste Ableitungen besitzen. v gehört der Klasse jener Vektoren an, die mit 
ihren ersten Ableitungen in D-- B stetig sind und in D der Gleichung genügen 


LES Bar DU 1 Ga a SB Be ER AL} 


Wir werden weiterhin die Unterklasse U, von U betrachten, die aus den Vektoren von U be- 
steht, die auf B, verschwinden und Lösungen von (1°) sind. 

Es ist möglich, ein Existenz- und Eindeutigkeitstheorem in der Klasse U, für das gegebene 
Problem aufzustellen ?). 

Der erste Schritt des Beweises besteht darin, zu zeigen, daß die Klasse U, angesehen 
werden kann als ein linearer, normaler und vollständiger Funktionenraum im Sinne der all- 
gemeinen Analysis. Das wird möglich, indem man die Funktionalnorm eines Vektors «in U, in 


folgender Weise definiert | 
all = | J Tel? dep. 


Der zweite Schritt besteht im Beweis der Existenz in U, eines orthonormalen und vollständigen 
Systems {vo}, wobei v') ein in D-+ B mit stetigen ersten Ableitungen versehener Vektor ist. 
Jeder Vektor « von U, kann dann in folgender Weise in eine im Innern von D gleichmäßig 
konvergente Reihe entwickelt werden 


oo 


uam, Bl tBub Diego 2'022 002 70.02 
D 


n=0 
Die Eindeutigkeit der Lösung in U, folgt leicht aus der Bemerkung, daß 
6, el {T [u] . I{w®]} de= — if (um) 5 8) d<=0 i 
B . 


wenn 2 =(0 auf B,. 
Wenn umgekehrt der Vektor % gegeben ist und 


| Er (v9 .» 2) dx 


gesetzt wird, beweist man unter der mehr als hinreichenden Bedingung, daß & auf B, stetig ist 
und hier eine Höldersche Bedingung erfüllt, daß die Reihe (2) die Lösung des Problems 
definiert, und in dieser Weise gelangt man zum Existenztheorem. 

Das Existenz- und Eindeutigkeitstheorem, das wir skizziert haben, sichert, daß die er- 
haltene Lösung folgende Bedingung erfüllt 


[rw aa<+ ©. 


Im Falle der Elastizität (= 3, «>1, A=1) wird diese Bedingung von der physikalischen 
Natur des Problems erfordert. Sie drückt in der Tat aus, daß die Deformationsenergie endlich 
ist. Aber von diesem Standpunkt wird weiterhin die Stetigkeit des Vektors u und seiner ersten 
Ableitungen in D-+ B erfordert. Es ist möglich, anhand von Beispielen zu zeigen, daß die 
bloße Stetigkeit oder Differenzierbarkeit beliebig hoher Ordnung des auf B, gegebenen % nicht 
die Stetigkeit der Ableitungen v in D-+ B nach sich zieht. Es ergibt sich daher die Notwendig- 
keit einer weiteren Untersuchung über die Bedingungen, denen 2 genügen muß, um diese Stetig- 
keit zu sichern. 

An einem Spezialfall werden wir den überraschenden Umstand beleuchten, nach dem für 
die Stetigkeit in D-+ B.der ersten Ableitungen von « notwendig ist, daß der gegebene Vektor & 

eine quantitative Bedingung folgenden Typs erfüllt 


je 90-0 a | 08; 
wobei r ein Vektor ist, der ausschließlich vom Gebiet D abhängt. Setzen wirr=23,x=/—0, 
d.h. betrachten wir den Fall des gemischten Problems von Dirichlet-Neumann für 


2) Vgl. G. Fichera, Sull’estisenta e sul calcolo delle soluzioni dei problemi al contorno relativi 
 all’equilibrio di un corpo elastico. Ann, della Sinola Normeli di Pisa 1950. 
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die harmonischen Funktionen in zwei unabhängigen Veränderlichen. Die Grenze B ist in zwei 
offene Kurven B, und B, zerlegt, die als gemeinsame Endpunkte z®U und x haben. Man be- 
trachtet die beschränkte harmonische Funktion », die auf B, den Wert z/2 und auf B, den 
Wert 0 annimmt. Es ist bekannt, daß eine solche Funktion in folgender Weise zerlegbar ist: 


va) =) — 5 9la)+ 5 Aa). 


[4 


9, (x) ist der Winkel (zwischen O0 und 27), den z, mit der posigiven Tangentenachse an Bin z, 
bildet, und in analoger Weise wird ©, (x) definiert; (2) ist eine harmonische und stetige Funktion 
in D-+ B, und man stellt fest, daß ihre Konjungierte & eine stetige Funktion in D+ B ist. 


Wir setzen 


1 1 
uU=W— — > 108 Iu— «| Sr log |, — |, 


— „42 

(2) = e* 05 9 = = ein e” cos®, 
RAN 

B(z) = e®sen #—= m, = A erisenv. 


Sei C eine einfache reguläre geschlossene Kurve im Innern von D. Man hat 


du da - ((« du Eu) as = 


Setzt man die ersten Ableitungen von « als stetigin D-+ B voraus und läßt man die Kurve © 
gegen B streben, so erhält man 


Lo — 12 | 
(nn) er 2-0. ee a er RE . (4). 
2 $ 
2.B. wenn B de is a < © < 2 ist, wobei 
in diesem Fall ® = 0 wird, hat man 
7 %, — a||112 ’ 
| sie) doo A 


Im allgemeinen Fall kann man berechtigterweise erwarten, daß eine Bedingung des Typs (3) 
besteht. In der Tat ist es möglich, die theoretische Existenz eines Vektors a zu beweisen, dessen, 
Komponenten r additive Funktionen mit endlicher Variation der Borelschen Mengen von B, 
sind, so daß für die Stetigkeit in D-+ B der Derivierten von % die folgende Bedingung not- 


- wendig ist: 


> do = On ee 
BEL. (5) 


Das Integral in (5) ist ein Integral im Sinne von Stieltjes. Es ist ein Problem der mathe- 
natischen Physik, die Bedingungen (4), (4°), (5) physikalisch zu interpretieren. ; 
Wenn man beweisen kann, daß derartige Bedingungen keinerlei physikalischen Sinn be- 
sitzen, muß man schließen, daß die mathematische Schematisierung jener physikalischen Pro- 
bleme, auf welche sich jene Bedingungen beziehen, — insbesondere des Problems des elastischen 
Gleichgewichts — illusorisch sind. i 


Zum Elastizitätsproblem des kreiszylindrischen Flüssigkeitsbehälters 
veränderlicher Wandstärke 


Von Karl Federhofer in Graz 


Die Lösung des Elastizitätsproblems des kreiszylindrischen Flüssigkeitsbehälters mit stetig 
veränderlicher Wandstärke erfordert die Integration einer linearen Differentialgleichung 4. Ord- 
nung für die radiale Verschiebung der Punkte der Mittelfläche ; sie besitzt veränderliche Koeffi- 
zienten, die durch das Gesetz der Wandstärkenänderung bestimmt sind. Geschlossene Lösungen 
durch bekannte Funktionen sind für konstante und linear veränderliche Wanddicke 6 bekannt; 
es scheint bisher nicht bemerkt worden zu sein, daß im Falle einer quadratischen Veränderlichkeit 


& 
s 
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von ö das Problem durch elementare Funktionen gelöst werden kann.!) Wählt man nämlich 
als Grundveränderliche die Schubkraft Q und die Veränderliche V=y ö2, wo x den Biegewinkel 
_ der Zylindererzeugenden bedeutet, so genügen Q und V mit Einführung des Operators 
AIR E 


den Differentialgleichungen 
E ) 2 (1—»% 
LO)+ v-(p—%9)5, un—2uy-ZLNg 


(p=yx&= Druck der Füllflüssigkeit in der Tiefe x unter dem oberen Rande.) Hochgestellte 
Zahlen bedeuten Ableitungen nach x. Das homogene Gleichungssystem für Q und V lautet 


dann: 
Lug zldg-o),  z1um-esn+ ld yo. 
Die Aufspaltung der a. 4. Ordnung in solche 2. Ordnung von der Kor 
LQ)+cQ=0 [BE COUSUW re Eh) 
liefert die Aufspaltungsbedingung 
12 (1—»?) 
@+2ö1c ————=(, 


wonach ö!! gleich sein muß einer beliebig wählbaren Konstanten k. Aus öl=% folgt für das 
Wandstärkengesetz 


= artbutb, RE BE RUE I, 


Die Konstanten %,b, b, sind bestimmt durch die Werte von ö am oberen und unteren Behälter- 
rande und durch die noch freie Wahl der Lage der zur Zylinderachse parallelen Scheiteltangente 
der nach (2) parabolischen Wandbegrenzung. 

Mit dem Wandstärkengesetz (2), das die Fälle konstanter und linear veränderlicher Wand- 
stärke als Sonderfälle enthält, läßt sich (1) bei beliebiger Lage der Scheiteltangente überführen 
in eine hypergeometrische Differentialgleichung. Verlegt man aber die Scheiteltangente der 
Wandbegrenzung in die Mittelfläche der Zylinderschale, was bei der Annahme 2%b,—b} zutrifft, 
dann vereinfacht sich die Gl. (1) in eine gewöhnliche Eulersche Gleichung, die durch einfache 
Potenzfunktionen von &in geschlossener Form integriert werden kann, so daß auch Biegungs- 
moment, Querkraft und Formänderung des Behälters durch geschlossene Ausdrücke ebenso 
einfach wie im Falle konstanter Wandstärke dargestellt werden können. 

Eine ausführliche Darstellung ist kürzlich erschienen im Anzeiger der math.-naturw. 
Klasse der Österr. Akad. d. Wiss., Jg. 1950, S. 275—287. Über numerische Ergebnisse für 
Behälter verschiedener charakteristischer Kennzahlen soll im Bd. V (1951) des Österr. Ing.-Arch. 
berichtet werden. 


Die mehrfach durchlochte Platte 
Von 0. B. Biezeno in Delft 


Einer Aufforderung der Tagungsleitung befolgend gibt der Redner eine Übersicht über die 
letzten Arbeiten seines während des Krieges ums Leben gebrachten Mitarbeiters Dr.ir.K.J.Schulz. 
Sie beziehen sich auf die durchlochte Platte und behandeln in einheitlicher Weise, und zwar 
mittels eines Iterationsverfahrens eine ganze Reihe alter und neuer Probleme, welche unter 
Verweis auf die Originalabhandlungen hier nur in aller Kürze aufgezählt werden: die einfach 
durchlochte Platte, die zweifach durchlochte Platte, die Platte mit einer unendlichen Reihe 
kollinearer Löcher, die Platte mit drei nicht kollinearen Löchern, die Platte mit zwei parallelen 
(unendlichen) Reihen von Löchern (welche auch gegeneinander verschoben sein können), Platten, 
bei denen die Lochränder in vorgeschriebener Weise belastet sind, die Halbplatte mit einer un- 
endlichen Reihe von Löchern parallel zu ihrem freien Rand, der unendliche Streifen endlicher 
Breite mit einer unendlichen Reihe von Löchern, entweder auf einfachem Zug oder auf reiner 
Biegung beansprucht. 

Siehe die Originalarbeiten: 


'_K.J.Schulz: Over de spanningstoestand in N oorde platen. Diss. 1941, Delft; On the State of 
Stress in Perforated Plates, Proc. Ban Ak. v. Wet. A’dam XLV, (1942), S. 232—239, 341—346, 457 —464, 
524—532; ferner: XLVIII (1945), S. 2832—291 und 292—298. 


1) Zu obiger Feststellung, und zwar erweitert auf die anisotrope Zylinderschale, gelangte bereits, wie 
ich erst nachträglich bemerkte, E. Steuermann: Z. angew. Math. Mech. 5 (1925), S. 463 auf Grund einer 
von der hier benutzten Meißnerschen Methode abweichenden Beweisführung. Eine bis zur numerischen 
Durchrechnung reichende Lösung des Problems des Flüssigkeitsbehälters mit quadratisch veränderlicher 

‘ Wandstärke ist dagegen in der erwähnten Arbeit nicht gegeben. 
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Zur Spannungsfunktion bei Membranschalen 
Von W. Zerna in Hannover 


Die Schalenmittelfläche # möge durch allgemeine Flächenkoordinaten 9,, d, bestimmt sein. 
Mittels einer Parallelprojektion werde das auf F gelegene: Koordinatennetz auf eine zu den Pro- 
jektionsstrahlen senkrecht gelegenen Ebene E projiziert gedacht, so daß das Kurvennetz auf F 
als Bild auf E& erscheint. Jedem Punkt auf F ist damit ein Punkt auf E zugeordnet. Es werde 
vorausgesetzt, daß diese Zuordnung ein-eindeutig ist. Die auf den Projektionsstrahlen gemessene 
Entfernung der einander entsprechenden Pinikte sei’2— 2(9,,3,)- 

Das Gleichgewicht des in Z gelegenen Bildes eines Flächenelementes läßt sich ausdrücken 


durch!) 
nr, + PP EINER RR WI DE 


Darin ist pP durch die Schalenbelastung gegeben und n°? ist ein in E gelegener kontravarianter 
Tensor, der mit dem in F gelegenen Längskrafttensor N*# verbunden ist durch 


n?P— \% N&#B", 


worin a und A die Determinanten des Mater in F und #& sind. 
Da n*P ein Tensor des euklidischen Raumes ist, sind die höheren kovarianten Ableitungen 
vertauschbar und (1) läßt sich durch den Ansatz 


neB— ern Ol,» — 4°? STR Er 2 Re A ER: (2) 


identisch befriedigen. In (2) stellt © eine Spannungsfunktion dar. 
IR, 1 { | 
Heel eat, und. .AeR, =, wobei A en 


gewählt werden kann. 
Zur Ermittlung der Spannungsfunktion steht die dritte Gleichgewichtsbedingung der 
Schalentheorie zur Verfügung, die in diesem Fall 


+0 EB LEN 


lautet, worin p? der Lastkomponente senkrecht zu E entspricht, und 5,3; sind die Koeffizienten 
der zweiten Fundamentalform von F. 
‚Einführung von (2) und (3) ergibt die Differentialgleichung für die Spannungsfunktion, 
sie lautet | 
ERET O RA 
mit - 


5 
hr VE ph). 


Ist das Koordinatensystem in E kartesisch, so gehen die kovarianten Ableitungen in (4) in 
die gewöhnlichen partiellen über. Die sich dann ergebende Gleichung ist in anderer Form und 
auf anderem Wege von Pucher?) gefunden worden. Hier zeigt sich die Überlegenheit des 
Tensorkalküls, der es gestattet, die allgemeinen Beziehungen sofort mühelos zu gewinnen. 


Berechnung der Druckverteilung eines belasteten Balkens 
auf elastischem Untergrund 


Von Th. Schade in Aachen 


1. Einleitung 


Die vorliegende Methode weist einen Weg, das Problem des belasteten Balkens auf einem 
elastischen Untergrund mit allen Randbedingungen und physiaklischen Bedingungen näherungs- 
weise zu lösen. Dabei wird ein „‚Bettungsmodul‘“ definiert, der sich aus physikalischen und geome- 
trisch bekannten Größen zusammensetzt und nicht durch Versuche ermittelt werden muß. Die 
Formulierung des Problems als lineare Integralgleichung gestattet dabei die Anwendung des 
Superpositionsprinzips der Fertigkeitslehre. 


!) Die Indizes durchlaufen die Zahlen 1, 2, über doppelt auftretende gleiche Indizes ist zu a 
der senkrechte Strich zeigt kovariante Differentiation an. an 
®2) A. Pucher: Beton u. Eisen 33 (1934), S. 298. 
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2. Analytische Grundlagen 


Die Integralgleichung, die sich aus der Bedingung gleicher Durchbiegung von Balken und 
Untergrund in y-Richtung ergibt, lautet in dimensionslosen Koordinaten 2=a£. 


0 


Kar RT 
PURE ; 2 3 
J pw) InlE — vldv— — ST ja [as[ae| (1E)-P®) deto+ +05 +5 A). 


Wenn E, der Elastizitätsmodul des Balkens, E7 der Elastizitätsmodul des Untergrundes, I das 
Trägheitsmoment des Balkens, b seine Breite und 2 a seine Länge bedeutet; q (£) bzw. p(£) sind 
als Auflastfunktion bzw. Bodendruckfunktion 
in [kg/cm?] angegeben. Die Integrationskon- 


— oo ur ELLI 
a ‚P (&) (re — 
stanten [02 und [23 werden aus der gegebenen 


p(&) 


Querkraft und Momentenbedingung an den ; 
Balkenenden bestimmt. ce, ergibt sich aus dem 
Wert der ersten Ableitung in Balkenmitte. 

6, bleibt unbestimmt. 


a — 


[) 


Bild 2. Bodendruckverteilung für gleichmäßige Auflastverteilung 
(Fall 1) 


Bild 3, Bodendruckverteilung für eine Belastung durch eine Bild 4. Bodendruckverteilung für eine Belastung von symmetrisch 
Binzellast in Balkenmitte (Fall 3) angeordnet und gleichen Einzelkräften für = 0,707a (Fall 4) 


Setzt man &= — cos 9, 9—= — cos &, schreibt für die rechte Seite der Integralgleichung 


*. 9(— cos p) und setzt für den Faktor der rechten Seite, den „‚Bettungsmodul“ 
J N - x E Ay ba® Ey 9 
i\ 9 S: ARagaRE u EEE 
a so heißt die Hamelsche Lösung der Integralgleichung 
WER NT ER / TR RER /R IS 
2 = — —— [| |) + — 0 Sr ub,cosuPp » 2: Hr...) 
r 70) el) De np ( 
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mit den Fourierkoeffizienten 
7 Aa} 


Rrie, 
Be g(— cos p)dp, b,= Vasen nesura EN 
Yr, ie 
0 
Mit dem Lösungssatz 
1 oo 
np Gy COS PD 2 ee Er 


werden die b, Linearkombinationen der a,. Setzt man (5) in (3) ein, so kann man mit der Gleich- 
gewichtsbedingung 


era 3 KICL: ae Re RE (6) 


aynn2 9% 
und der Spannungsbedingung an den Balkenenden, 


HE Pe a cn Fi EEE 


v=0 


wenn. man bei einemgewissen u — N abbricht, mit derMethode der kleinsten Quadrate das damit 
entstehende Gleichungssystem für die Koeffizienten der Druckfunktion lösen. In den folgenden 
Bildern (Bild 2—4) sind die Druckverteilungen für einige Belastungsfälle (N = 10) ausgerechnet. 


Dabei ist = 
7T IT 

ee SE N ER ER Do ee a N 

SER y3 Ss 


Es sei noch bemerkt, daß, wie man aus (1) erkennt, folgende Grenzkurven gelten: 


Ay 


ie 


a>0 p>g, REDE 


Zum Problem der schiefen Platte 
Von W. Swida in Karlsruhe £ 


Das Problem der schiefen (parallelogrammförmigen) Platte hat eine praktische Bedeutung 
für den Stahlbetonbrückenbau, weil die meisten Brücken schief sind. Die Lösung dieser Aufgabe 
in schiefwinkligen Koordinaten ist in den Arbeiten von Favre!), Lardy?) und Nielsen?) 
angegeben worden. 

Im folgenden werden die Grundlagen eines Näherungsverfahrens entwickelt, bei dem die 
Platte durch ein unsymmetrisches hexagonales Gewebe ersetzt wird. 

Die Gleichgewichtsbedingung für einen beliebigen Maschenpunkt % lautet (Bild 1): 


” tg (0, — 20, + w,) +18 Bl — 20; + w) Hg y (m — 2m + w,) + Pie »=0 (1) 


wobei 


4 |t ts B+ "SP (sin2&+sin2y) 


ist und w,, w;, w; usw. die Durchbiegungen des Gewebes in den Maschenpunkten sind. Im Bild 1 
sind R,, R,, R, usw. die Spannkräfte der Ersatzdrähte. Weiterhin ist 8, die Beanspruchung der 
elastischen Haut je Längeneinheit. 

Die siebengliedrigen Differenzengleichungen nach der Art der Gl. (1) müssen für jeden 
Punkt des Maschennetzes angeschrieben werden. Die Auflösung der Gleichungen liefert die 
Ordinaten der Verformungsfläche des elastischen Gewebes. 

Wenn man nun die Verformungsfläche des elastischen Gewebes als Belastungsfläche für 
ein neues zweites Gewebe betrachtet, wenn man also die elastischen Gewichte p,—=w; an den 


1) H. Favre: Schweiz. Bauztg. 120 (1942), "8. 35. Abhandl. der Intern, Vereinigung für Brückenbau 
und Hochbau. Zürich 1943/44, S. 91. 

®) P.Lardy: Schweiz. Bauztg. 67 (1949), Nr. 15. 

s) N.I. Nielsen: Skaevvinklede Plader, Kovenhavn, Danmark 1944. Akademiet for de Tekniske 
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Knotenpunkten eines Gewebes mit der spezifischen Beanspruchung S, anbringt!), so erhält man 


folgendes Gleichungssystem für die sieben benachbarten Punkte, i, k, 1, o, p, 9, r (Bild 2) 2 
ala —m— 2) HE Paa—n— 2) Hg y@a nn), | 
N 2 


EBENE NE ES Zune ae a N Rn et en en er aan 


= 
Di \ 
Bild 1 Bild 2 


w-Werte aus diesen Gleichungen und setzt man sie in die G]. (1), so erhält man folgende 19-gliedrige 
Differenzengleichung: 


2%" N 52% (2) Ar 24) Ö, er (2; es 2) Ö, Br (2o % 2,) Ö4 m 2 tg & tg B(z +24) > 2 tg & tg Y(Zm ar | / 
Hat Bigyat a) +++ Ban ta) He + = F 


SS 
,—=Altgartgß-+tgy)+2(terattefttery), 
,—=A4tgaligattgßttgy)—2tgPß-tgy, 
‘ = 4tgßligattsß+tgy)— 2tga-tgy, 
r 4, —=4dtgy(ligattsßrtgy)— 2tga-tgß. | 
) Aus (2) kann als Sonderfall die Gleichung von Göttlicher?) für ein symmetrisches hexa- 
gonales Gewebe erhalten werden. . 


mit 


Die Ordinaten der Biegefläche der Platte werden &, = 2;: en ®. 
Die Ordinaten w und z verschwinden an den Rändern einer frei aufgelagerten Platte. Die 


Erfüllung der Randbedingungen für eine festeingespannte Platte bereitet einige Schwierigkeiten. ‘ 


Über die physikalische Bedeutung der Mohrschen Hüllkurve 
Von ©. Torre in Wien 


Otto Möhr[1] hat die Hüllkurve aller Grenzspannungskreise aus rein gedanklichen Er- 
wägungen aufgestellt und ihren Verlauf als eineglatte Kurve schematisch angegeben. Ausseiner 
Zeichnung läßt sich entnehmen, daß die Hüllkurve stetig und stetig differenzierbar sein muß. 
Erste Versuche zur Klärung der physikalischen Bedeutung der Mohrschen Hüllkurve haben 
die Schüler von Ludwig Prandtl: Th. v. Kärmän[2] und Robert Böker [3] durchgeführt. 
N Sie haben u. a. die Übereinstimmung der wirklichen mit den aus der Hüllkurve ermittelten Bruch- 

% winkel untersucht. Es ist nun die Normal- und Schubspannung o und z in der Bruchfläche aus 


1) H. Marcus: Theorie des elastischen Gewebes, 1924, S. 17. 
2%) Göttlicher: Ing.-Arch. 9 (1938), S. 12 (s. auch H. Vogt: Beitrag zur Berechnung schiefwinkliger 


EI, 4 


re Platten. Dissertation, Techn. Hochschule Hannover, 1938). 
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den versuchsmäßig ermittelten, an der Bruchgrenze liegenden Hauptnormalspannungen 9], 03 
und dem Bruchwinkel & leicht zu berechnen. Sie lauten (aus o=o, c0os?0-+-o, sin? a, wenn 


sinus und cosinus durch tangens ausgedrückt werden) 


o-+ o3tg8?& 
1+tg2a ° 


Oy — 


tg en Kl 


Der Index ‚‚w‘‘ bedeutet die wahren Spannungen in der Bruchfläche gegenüber jenen, die sich 
aus der Mohrschen Hüllkurve (Index ‚,h‘“) ergeben [4], [5]. \ 


+90 


1+a ° 


0, — 


0] 03 | ’ 
——— ° * * . . ” * * * ” 2 [3 
177 = 1 01 Yoi ( ) 


Dabei ist 0; = do,/do, die erste Ableitung der Grenzkurve 


> 0, @s)\...:» order in Be a in 
Gl. (2) ergibt sich aus einer Berührungstransformation der Gl. (3). Mit Gl. (1), (2) sind die 


Kurven r—=r(o) in Parameterform ausgedrückt. 


(G 


0, 20/reell) 


Bild 1 


Gl. (4) besagt, daß die reelle Berührung der Grenzspannungskreise mit der Hüllkurve 4 > 0 


’ 
| (imaginär) 
6;:0 /Krummungs- 
‚kreis ım Scheitel) 


Wie man sieht, die Hüllkurve [Gl. (2)] ergibt 

sich aus einer rein formalen Transformation 
der Gl.(3), so daß Gl.(2) zwar mit der 
Gl. (1) übereinstimmen kann, dies aber 
nicht muß. Das würde dann und nur dann 
zutreffen, wenn laut Gl. (1) und (2) 


{gr ao 9 Rue Den, 


Bild 2 


erfordert. Da erfahrungsgemäß die Hüllkurve der o-Achse die konkave Seite ständig zuwenden 


muß, so muß 


o<0O 


ständig sein. Dann hat die Hüllkurve keinen Wendepunkt (Bild 1). Das hat auch Mohr 
vermutet und das mit einer Skizze zum Ausdruck gebracht. Die Gleichung des Krümmungs- 


halbmessers der Hüllkurve lautet 


Pr ld DEE 


2 


[4 


2) 


Ita, 1 BEE VERS ee PETE 


020 


"Da dieser bei reeller Berührung größer als die Normale der Hüllkurve (0, —03)/2 ist, so muß 


das zweite Glied der G1.(6) positiv sein. Das ist wegen o} > 0 und oj < 0 auch tatsächlich der 


Fall. Im Scheitel der Hüllkurve gilt n—=0 bzw. R=(o,—0;)/2; hier ist der Krümmungskreis 
gleich dem Spannungskreis. Ist 04 < 0 (beioi < 0), so wird R< (o,—03)/2. Innerhalb des Krüm- 
mungskreises im Scheitel der Hüllkurve treten hier zum erstenmal die die Hüllkurve imagi- 
när berührenden Grenzspannungskreise, die wahrscheinlich alle den Trennbruch 
ergeben, auft). Sie schrumpfen in einem Punkt D zusammen (Bild 1), der die allseitig gleiche 
Zugfestigkeit, d.h. die Kohäsionsfestigkeit des Materials ergibt. Er ist ein isolierter 


Punkt der Hüllkurve und wahrscheinlich der A 


für alle Belastungsarten (Bild 1). 


Zwecks Überprüfung der Übereinstimmung dieser Gleichungen haben wir die Versuchs- 


usgangspunkt aller wahren Grenzkurven [G1.(2)] 


ergebnisse von Ros und Eichinger [6], [7] im K’S(o,r) aufgetragen (s. [4], Bild 14 oder [5], 


Bild 12). Es ergibt sich, daß die Ergebnisse der Gl. (1) und (2) nicht identisch sind. Im Druck- 
gebiet stellt die Hüllkurve eine gute Näherung dar, während sich im Zuggebiet beträchtliche 
Unterschiede von der wahren Grenzkurve [Gl]. (2)] ergeben. Hier geht diese Grenzkurve l- 


nicht gleich groß sind. 
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mählich vom Verschiebungs- zum Trennbruch über, wobei der Krümmungskreis im Scheitel 
der Hüllkurve jener Spannungskreis ist, der als erster den Trennbruch liefert (nach Leon [8]). 
Weitere Trennbrüche ergeben sich — wie schon erwähnt — aus den der Hüllkurve imaginär 
berührenden Grenzspannungskreise. 

Aus den Versuchen wissen wir, daß die gemessenen Bruchwinkel mit den theoretischen 
nach der Druckbelastungsart besser, mit jenen nach der Zugbelastungsart schlechter überein- 
stimmen. Um dies zu erklären, müssen wir uns den Schnitt der Grenzfläche F(o;, 05, 0,)—=0 mit 
der Ebene ,=r0,—+-.003, r-+-g=]1, räumlich vorstellen. Die Projektion dieser Schnittkurve in 
die Ebene (03, 0,) liefert die Grenzkurve Gl. (3), aus der sich durch Transformation die Hüllkurve 
[Gl. (2)] ergibt (s. Bild 2, 6, 7 des Verfassers [5]). Je nachdem, ob die Ebene o,—=ro,—+-g90, mehr 
oder weniger schief zur Projektionsebene (0,,0,) steht, werden die Grenzkurven [Gl. (3)] mehr 
oder weniger verzerrt ausfallen. Ebenso werden auch die Hüllkurven mehr oder weniger verzerrt 
sein. Diese Verzerrungen verursachen die besseren bzw. schlechteren Übereinstimmungen der 
wirklichen gegenüber den theoretischen Bruchwinkeln, d.h. jenen Bruchwinkeln und Spannungen, 
die sich aus der Hüllkurve [Gl]. (4), (2)] ergeben ?). 

Der isolierte Punkt D fällt mit dem Scheitelpunkt A der Hüllkurve nur dann zusammen, 

wenn die beiden Punkte ins Unendliche rücken. Dann fällt auch die wahre Grenzkurve mit der 
Hüllkurve in einem Geradenpaar zusammen. Dieses liegt parallel und symmetrisch zur o-Achse. 
Physikalisch bedeutet dieser Fall den Fließbeginn von bildsamen Metallen. Dann fallen auch 
die Kurven 7,(0,,) — die alle Scheitelpunkte der Spannungskreise verbinden — mit den erwähnten 
zwei Kurven zusammen. Hieraus ersieht man, daß das Fließen der bildsamen Metalle nur ein 
Sonderfall der allgemeinen Bruchtheorie ist. Es wird daher die Untersuchung der Grenzzustände 
spröder Stoffe von grundsätzlicher Bedeutung für die Aufstellung der allgemeinen Bruch- 
theorie sein. 
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Zur Berechnung von Gewölbestaumauern als Schalen 
Von H. Parkus in Wien 


Es sei vorausgesetzt, daß die Mittelfläche der Schale Kreisbogen als Horizontalschnitte 
besitzt, deren Mittelpunkte auf einer ebenen Kurve liegen. Diese wird meistens so gewählt, daß 
sich annähernd der gleiche Zentriwinkel für alle Bogen ergibt (,Gleichwinkelmauer‘‘). Be- 
zeichnet den von der Ebene der Mittelpunktskurve aus gezählten Bogenwinkel und z den 
nach abwärts positiv gezählten Abstand von einer durch die Mauerkrone gelegten Horizontal- 
ebene, so wählt man zweckmäßig das Koordinatennetz 9,=9, 92 auf der Schalenmittelfläche, 
welches aus den Kreisbogen z= konst und aus den diese schiefwinklig schneidenden Kurven 
9 = konst besteht. Aus der Parameterdarstellung der Mittelfläche: 


= r(2) cos P + f(2) | 
% — r(2)sin @ RE ne 2 a CE 
= 2 | 


wo z,—=f(z) die Gleichung der in der x,—2z;-Ebene liegenden Mittelpunktskurve bedeutet, 
lassen sich dann die geometrischen Kenngrößen der Mittelfläche, also der Maßtensor 9,5, der 


2) Hier kommt hauptsächlich das Gebiet zwischen der allseitig gleichen Zugfestigkeit a 4z , der Zug- 
festigkeit oz und Umschlingungsdrucktestigkeit opp in Frage, wo die Winkel deutlich schlecht übereinstimmen. 
- Weiter im Druckgebiet stimmen die Winkel für alle Belastungsarten gut überein. Einen Sonderfall wollen 

wir hier noch besprechen: Die Ebenen der Zugbelastungsart (o,—=0,, r—=0, enthält die o,-Achse) und der 
Druckbelastungsart (o,=0o,, r=1, enthält die o,-Achse) schließen mit der Projektionsebene den gleichen 
Winkel von 45° ein. Über den Grad der Verzerrung der Grenzkurve kann hier nur die relative Lage der 
Kurven zur Projektionsebene entscheiden. Die Kurve der Zugbelastungsart schneidet die o,-Achse (=07z) 
und läuft in weitem Bogen zur Spannung 'opp (schlechte Übereinstimmung der Winkel), während die Kurve 
der ee nstungesss ‚schneidet die A Rtchee (op) und läuft zur Spannung 077 (Umschlingungszug) mit 
einem der og-Achse nahe liegenden Bogen (gute Übereinstimmung der Winkel). (Abb. 2.) 
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Haupttensor b,;und der Maßtensor der sphärischen Abbildung c,;z sowie die Christoffelschen 
Dreiindizessymbole berechnen). Die Dehn- und Biegesteifigkeit der Schale folgen aus dem 
vorgegebenen Verlauf der Schalendicke. 

Eine analytische Lösung der Schalengleichungen dürfte auch bei radikaler Vereinfachung 
kaum möglich sein. Man wird daher zweckmäßig durch Einführung finiter Differenzen von den 
Differentialgleichungen auf Differenzengleichungen in den drei Verschiebungskomponenten u* 
(in der Mittelfläche) und w (senkrecht zur Mittelfläche) übergehen und diese Differenzenglei- 
chungen dann entweder direkt oder mittels eines Iterationsverfahrens (Relaxationsmethode) lösen. 

Die in jedem Netzpunkt anzuschreibenden drei Grundgleichungen lauten ?), ®) 


Open gen + SE tR=0 | 
ER N 3 
beim ot Ca Auer Be Pe 


Das Zeichen 9/dg° bedeutet die absolute (kovariante) Differentiation. n.s und &*P sind die 
Komponenten des ‚Krafttensors“ bzw. ‚Momententensors‘, die mit den Schnittkräften und 
Schnittmomenten gemäß 

DErP 


P,—= na» Ms — &Pay,, a ra Yayy au Bee 


zusammenhängen. »z ist der kovariante Normalenvektor der die Nee der Mittelfläche 
enthaltenden Schnittebene. 

Die Ausdrücke für den Krafttensor und den Momententensor als Funktionen der Kompo- 
nenten u°, w des Verschiebungsvektors der Mittelfläche gewinnt man am einfachsten mit Hilfe 
des Satzes vom Minimum der potentiellen Energie®?). Einen von den üblichen Annahmen der. 
Schalentheorie abweichenden Weg sind Green und Zerna gegangen®). 

Die Mauer besitzt an der Krone einen freien und längs der Talflanken sowie der Talsohle 
einen elastisch eingespannten Rand. Die elastische Einspannung berücksichtigt man nach 
Vogt?) zum Beispiel in der Weise, daß man sich die Schale um ein bestimmtes Stück in den 
Fels hinein verlängert denkt und dort starre Einspannung annimmt. Es gelten dann längs der 
Flanken und der Sohle die vier geometrischen Randbedingungen 

ow 


BR ee — a ee, ER he 
u uUa——W 0% Sp 0 (4), 


0 
wobei die letzte Bedingung auch durch 0 an den Flanken bzw. a an der Sohle 
ersetzt werden kann, da längs dieser Ränder w identisch verschwindet. An der Mauerkrone 
gelten die vier dynamischen Randbedingungen 
ar 0, aM?—0, a9 = ae A, 


Der geklammerte Index soll andeuten, daß es sich um die ri 950 handelt, so 
daß also für den Normalenvektor (5%, 


l....222.. 0). 


1 


pP: 


Er: em 0; (8% = . (6) 


gilt. Q ist die Kirchhoffsche Ersatzquerkraft. 


Spannungs- und Formänderungszustand in Hochdruckrohren 
mit ausmittiser Bohrung 


Von W. Wuest in Göttingen 


Im Schrifttum existieren voneinander abweichende Lösungen für den Spannungszustand 
in exzentrisch gebohrten Rohren unter Innendruck. Dieser wird daher erneut abgeleitet, wobei 
es sich als wesentlich erweist, die bei zweifach zusammenhängenden Flächen geforderten Umlauf- 
relationen zu erfüllen. Die Gültigkeit der Weinelschen Formeln für die Größtwerte der 


1) Griechische Indizes durchlaufen stets die Werte 1, 2. 

2) E.A. Deuker: Z. angew. Math. Mech. 23 (1943), 8. 81, 169. 

®) H. Parkus: Österr. Ing.-Archiv 4 (1950), 8. 160. 

*) A.E. Green und W. Zerna: Quart. J. Mech. appl. Math. 3 (1950), 8.9. 
5) F. Vogt: II. Weltkraftkonferenz. Berlin 1930. Bd. IX. 
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Tangentialspannung an der Innen- und Außenseite der schwächsten Stelle des exzentrischen 
Kreisquerschnittes wird auf diese Weise bestätigt. Die so gewonnene Lösung wird dazu benutzt, 
die Krümmungsänderung exzentrisch gebohrter Rohre unter Innendruck zu berechnen. Die 
Rechnung liefert das einfache Ergebnis, daß die Querkontraktion die Krümmungsänderung im 
Verhältnis 1—2/m herabsetzt!). Der Krümmungseffekt wird zur Messung sehr hoher Drucke 
(im Bereich von 1000 bis 10000 at) technisch angewandt. Die Theorie wird in Übereinstimmung 
mit Meßergebnissen befunden. (Die ausführliche Arbeit ist im Ingenieur-Archiv 19 (1951) 
H. 1, S. 12—21 erschienen.) 


Druckverteilung unter elastisch gelagerten Kreisplatten 
Von A. Jung in Stuttgart 


Wird der elastische Halbraum 2> 0 durch eine rotationssymmetrische Belastung bean- 
sprucht, so ist die Verschiebung der Oberfläche z=0 mit der Abkürzung 


& 
9,(A) = [ Pı(e) ET SSR (1) 
0 
gegeben durch?) 


w()= u] TREE N AR BETT (2a). 


Die Druckverteilung zwischen Platte a 
Halbraum ist 


p(r)= IRAQ Adı(Ar)dA :. (2b). 
Die Durchbiegung der Kreisplatte (Bild 1) ge- 
nügt der Differentialgleichung 
1 
AAw= N (Pe(r) — Pı(lr)) ; 
a 1d @). 


Aare Tr Bild 1 


Die Lösung von (3) läßt sich in der. Form angeben 


| | Pelo)-eI(Ao)de— (A) + a Tı(Aa) + 004, (Ra) 
N: b AN 
| Herne) tern), Jo(Ar) dA 


Die Konstanten c, und c,lassen sich aus den Randbedingungen bestimmen. Damit das Integral (4) 
existiert, muß die Funktion 


1.) = 4 | | Pa(e) : oe Jo(Ao) de — H(A)+ 09a Jo(Aa) +60 A J,(Aa) | 


N 
ö 6) 
+ a4? J,(Aa) +03 aA? J, (Aa) | 
beschränkt und quadratisch integrierbar sein. Daraus ergibt sich, daß der Grenzwert 
iu REN N Pas ee WO Se (6) 


. ken muß, Al dieser ER Et man gerade zwei Gleichungen zur Bestimmung 


von c, und 4. - 
Macht man nun die Kenaktne: die Platte liege reibungsfrei auf dem Halbraum auf, so muß 


w(r) = u,(r) für rsao, 
?ı(r) =0 BER & 


1) Inzwischen konnte nachgewiesen werden, daß unter der Voraussetzung reiner Druckbelastung diese 


Beet auch bei beliebi en RT ok ahetgn gilt. 


FR H. Jung: Ing.-Archiv I (1950), 8. 178. 
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erfüllt sein. Mit (4) und (2a) erhält man die Integralgleichung?) 


En , | 91(4) Jo(Ar) dA — za fr pP Jo(A 0) RAN; (ae) 


ö ob (7) 
+05aA2J (Aa) + 3a A J, 0) Ju(Arpdi; r<sQa 
mit der Nebenbedingung 
[n@ aan) dr =0 für 0 
Macht man für 90) den Ansatz 
1 ; 
(A): ah J (aA) ana N EP IR IN 
so ist (8), wie einfach nachzuweisen ist, erfüllt. Setzt man (9) in (7) ein, so wird 
v Gm] pl I M ‚ı (aA) Jo(Ar)di 
Aa 1 
Br. 1 ale (10 
Bi -1 | / Pelo)- 0 FAo)de—ma N "I 1a dt anarhla)g - - (10). 


0 |_0 


—+ aA? J, (Aa) -+ 63 aA° J (Aa) Jar) da 


Aus den Randbedingungen und der Gleichgewichtsbedingung erhält man drei Gleichungen für A,, 
c, und c,. Die Konstanten A, und c, werden nun dadurch bestimmt, daß man (10) auf den 


Kreisen r—° («= 1,2, ...) streng erfüllt. Man erhält gerade n Gleichungen zur Bestimmung 
2 : 
der Ar (n== 1,12... % Junllch, 
Liegt die Platte im belasteten Zustand auf der Fläche r<« auf, so muß auf dem Kreis 


1 


r— die Ableitung ee endlich bleiben. Dies führt auf die Bedingung ; 
Aula) ED a are U TEE 3 


wobei zu beachten ist, daß A, als Funktion der A, und a durch die Gleichgewichtsbedingung 
gegeben ist. 


Stoßwellen in isotropen elastischen Medien 9) 
Von Hans Richter in Haltingen/Baden | 


Die üblichen Rechnungen über elastische Wellenfronten in festen Körpern sind akustischer 
Natur mit Hilfe des Hookeschen Gesetzes. Das letztere genügt nicht mehr zur Beschreibung 
echter Stoßwellen mit Entropiesprung; hierfür ist eine Zustandsgleichung erforderlich, die für 
endliche Deformationen und endliche Entropieänderungen den Spannungszustand angibt. Die 
folgenden Ausführungen beschränken sich auf ebene Stoßwellen, die in entspanntes isotropes 
Material hineinlaufen oder Entspannung desselben bewirken. »Zugrunde gelegt wird das isotrope 
Elastizitätsgesetz?) in seiner allgemeinsten Gestalt?), wobei der Spannungstensor P als Polynom 


2. Grades des ISRRNEHBESDEN Deformationsdeviators Z erscheint mit den Ableitungen der 


” Es ist 9 =0, da am Aubenätd der Platte die Querkraft = 0 sein muß, : Br 
2) Eine ausführliche Wiedergabe ist beabsichtigt. ö ET 
®) Die bei starken Stoßwellen entstehenden plastischen Deformationen treten wegen. der eringen zur ah, 
Verfügung stehenden Zeit erst hinter der Stoßwelle auf und wirken auf dieselbe im Sinne einer Schwächung 
im Laufe des Fortschreitens ein. Auf die Stoßwellenberechnung selbst sind sie ohne Einfluß. ne \ 
*) H. Richter: Das isotrope u a at Z. angew. Math. Mech, 28 (1948), 8. 200-2095 
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spezifischen inneren Energie u = u (7, 9, 2, s) nach den Invarianten 7, y, 2 von L als Koeffizienten; 
s = Entropie. Die lineare Näherung liefert das Hookesche Gesetz mit vier Materialkonstanten, 
die sich z.B. durch den adiabatischen Volumkompressionsmodul, den Schermodul, den Wärme- 
ausdehnungskoeffizienten und eine spezifische Wärme ausdrücken lassen. Bei der quadratischen 
Näherung treten sieben weitere Konstanten auf. 


Bei stationärer Betrachtung der Stoßwelle wird die Stoßwellengeschwindigkeit », zur 
senkrecht zur Front stehenden Anströmgeschwindigkeit im spannungsfreien Ausgangszustand 7, 
die Abströmgeschwindigkeit besitzt die Komponenten v, senkrecht und %, parallel zur Front; 
d sei der Abknickwinkel der Stromlinien. Auf einer Kontrollfläche parallel zur Front im End- 
zustand ZI herrschen die Normalspannung o, und die Tangentialspannung r,. Die vier Erhaltungs- 
sätze für Masse, Impulskomponenten und Energie liefern drei Gleichungen zwischen ®,, ®, %,, 
O3, Tg, Wo. Bei Bildung von L in II lassen sich alle Verzerrungen und der Winkel @ zwischen 
Hauptachsen und », durch die „‚natürlichen Dehnungen“ £ und n ausdrücken. Die Verwendung 
des Elastizitätsgesetzes liefert dann drei Gleichungen für £&, 7, s, und den ‚inkompressiblen 


1 BR - 
Staudruck“ q= > v7. Bei Lösung durch Potenzreihenansatz ergeben sich wie im akustischen 
De 
Falle zwei Lösungstypen: 
Typ K. Reine Kompressionsfront. Maßgeblich ist neben der akustischen Longi- 
tudinalwellengeschwindigkeit c, eine Materialkonstante F, der quadratischen Näherung des 
0a 


Fon c 
-eı fest, so wird 


Elastizitätsgesetzes. Legt man die Stoßstärke durch die Verdichtung & = 
eı 


für Stoßwellengeschwindigkeit und Entropiesprung 


\ 1. vs Fk 5 
vla=14+fFk:w+:--, a 1=zä' 777 Be 
c 1 


mit ©, = absolute Temperatur in /. Bei Fx>0 ist die Front stabil, und v, wächst mit der 
Stoßstärke; bei Fz< 0 ist eine Verdichtungsfront unmöglich. Versuche von Pack-Evans- 
James!) ergaben für Stahl keine meßbare Erhöhung von v, gegenüber c,: Fx ist also, wenn 
überhaupt positiv, sehr klein. Bei Blei dagegen zeigte sich eine Erhöhung von ce; = 2150 m/see 
auf v, = 3020 m/sec bei einem Stoßwellendruck von 93 = 270 000 at. Hieraus ergibt sich 
© =35%, und = 900° C für die Endtemperatur, was in Anbetracht der Größe von 9, kein Wider- 
spruch zur Voraussetzung des festen Aggregalzustandes in /I zu sein braucht. 


Typ S. Front mit Scherung. Neben der akustischen Transversalwellengeschwindigkeit «, 


£ ; E 2 r OR 8 $ 
tritt noch eine Materialkonstante F, der 2. Näherung auf. —" — 1, und o, sind hier von 2. Ord- 


6 
nung in der Stoßwellenstärke, die durch ® oder r, charakterisiert wird. Die Richtung der Haupt- 
spannungslinien weicht in 1. Ordnung von 45° ab. Der Entropiesprung ist von 4. Ordnung und 
unabhängig vom verwendeten Material stets negativ. Diese Fronten existieren also nur dann, 
wenn sie Entlastung hervorrufen; doch ist die Auffächerung bei Belastungsschubironten so 
gering, daß genügend genau mit einer Front gerechnet werden darf unter Berücksichtigung des 
Auftretens von o,, ® und der Korrektur an @. 


Ein Referat des Vortrages: 
W. Weibull, Bofors: Die Festigkeit von festen Körpern als Wahrscheinlichkeitsfrage 
ist nicht eingegangen. 


Ein ausführlicher Bericht über den Vortrag von M. Ros, Zürich: Die Ermüdung 
der Metalle wird in einem späteren Hefte der Zeitschrift erscheinen. 


1) D.C. Pack, W.M. Evans und H. J. James: The propagation of shock waves in steel and lead. 
Proc. Phys. Soc. 60 (1948), p. 1—8. 
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G. Strömungslehre 


Inkompressible Potentialströmungen durch Serade, unendliche 
Schaufelgitter 


Von M. Strscheletzky in Friedrichshafen am Bodensee 


Die klassische Methode behandelt das Problem der Umströmung von geraden Gittern durch 
konforme Abbildung. Für praktische Zwecke ist es wesentlich bequemer und anschaulicher, die 
Umströmung von Schaufelgittern aus Singularitäten aufzubauen, die innerhalb der Schaufeln 
oder auf ihrer Oberfläche kontinuierlich verteilt sind. 

Bei praktischen Berechnungen treten zwei Probleme auf: 

1. Ermittlung der Schaufelform des Gitters, wenn die Strömung vor und hinter dem Gitter 
vorgegeben ist, und 

2. Ermittlung der Wirkung des vorgegebenen Gitters auf die Strömung, die weit vor dem 
Gitter bekannt ist. 

In der vorliegenden Mitteilung werden diese beiden Probleme behandelt und zwar unter 
der Voraussetzung, daß die Schaufeln unendlich dünn sind und die inkompressible Flüssigkeit 
reibungsfrei ist. 

Die Länge Z der Schaufelsehne, der Neigungswinkel ß und die Teilung i kennzeichnen die 
Bauart des Gitters (Bild 1). Der Verlauf des Winkels ß(s) längs der Schaufellänge I, bestimmt 


Bild 1 


eindeutig die Schaufelform. Weit vor dem Gitter herrscht die gleichförmige Strömung t,. Unter 
dem Einfluß des Gitters ändert die Geschwindigkeit iv, ihre Größe und Richtung bis auf i,, 
die weit hinter dem Gitter ein gleichförmiges Geschwindigkeitsfeld bildet. 

Die Umströmung der Schaufeln erhält man dadurch, daß man einer kontinuierlichen Ver- 
teilung der gebundenen Wirbel y(s) auf den Schaufelkonturen eine Parallelströmung überlagert, 


die auf der betreffenden Ebene (x, y) in Abwesenheit des Gitters herrsche (m). Die Zirkulation 


der Geschwindigkeit um eine jede Schaufel ist dann 


L.= a yls) » den in ee De (1). 


Die lokale Geschwindigkeit w(w,, w,), 2 n een des Gitters im Aufpunkt BE. Yo) 
der Schaufelkontur entsteht, ist 


en) +40 RE SEN (2), 


a 


Se [ Ri; 
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v() 
2 


wobei die Komponente d + von den gebundenen Wirbeln induziert wird. Die Geschwindig- 


keit an ist zu der gebundenen Wirbelschicht bzw. Schaufelkontur tangential gerichtet. 


Da die Schaufelkontur mit einer Stromlinie der Umströmung zusammenfällt, soll die 
Normalkomponente der lokalen Geschwindigkeit an der Schaufeloberfläche verschwinden. Daraus 
folgt die Gleichung: 


te = -@). 


die den Zusammenhang zwischen der Schaufelform und dem Geschwindigkeitsfeld darstellt. 
Bei der Ermittlung der induzierten Geschwindigkeit v(v,, v,) pflegt man, der Methode von 
A. Betz folgend!), zuerst das Geschwindigkeitsfeld zu betrachten, welches an der Stelle der 


„0“-Schaufel unter dem Einfluß der übrigen Schaufeln entsteht, und dann dieses Feld mit den 
von der ‚‚0“-Schaufel induzierten Geschwindigkeiten zu überlagern. Damit erhält man: 


| 1 (Hy) y(S): ds ' 
%= — | as, 5)‘ y(s)- ds ; ee Yard Unı 
en unit tt 


(25) 
Kerle) 
ep A ae m B)rleo)sde rn 9). 
Yy = t IRZE So) (8) dis I Y a tt 
f R (25) (25) 
wobei: 
Sn HZN:R 
i Ü %—Y 
als, H)=—° ae Bi 
Er 2 meinen: 27. (y — Y) + (m — ©)? 
1 Enz t % % 
b(s,s)=— —: 22 \ BE 2 
(#50) 2 FREIES ECHEREIC EEE (Y — Y)? + (2 — 2)? 


bedeuten. Die Funktionen «a(s, s,) und b(s, s,) sind im ganzen Integrationsbereich stetig. Daraus 
kann man folgern, daß eine numerische Auswertung der ersten Integrale (4) immer möglich 
ist, sobald die Wirbelverteilung und die Kontur /, bekannt sind. Die Funktionen] «a(s, s,) und 
b (s, s,) hängen nur von der relativen Lage der Wirbel und des Aufpunktes ab und können deshalb 
ein für allemal gerechnet und in Form von Berechnungskurven dargestellt werden ?). 

Für die Ermittlung der von der ‚0“-Schaufel induzierten Geschwindigkeit (v”) ist es 
zweckmäßig anzunehmen, daß die Wirbelbelegung der Formel entspricht 


N 
& R,:cosnp 


Aber N ERIR ONE ER al; 


wobei 9 = arc cos (> 2 bedeutet. Für die Schaufeln, die genügend kleine relative Wölbung 
s 


haben und genügend breite Gitter bilden, kann man die gebundenen Wirbel nicht über die 


N 
A 
r 


£ 
E 
er 
R 
2 


Schaufelkonturen, sondern längs der Streckenprofile (8, l) kontinuierlich verteilen. Damit 
erhält man: 


R.» 
dl sin ß De sin ß = 77 sın NY 
2 TREE oT DREIER ARE 
| u 29078 2 sin y 
1 e, a See (6), 
weh. ost SE 
ar Mi 5 2 sin y 


wobei y == arc cos e = ist. Bei größeren relativen Wölbungen und engeren Gittern kann 


_ man aus dem Integrationsbereich (l,) eine genügend kleine Umgebung des Aufpunktes aus- 


1) A. Betz: Diagramme zur Berechnung von Flügelreihen. Ing.-Arch. II. (1931), S. 359371. 

« 2) Die vollständigen Berechnungskurven a (s, s,) und b (s, s,) können vom Verfasser als Pausen 
bzw. Fotokopien in der Größe 45 x 66cm (17200 mm) bezogen werden. 

Ih 
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schließen, für welche eine streckenartige Wirbelbelegung ohne weiteres angenommen werden 
kann. Die entsprechenden Berechnungsformeln für v; und v, werden im Vortrag angeführt. 

Bei dem zweiten Problem ist es nötig, die Wirkung eines vorgegebenen Gitters auf die 
Strömung zu ermitteln. In diesem Falle ist sowohl die Schaufelform als auch die Bauart des 
Gitters bekannt. Außerdem wird die Geschwindigkeit weit vor dem Gitter vorgegeben. Es wird 
die Verteilung der gebundenen Wirbel, die Zirkulation und die Auftriebzahl bzw. die Charakte- 
ristik des Gitters c,„ = f(a) gesucht). 

Aus der Gl.(3) erhält man unter der Voraussetzung & die Integralgleichung: 


wen [ronant 1 [taste 
(7) 


1 ; (n—X)- y(s) ds — 
K ie +7: [8607085 = nn Ce) = | 


tg P(s) 


di, 
oder, nach Auswertung der Integrale für RN Wirbelbelegung (f, 2), die lineare 
Gleichung: 


I sin = 
tg ß(s) - = 5, 5 Bo” 1 Du: Kınt 5, “ ru | 


I x cos ß va 
Br: et: Au Kong, Er) —0 
t = 0 m n=1 


u 


n= 
wobei 


E17 


a\S, $ 
NE b)n = b : ® » COS No 2 dp 
ö En 


und 


bedeuten. 

Die gesuchte Verteilung der gebundenen Wirbel soll außer der Gl.(7) bzw. (8) noch der 
Bedingung der tangentialen Umströmung der scharfen Austrittskante entsprechen. Daraus 
ergibt sich die Gleichung: : 


VOR RL 
n=0 - 


Um die Werte von R, zu ermitteln, pflegt man die Gl.(8) für » Aufpunkte der Schaufelsehne 
zusammenzustellen. Diese n Gleichungen bilden zusammen mit der Gl.(9) ein System von 
n-+1 linearen Gleichungen mit n+1 "Unbekannten, welches leicht zu lösen ist. 

Im Vortrag werden einige Zahlenbeispiele angeführt. Die Berechnungen können N: 
diesem Verfahren mit praktisch beliebiger Genauigkeit durchgeführt werden. Das Verfahren 
gilt für beliebige Bauart des Gitters und für beliebige Form der Schaufeln. 


Berechnung von Strömungsfunktionen 
Von Guntram v. Gorup in Darmstadt 


| Die Kontur © eines Zylinders, dessen Außenraum mit idealer nern ruhender 

Flüssigkeit erfüllt ist, sei zeitlich veränderlich. Gesucht wird die Strömungsfunktion nr 2 der 
im Außenraum entstehenden ebenen Strömung. S 

Voraussetzungen: 1. Kenntnis der Abbildungsfunktionen 2(z,t) zwischen C und. as 
Einheitskreis (E.K.) der r- -Ebene; die Funktion soll die beiden Bar ausnahmslos konform 


ineinander abbilden und das Unendliche fest lassen. 2. z(7, t) und =2(%1) seien meromorph in N 


7 
Zur Berechnung von W(z,1) = w(z,t) wird die Randbedingung 


E),- (En 


dr ER ES APE 


)a= Po —P. 
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ER UNER ROOT RER 
( „= in eine Funktionalgleichung für w umgeformt. Dabei wird der folgende Kunstgriff 


benutzt, der es erlaubt, die konjugiert komplexen Werte einer beliebigen analytischen Funktion 
/(@,t) längs E.K. durch die Randwerte einer anderen analytischen Funktion f(r, t)* aus- 
zudrücken: 


IE9=f® en 1) EIN EN. a RESER URS Be} 


Man erhält so die Funktionalgleichung: 


1 RL, 12" 
Be u 0 ee I RR RO 1 RE (RE 


Mit Hilfe der (wegen der Meromorphie möglichen) Zerlegung: 


= Rötz) +SHr). 2... a a 5 


(r im Außenraum regulär, im Unendlichen wie > verschwindend)kann man die Lösung von(5) 
anschreiben: 


w— R(ötz) +, SE ERRUENN IE Sean (5). 


Die darin enthaltene Zirkulation läßt sich durch Addition eines Termes IE . einer Vor- 
cV 


gegebenen Bedingung anpassen. 
Beispiel: Pulsierende Ellipse. 
Pr ERTL N E* 1 IHm=a\ 
| Abbildungsfunktion: ‚2 =Ir+m Be \ BUN ) 
Folglich: z =iIt+tm = 
| Au je + mr 
= : 


2 im, } gibt - S(2*z)=Imr. 


; , ; el 
Somit nach (5): w’ = (ll — mm) = — (Im — Im) er 


3 / ; | 
w=(lU— mm)Int—+ : (Im — Im) ee 


Integralgleichungen erster Art in der Strömungslehre 
und ihre Auflösung 


Von K. Nickel in Stuttgart 


Bei stationären, inkompressiblen und reibungsfreien Potentialströmungen läßt sich die 
Strömung um beliebig geformte Konturen auf folgende Art erhalten: Man belegt die Randkontur 
mit Wirbeln und bestimmt deren Stärke y so, daß überall auf dem Rand die Normalkomponente 
der Geschwindigkeit verschwindet. Bei der rechnerischen Durchführung kommt man dabei 
häufig auf Fredholmsche Integralgleichungen erster Art für y. In dem Vortrag wurden zwei 
Probleme behandelt: 


1. Die zweidimensionale Strömung um ein dünnes, schwachgewölbtes Flügelprofil führt 


nach Birnbaum (diese Zeitschrift, Bd. 3 (1923), S. 290 ff.) auf eine Integralgleichung, deren 
Lösung seit langem bekannt ist. Betrachtet man n Tragflügel, so kommt man auf ein System 
von n Integralgleichungen erster Art. Dieses läßt sich streng lösen, falls alle Flügel hintereinander 
liegen (abgedruckt in Mathem. Ztschr. Bd. 54 (1951), S. 81—96; für n—2 war das System schon 
von H. Söhngen: Mathem. Ztschr. Bd. 47 (1941), S. 16—33 gelöst worden), es lassen sich 
weiter geschlossene Ausdrücke für den Gesamtauftrieb und das Gesamtmoment angeben. — 
Falls die n Tragflügel beliebig angeordnet sind, konnte das Gleichungssystem zwar nicht gelöst, 
_ aber wenigstens auf eine einzige Fredholmsche Integralgleichung zweiter Art zurückgeführt 
"werden. — Da in der Prandtlschen Theorie der tragenden Linie Integralgleichungen derselben 


ER NR REN 


Bi: a 
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Art vorkommen, kann die gefundene Auflösungsformel auch auf solche dreidimensionale Probleme 
angewandt werden. Als Beispiel wurde die Auftriebsverteilung geringsten induzierten Wider- 
standes zweier benachbarter Tragflügel behandelt und der zugehörige Widerstand bestimmt 
(ausdrückbar durch elliptische Integrale). E | 

2. Analog zu dem Ergebnis von Birnbaum im Falle eines einzigen Flügels ergibt sich im 
Falle der Staffelwinkel B=0° oder ß=90° für die Zirkulationsverteilung y(&) eines unendlichen 
Gitters aus dünnen, schwachgewölbten Flügeln in zweidimensionaler Strömung die Integral- 
gleichung: 


df@a)_ 


dz ı, 2ru, 


b 
1 y(y) eig ka y)dy. 


a 


J 


Dabei ist nach Bild 1 die Skelettlinie eines der Gitterprofile auf der Strecke a <x <b durch / (2) 


gegeben, weiter ist a ei gesetzt, die Anströmgeschwindigkeit im Unendlichen sei %,. 


Die Auflösung dieser Integralgleichung lautet): 


BEN RR 
N ae tere nt we | k(aTa) de. 


IT 


/ sin k(y—a) sin k(b— x) 
[77 


«1» 
— — _— is 
—i— 


ee halt au 


2 
V3 


/ © it: 
Cq,= J 2 008,- 1 (cos ENT ae N 


6 y= — 60828 + etg2 n(sin?2 0 —2 ©cos’n) 
de sin am &a, = Ysin (nr + ©) sin (a— x) (cos —2 ctg2 nsin ®) 
& 6 l 1 y p EA NE ie 
re ET en Kay a 
sin —,— R da sin(r+a)’t 4 
. A 


Bild 2 i r Bild 3 


Als Gesamtauftrieb eines Flügels aus dem Verband ergibt sich j 


. b b i E 
FRes: _ dom (een er 
4= un YW) day b—_a | ds Ysink(b—a) 23 NE 2 
cos k | | 
Ye a 2 a - c Be 
b x Y „rg 

Für = “ —= & (Streckenprofil) findet man hieraus z. B.: 4 re! 


2 8U, sinkb—y) | { a 
„W)=-—— an A= FR gr ik? ; 5% 
 cosk sin k(y— a) ; nah ER 
in Bild 2 und 3 sind noch zwei weitere Beispiele aufgezeichnet. i En, | sh 


1) Auf der Tagung erfuhr ich von Herrn Dr. J. Dörr- Darmstadt, daß ihm die s. 
bereits seit drei Jahren bekannt ist; eine Arbeit darüber erscheint im nächsten Heft des Ingeni 
Te REN 


% 
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Zur Theorie der unstetigen Strömungen 
Von Richard Eppler in Stuttgart 


Die wichtigste Methode der Hydromechanik, die es bisher ermöglichte, auf rein theore- 
tischem Weg Widerstandsbeiwerte zu berechnen, ist die Theorie der unstetigen Strömungen, die 
auf Helmholtz und Kirchhoff zurückgeht. Sie liefert aber nur dann einen (viel zu niedrigen) 
Widerstand, wenn die Totwasserbreite parabolisch unendlich wird, was mit der Wirklichkeit 
schlecht übereinstimmt. Es wurde schon vielfach versucht, diese Mängel zu beseitigen. Dabei 
wurden aber immer gewisse Daten aus dem Versuch entnommen, meist der Totwasserdruck, der 
ja bekanntlich beim Experiment immer tiefer liegt als der in der Theorie benützte statische 
Druck. Es wird hier über eine Methode berichtet, die rein theoretisch zu ausgezeichneter Über- 
einstimmung mit dem Experiment gelangt. 

Zunächst wurden mit Hilfe der bekannten Methoden der konformen Abbildung Totwasser- 
strömungen dargestellt, bei welchen längs der sich von der Kontur ablösenden freien Stromlinien 
die vorerst willkürliche Geschwindigkeit v, herrscht 
(Bild: 1). Sie haben endliche Länge und münden 7 mare" 
(bei Symmetrie zur Staustromlinie) in zwei zur , 
Anströmrichtung parallele Geraden ein. Diese 
geradlinigen Stromlinien können theoretisch nicht y 
existieren, da auf ihrer einen Seite der konstante "> > -- 
Totwasserdruck ?ı herrscht, der nach der Ber- Bild 1. Die berechnete Strömungsform 
noullischen Gleichung der Geschwindigkeit v, 
entspricht, auf ihrer andern dagegen die Geschwindigkeit von v, nach dem im Unendlichen herr- 


schenden Wert v, geht. Dies schadet aber nichts, da sich im Bereich dieser Geraden ohnehin 


das Totwasser in Wirbel auflöst. 

Es wird nun nur angenommen, daß die Totwasserauflösung an der Strömung bis zum 
Beginn der Geraden nichts ändert. Man kann sich theoretisch überlegen, wie weit diese Annahme 
gut und berechtigt ist. Im Endeffekt wird darüber aber der Vergleich mit dem Experiment 
entscheiden. Die Folgen dieser Annahme sind wichtig. In ihr stecken nämlich Anfangsbedin- 
gungen für die turbulente Bewegung bei der Totwasserauflösung, mit denen man diese Ver- 
mischung weiter behandeln kann, um eine Beziehung für die noch freie Größe v, zu erhalten. 
In dem einfachen Fall, in dem die Vermischung so erfolgt, wie es etwa Tollmien !) annimmt, 
in dem also die „‚Nachströmung‘“ effektiv auf der ganzen Totwasserbreite erfolgt, reichen die 
Anfangsbedingungen in Verbindung mit der immer zu erfüllenden Energieforderung, daß die 


von der Verwirbelung verzehrte Energie gleich der vom Widerstand geleisteten ist, schon aus, um 


ve Van elll) 
festzulegen. Man kann dabei noch zeigen, daß 
ısch dieser Wert höchstens etwas erhöht, der 
Druck im Totwasser also erniedrigt, wenn die 
Nachströmung nur im Mittel so erfolgt, wie wir 
es angenommen haben. 

Mit dem Wert von v, nach (1) können 
nun, genau wie mit v9, —=v,, rein theoretische 
Strömungen berechnet werden, und es zeigt 
sich, daß sie bei relativ breiten Konturen 
überraschend gut mit den Experimenten über- 
‚einstimmen. Bei der ebenen Platte senk- 
recht zur Anströmrichtung ergibt sich bei- 
spielsweise nicht nur mit c„=1,79 ein mit | | 
den Experimenten (0, Zu 1,80 bis 2,01) gut Bild 2. ee Or area ebenen Platte 
übereinstimmender Widerstandsbeiwert, son- 
dern auch die berechnete Strömungsform stimmt ausgezeichnet überein mit dem zum Vergleich 
herangezogenen Experiment von F. Walter?). Bei ihm wurde an der Vorderseite der Platte 
Farbe ausgesondert, die sich in die freien Strahlgrenzen fortpflanzt. Die berechneten freien 
'Stromlinien verlaufen ganz innerhalb der Farbspuren (Bild 2). Die punktierten Linien stellen 
zum Vergleich die freien Stromlinien der Helmholtz- Strömung dar. 


1) W.Tollmien: Berechnung turbulenter Ausbreitungsvorgänge. Z. angew. Math. Mech. 6 (1926), 
S. 468—478. 

2) F. Walter: Experimentelle und theoretische Untersuchungen über die Strömungsformen hinter 
scharfkantigen Widerstandskörpern, sowie Beziehungen zum Widerstandsproblem. Berichte der Mathe- 
matisch-physikalischen Klasse der Sächsischen Akademie der Wissenschaften zu Mira: Bd. XCII (1940), 
S. 139-234. 
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Auch die Kontur des Kreiszylinders wurde behandelt, wobei unter Heranziehung der 
Grenzschichttheorie das unterkritische Experiment in Widerstand und Druckverteilung sehr 
gut angenähert wurde. Beim überkritischen Experiment liegt der Totwasserdruck wesentlich 
höher als der (1) entsprechende Wert. Hier muß die Totwasserauflösung anders erfolgen, als 
zur Herleitung von (1) angenommen wurde. Man könnte aber vielleicht auch hier zu rein theore- 
tischen Ergebnissen kommen, wenn man die Strahlvermischung noch genauer untersucht. 


Ermittlung der Druckverteilung mittels Integralbedingungen 


Von K. Oswatitsch in Stockholm _ 


Zunächst sei nur die ebene, wirbelfreie, stationäre Strömung beliebiger Machzahl M be- 
handelt. Verdichtungsstöße sind in der Näherung wirbelfreier Strömung, wie sie in Schallnähe 
oder an schlanken Profilen gerechtfertigt ist, zugelassen. Die abhängigen Veränderlichen müssen 
dann längs geschlossenen Kurven gewisse Bedingungen erfüllen, wie etwa diejenige, daß das 
Integral der Tangentialkomponente der Geschwindigkeit oder der Normalkomponente der Strom- 
dichte verschwinden muß. Aus solchen Integralbedingungen läßt sich grundsätzlich eine Strö- 
mung beliebig genau ermitteln, wenn beliebig viele Bedingungen aufgestellt werden können. Dies 
ist bei ebener Strömung der Fall. 


Mit W, d als Geschwindigkeitsbetrag und -winkel lautet die gasdynamische a und 
die Gleichung der Wirbelfreiheit: 


| od r DAWN, 
mat ap 3 Wu Er (1). 
‘olgende Bedingungen lassen sich nun erfüllen: 
0 ’ ß) 
% an) ar) =0| 5 
Re | 2)» 


ID 


a (fecos v9)) + ee (g sin (v #)) = 0 | 


für beliebiges v, wenn die Geschwindigkeitsfunktionen f(W) und g(W) gewisse Differential- 


. bedingungen erfüllen, die sich nach Ausführen der Differentiation von (2) aus (1)ergeben. Durch 


Anwendung des Gaußschen Satzes auf (2) folgt schließlich mit ds und y als Bogenelement mit 
Riehtungswinkel der geschlossenen Kurve für #—= Wf und @= Wg: 


1, [ [Fin (#9) sin (d — y) — @ cos (vB) cos (#— y)] ds — 0: | (3) 
2. [| [F cos #9) sin (d — y) + @ sin (#9) cos (d — y)] ds = 0 J 3 
wenn @ und F folgende gewöhnliche Differentialgleichungen erfüllen: 
„d@G dF 
Wr (1 — M?) Won +1—-v?)1—M)G=0; Wil M)F+rO) (4). 
Explizit lassen sich folgende Lösungen angeben für: : 
EN a ARE EL für »—=1: — F=p+toW:; GENE Ewa 


was der Kontinuität, der Wirbelfreiheit und den Impulssätzen entspricht. Für M=0 lautet die 
Lösung von (4): 
F=6=Wrtt!, . für. beliebigesu? TEE Rn 


Die physikalische Bedeutung von (3) für v >1 ist noch unbekannt. Die Verallgemeinerung 
auf instationäre Strömung oder auf drei Dimensionen, wo Bedingungen auf Kurven und auf 
geschlossenen Flächen zu unterscheiden sind, ist nicht ohne weiteres möglich. Durch Abzählen 
der Variablen und der Bedingungen kann gezeigt werden, daß es im Raume für M=0 undv>1 
keine so einfache Bedingung wie Gl]. (6) gibt. 

Die Beispielrechnung steckt noch in den Anfängen. Mit einem einfachen 2-Parameter- 
ansatz erhält man folgende Maximalgeschwindigkeit an einem Kreiszylinder bei M „= 0,40: 


Wazlu, = 2,26 gegenüber einem Wert von 2,28 beim Janzen-Rayleigh- „Verfahren 4. Nähe- 
rung und 2,00 bei M,=0. 
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Berechnung von schallnaher Strömung an Profilen 
Von Tore R. Gullstrand in Stockholm 


Das Problem der Berechnung schallnaher Strömungen an symmetrischen, nicht angestellten 
Profilen ist mit Erfolg von Oswatitsch), ?2) behandelt worden. Eine Verallgemeinerung auf 


flache Flügel ist leicht und führt zu folgender Integralgleichung: (v- u Se): DU 
C 5 ee z 
U, Lösung mittels Prandil-Regel): 
7% 1 (8 Pte —2C- 2) 
. 2 2 SER rare, E— nz)? 
Gun -4+|lSc en ER ITEREEH 
ante) a Eee pr de ind: (1). 


Durch einen Ansatz für U? in y-Richtung kann eine Integration durchgeführt werden. Dieser 
wird abweichend von !) und 2) als eine Summe von y-Funktionen mit Funktionen von &, z 
als Koeffizienten gemacht: 


ul ehle)t 
mit: | 
AO)=1; fi); %0)=0;. f,0)=1; 


a hat die Dimension einer Länge und könnte noch von x, z abhängen. Aus (1) und (2) folgt: 


2). 


+ \ 
U 52 Ben L—g 
u) nl. „jasae 
“ (3). 
+ 
Bf RR € See G—2 & 5 g 
eh D de Br| RE a Un Blur ah: 


Bei ebener Strömung fällt die Integration über £ fort. Eine Integration ausgehend von U, 
und Einsetzen dieser Größe im Integral führt wie das Janzen-Rayleigh-Verfahren nur bei 
unterkritischer Strömung zum Ziel. Beilokalen Überschallgebieten erwies sich folgende Iteration 
als stark konvergent. Angenommen wird eine bestimmte, bei jedem Schritt gleichbleibende 
Stoßlage, nicht aber M.. Aus einer angenommenen Geschwindigkeitsverteilung, welche in die 
Integrale eingesetzt wird, wird M, und eine neue Geschwindigkeitsverteilung mittels (3) be- 
rechnet, mit der weiter iteriert werden kann. 

Als Beispiel diente vorläufig das Profil NACA 0006, wobei sich sehr gute Übereinstimmung 
mit den Versuchen und etwas kleinere Widerstände als bei ?) ergeben. 

Die Methode soll bis M»—=1 verwendet werden, wobei U,, nicht mehr mit M.,, sondern 
mit einer kleineren Machzahl zu linearisieren ist, damit es mit M&> 1 nicht über alle Grenzen 
wächst. 

Bei angestellten Profilen gilt die Integralgleichung nicht mehr in gleicher Form, weil die 
unsymmetrische Geschwindigkeitsverteilung die Randbedingung auf der .x-Achse beeinflußt, 
weshalb dort Zusatzquellen anzubringen sind. Hingegen heben sich die Einflüsse auf die Ge- 
schwindigkeitsverteilung von oben und unten gerade auf. 


Ausbreitung kugelförmiger Stoßwellen auf große Entfernung 
Von Fr. Wecken in Haltingen 


Zur Beschreibung kugelsymmetrischer Druckwellen, die in homogener Atmosphäre aus 
idealem Gas (mit p — 0?) auswärts verlaufen, wird eine zugleich strenge und allgemeine Lösung 
der isentropischen Grundgleichungen gesucht. Um der Störungsfortpflanzung längs Charakte- 


ristiken (C+ auswärts, CO einwärts) Rechnung zu tragen, werden charakteristische unabhängige + 


Variable eingeführt: &=const auf C+, y=const auf 07. r, t, w (Stoffgeschwindigkeit), 
a (Schallgeschwindigkeit) genügen als Funktionen von & und y vier Differentialgleichungen, die 
gelöst werden durch eine Potenzreihenentwicklung nach 1/s und A/s mit von @ abhängigen Koeffi- 
zienten; dabei ist s=y—z, A=In (s/s,), so eine Integrationskonstante. Die Entwicklungs- 


1) K. Oswatitsch: Z. angew. Math. Mech. XXX/l, 2 (1950), 8. 17—24. 
2) K.Oswatitsch: Acta physica austr. IV/2, 3 (1950), 8. 228—271. 
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"koeffizienten sind rekursiv bequem zu ermitteln. Eine darin auftretende willkürliche Funktion f(x) 
ermöglicht die Anpassung an vorgegebene Ausgangswerte. Die Anfangsglieder der Entwicklungen 
ergeben die bekannte akustische Lösung, wobei &—= a, 1—r, s=2r wird. Die Reihen konvergieren 
nur dort, wo die Wellenerregung hinreichend schwach ist, insbesondere für großes s; sie liefern 
zunächst nur eine stetige isentropische Welle. Die Abweichung von der akustischen Lösung 
kommt besonders in einem Glied const - A f’(z) zum Ausdruck, das die Koordinatenverzerrung 
bestimmt. 

Für große s (oder r) werden « und a als Funktionen von r und i bei jeder Wahl von f(x) 
mehrdeutig — ein Anzeichen für die tatsächliche Entstehung von Stoßwellen, wodurch die 
gewonnene bezüglich x und % eindeutige Lösung z.T. ungültig wird. Ist A= Fr r 


ist im „normalen“ Gebiet A>0, und Stoßwellen W,, W;, ... beginnen an denjenigen Punkten 
SS» Sg, . . ., wo auf der durch A=0 gegebenen Kurve E der (&, y)-Ebene t (und zugleich %) ein 
Minimum annimmt. Die Koordinaten der S; sind formelmäßig angebbar. 

Der weitere Verlauf jeder Stoßwelle W wird durch gewöhnliche Differentialgleichungen 
bestimmt. Werden Glieder dritter Ordnung in der Stoßstärke vernachlässigt, so ist auch hinter W 
die isentropische Lösung verwendbar. — Auf dem vorderen und hinteren Ufer (W®, W%) haben 
&, u, a verschiedene —, %, r, t gleiche Werte. — Die Stoßwellengleichungen lassen sich in erster 
überakustischer Näherung geschlossen integrieren und liefern die Beziehung 


Far) — f(@) = (ar — a , 


die sich an der f’-Kurve als ‚‚Flächenregel‘ deuten läßt. Da weiter 
a eye j\ 
Fe)—f(@) 2 | 

ist, läßt sich aus dem Verlauf von f’(z) leicht die Entstehung und Entwicklung der auftretenden 
Stoßwellen, insbesondere auch ihr Verhalten beir > oo, ablesen. Der Druckverlauf längs einer 0 
quer durch den Wellenzug ergibt sich aus der f’-Kurve durch affine Verzerrung, wobei gemäß 
der Flächenregel bestimmte Bögen durch Sehnen zu ersetzen sind. — Bei einer Explosion ent- 
stehen stets mindestens zwei Stoßwellen, bei einer Implosion mindestens eine. 

Nach hinreichend langer Laufzeit nehmen Explosionswellen stets die von Dumond (1946) 
und anderen vorausgesetzte Sägezahnform an, jedoch tritt dies in vielen Fällen erst außerhalb 
des technisch interessierenden Bereiches ein, so daß die Anwendung des hier angegebenen Ver- 
fahrens gerechtfertigt erscheint, das mit nur mäßigem Rechenaufwand eine wesentlich genauere 
Beschreibung liefert. 

An einem für die analytische Behandlung bequemen Beispiel (f(x) = sin x) lassen sich 
die Verhältnisse im einzelnen verdeutlichen. 


Auf eine Arbeit von Whitham (Proc. Roy. Soc. 1950), worin er dasselbe Problem in 


ähnlicher Weise behandelt, wurde Vf. erst nach dem Vortrag aufmerksam. 
Eine ausführlichere Veröffentlichung der Untersuchung ist beabsichtigt. 


Zur Algebra des gegabelten Verdichtungsstoßes 
Von W. Kofink in Karlsruhe 


Für den einfachen Fall des gegabelten Verdichtungsstoßes, bei dem über das ganze Strö- 
mungsfeld dasselbe x gilt und sowohl der Hauptstoß als auch der vordere Nebenstoß mit der- 


. selben dimensionslosen Anströmgeschwindigkeit w, angeströmt werden, läßt sich für eine passende 


Unbekannte weine Gleichung 6. Grades ableiten. Aus ihr erhält man alle am Gabelstoß interessie- 
renden Größen. Bekannte Variable sind die drei 


—1 oo —1 ,— 1 
RR: ee a Mi ee ee 2? 


_ die Unbekannte ist 
i On Se 1 ; 4 ‘ 
met! ma=M -1= Ho; > } 
ir 1— kon, 


die Variablen !, m, © laufen im Überschallbereich von 0 bis ©. Aus den beiden Gabel- Y 


bedingungen von A.Weise, nämlich Druck- und Richtungsgleichheit im Abstrom hinter dem 


Hauptstoß und dem hinteren A folgt zunächst eine Gleichung 12: Grades nt 2; aus. ©) 
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der sich zwei Paare trivialer Doppelwurzeln 2—=0 und x&—=m (vorderer bzw. hinterer Neben- 
stoß ist Machsche Linie) herausheben, so daß eine Gleichung 6. Grades übrig bleibt: 


6 
S A,(k, l, m) © —=0. Die Koeffizienten A, enthalten auch Z und m höchstens in 6. Potenz. 
N 


In Sonderfällen reduziert sich der Grad weiter: 1.) für  —=@ya,, 1= ©: Zerfall in eine Glei- 
chung 4. Grades und das Quadrat einer Gleichung 1. Grades; 2.) senkrechter Hauptstoß, I=m, 
Zerfall in das Quadrat einer Gleichung 3.Grades in x, 2. Grades in ]; 3.) <&=0, nichttriviale 
Lösung für vorderer Nebenstoß — Machsche Linie, Gleichung 4. Grades in m, 3. Grades in ]; 
4.) = m, nichttriviale Lösung für hinterer Nebenstoß — Machsche Linie, Gleichung 3. Grades 
in m, 2. Grades in ]; 5.) hinterer Nebenstoß ist senkrecht, Gleichung 4. Grades in m, Gleichung 
2. Grades in Z und lineare in x, nachdem } und m bestimmt sind. Die Lösungen der Sonderfälle 
sind die geeigneten Ausgangspunkte für die Lösung der allgemeinen Gleichung 6. Grades. Eine 
ausführliche Darstellung erscheint in den Annalen der Physik Bd. 9, Heft 2/4, 1951. (Aus 

dem Inst. f. Gasströmungen der T. H. Stuttgart.) 


Interferometrische Messung von Grenzschichten an ebenen 
Platten bei Überschallgeschwindiskeit 


Von R. Ladenburg in Princeton, N.J. 


Mittelst eines Zehnder-Mach- Interferometer wurde die Dichteverteilung von Grenz- 
schichten an ebenen Platten in einem intermittierenden Windkanal bei der Machzahl M—=2,3 
gemessen. Durch Änderung des Druckes im Reservoir konnte die Dichte und die Reynoldssche 
Zahl Re in der Meßkammer in weiteren Grenzen geändert werden. Da Druck und Temperatur 
im Reservoir bekannt sind und da der Druck innerhalb der Grenzschicht konstant ist, konnte 
die Temperaturverteilungin der Grenzschicht aus der Dichteverteilung des Gasgesetzes berechnet 
werden. Die Geschwindigkeitsverteilung wurde im Fall laminarer Strömung aus der Temperatur- 
verteilung und der gemessenen Oberflächentemperatur der Platte nach einer von L. Crocco!) 
herrührenden Beziehung für die Prandtlsche Zahl 0.72 berechnet und mit den Resultaten der 
Croccoschen Theorie!) verglichen. Mit wachsendem Abstand x von der Schneide der Platte 
weicht die gemessene Dichteverteilung und die berechnete Geschwindigkeitsverteilung bei Re 
zwischen 3 und 9. 10° mehr und mehr von der theoretischen Verteilung für laminare Strömung 
ab und nähert sich der Verteilung für turbulente Strömung, wie sie für Re> 2-10% erhalten 
wird. Die vorläufigen Resultate lassen sich durch transversale Ansteckung (,‚,‚transversal contami- 
nation“) nach Charters-Liepmann?) deuten, die von der turbulenten Grenzschicht an den 
Wänden der Meßkammer herrührt, in die die Ränder der ebenen Platte hineinragen. Die Unter 
suchungen sind noch nicht abgeschlossen und werden fortgesetzt. 


Die ausgebildete Kanalströmung in einem rotierenden System 
Von H. Ludwieg in Göttingen 


Es wurde die ausgebildete Strömung in einem Kanal mit quadratischen Querschnitt unter- 
sucht, der um eine Achse rotierte, die senkrecht zur Kanalachse steht. Dabei zeigte sich, daß 
man bei genügend hoher Rotationsgeschwindigkeit des Kanals das Strömungsfeld in zwei Gebiete 
aufteilen kann. Nämlich in eine Grenzschichtströmung an denjenigen Seitenwänden des Kanals, 
die in Richtung der Normalen zur Rotationsachse verlaufen, und in eine Kernströmung, die 
das übrige Gebiet des Kanals ausfüllt. Dabei sind in den Gleichungen für die Grenzschichtströ- 
mung nur die Coriolis-, Zähigkeits- und Druckkräfte wesentlich, während man bei der Kern- 
strömung für eine erste Näherung nur die in Kanalrichtung fallende Komponente der Corioliskraft 
und das Druckgefälle zu berücksichtigen braucht. Unter diesen Vereinfachungen und unter 
Berücksichtigung der Anschlußbedingungen zwischen den beiden Gebieten konnte für laminare 
Strömung das Strömungsfeld und der Reibungsbeiwert für den Kanal ausgerechnet werden. Es 
zeigte sich dabei, daß der Reibungswiderstand in einem Kanal in einem rotierenden System 
sich ganz erheblich von den Werten unterscheidet, die für die Kanalströmung im ruhenden System 
bekannt sind. Die erhaltene Formel für den Reibungsbeiwert konnte noch dadurch verbessert 
werden, daß auch das Reibungsglied in der Gleichung für die Kernstörung mitberücksichtigt 
wurde. Durch eine numerische Integration der Differentialgleichung für die Kernströmung 


1) L. Crocco jr.: Monografie Scientifiche di Aeronautica, No. 3 (1946). 
2) A.C.Charters: NACA TN 891 (1943). 
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wurde dann noch der Einfluß des Beschleunigungsgliedes auf den Reibungswert abgeschätzt, 
wobei sich ergab, daß dieser Einfluß im allgemeinen nur gering ist. 

Zur Prüfung der abgeleiteten Formeln für den Reibungsbeiwert wurden Druckabfall- 
messungen jn einem schraubenförmigen Kanal vorgenommen, der um die Schraubenachse 
rotierte. Dabei ergab sich in dem Re-Zahlbereich, in dem die Strömung noch laminar war, eine 
sehr gute Übereinstimmung mit den Rechnungen. Aber auch in dem Re-Zahlbereich, in dem 
die Strömung bei der Messung bereits turbulent war, stimmten die gemessenen Reibungsbeiwerte 
noch gut mit denen aus der Rechnung für laminare Strömung ermittelten überein. Der Grund 
für diese ungewöhnliche Tatsache liegt darin, daß zunächst nur der Kern turbulent wird, was 
keinen großen Einfluß auf den Druckabfall hat, da ja bei der Kernströmung die Wirkung der 
Zähigkeit und damit der Turbulenz nur gering ist. Ein wesentlicher Einfluß der Turbulenz ist 
erst dann zu erwarten, wenn auch die Grenzschicht an den Seitenwänden turbulent wird, was 
aber erst bei erheblich höheren Re-Zahlen eintritt als das Turbulentwerden der Kanalströmung 
- im nicht rotierenden System. Die Messungen geben ferner Aufschluß über den Reibungsbeiwert 
in denjenigen Gebieten, in denen die Voraussetzungen der Rechnung nicht mehr gelten. 


Zur rationellen Berechnung der turbulenten Reibungsschicht 
mit Druckgradient an rauhen Wänden 


Von N. Scholz in Braunschweig 


Für die Grenzschichtströmung längs einer Wand ohne Druckgradient (ebene Platte) sind 
die turbulenten Reibungsschichtgrößen sowohl bei glatter als auch bei rauher Wand durch die 
Arbeiten von L. Prandtl [1] und L. Prandtlu.H. Schlichting [2] für beliebige Reynolds- 
sche Zahlen und beliebige Rauhigkeitsgrößen als bekannt anzusehen. Und zwar kann hiernach 
das turbulente Geschwindigkeitsprofil in der Grenzschicht bei glatter und rauher Wand durch 
ein und dasselbe universelle Geschwindigkeitsgesetz beschrieben werden (logarithmisches Ge- 
schwindigkeitsverteilungsgesetz der Grenzschicht). 

Es läßt sich nun zeigen, daß die bisher zur Anwendung kommenden rationellen Verfahren 
zur Berechnung der turbulenten Reibungsschicht mit Druckgradient bei glatter Wand aus 
den bekannten Verhältnissen an der ebenen Platte bzw. aus der Rohrströmung abgeleitet werden 
' können (vgl. hierzu H. Schlichting [3]). Es liegt daher nahe, für die turbulente Grenzschicht 
an der rauhen Wand durch sinngemäße Übertragung der Ergebnisse der rauhen ebenen Platte 
ganz analoge Beziehungen aufzustellen, wie sie für die glatte Wand mit Erfolg angewandt 
worden sind. Dieses Vorgehen soll im folgenden in Form einer Gegenüberstellung der Formel- 
systeme für glatte und rauhe Wand kurz angedeutet werden. . 

Das universelle logarithmische Geschwindigkeitsverteilungsgesetz für die Grenzschicht an 
der ebenen Platte lautet: 


| rauh 


“ey _ 5,7510 *% +5,56 la); “= 5,7510 7 RE TBRENEEN 

Ux 
(u(y) = — Geschwindigkeit in der Grenzschicht; v; = Yr,/o ee, 
keit; » = kinematische Zähigkeit; k = Korngröße der ‚Göttinger Sandrauhigkeit‘ al Zur rech- 
nerischen Vereinfachung können diese Gesetze durch Potenzgesetze angenähert werden von der 


Form: 
uy) _ 0 Y Y uy) _ Fr 
RR o N AAN n ol BEP) RE ab): 
Schreibt man diese Gesetze für den Rand der Reibungsschicht, also für y= ö:u= U auf, so 
erhält man die bekannte Form des Potenzgesetzes für die Geschwindigkeitsverteilung in der 
turbulenten Grenzschicht (z. B. 1/7-Potenzgesetz): 


“in _ 2 BEN TN ME a RE (3). 


Der Exponent » ist hierbei für glatte Wand eine Funktion der Reynoldsschen Zahl Uö/v bzw. 
für rauhe Wand eine Funktion der relativen Rauhigkeit %k/6. Die Konstanten O und » können so 
bestimmt werden, daß die Impulsverlustdicke # und die Wandschubspannung 7, für das exakte 
Gesetz (1) und die Näherung (2) übereinstimmen. Da wir letztlich für die Grenzschicht auch nur 
die beiden. Kenngrößen ® und r, berechnen wollen, bedeutet die Verwendung des Potenz- 
gesetzes (2) gegenüber dem logarithmischen Gesetz (1) überhaupt keine Vernachlässigung. 
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Für die Wandschubspannung können aus Gl. (2) die Ansätze abgeleitet werden: 


T & T & 
sUE 1UB18 ef era fie bie (4a); Ta 1910 1 ee RE ae (4b). 


Diese Ansätze werden für die Grenzschicht mit Druckgradient übernommen, wobei die Hypo- 
these zugrunde liegt, daß für ein Grenzschichtprofil mit gleichem U #/v bzw. gleichem #/k auch 
die gleiche Wandschubspannung vorhanden ist wie an der ebenen Platte (diese Annahme ist 
für glatte Wand bei beschleunigter und mäßig verzögerter Strömung experimentell als gesichert 
anzusehen). Zur Berechnung der Impulsverlustdicke 9 kann nun die Impulsgleichung für die 
Grenzschicht mit Druckgradient herangezogen werden. Eine geschlossene Integration der 
Impulsgleichung gelingt durch Einführung der Ansätze (4) für die Wandschubspannung und 
eines konstanten Wertes für das Grenzschichtdickenverhältnis 4 = ö*/8 (vgl. |3]). Das Ergebnis 
erhält die Form: 3 


(Toy, ae, . 5a); RL. eorae NEN ORE 


0 


wobei die auftretenden Konstanten a, x und c Funktionen von Ud9/v bzw. 9/k sind, die sich aus 
den Ergebnissen für die ebene Platte nach [1] und [2] ermitteln lassen. 

’ Wir vermerken hier noch, daß sich durch ähnliche Analogiebetrachtungen zwischen glatter 
und rauher Wand auch die Näherungsverfahren zur Bestimmung eines turbulenten Ablösungs- 
punktes auf die rauhe Wand übertragen lassen. Ferner ist erwähnenswert, daß der Übergang 
von der ausgebildeten Rauhigkeitsströmung zur hydraulisch glatten Grenzschicht mit dem vor- 
liegenden ‚Verfahren zwangsläufig erhalten wird. Als Beispiel wurde bisher der ebene Diffusor 
bei hydraulisch glatter und verschieden rauher Wand berechnet. 


Literatur 
[1] L. Prandtl: Zur turbulenten Strömung in Rohren und längs Platten. Ergebnisse der Aerodynamischen 


Versuchsanstalt Göttingen, IV. Lieferung (1932). 
[2] L. Prandtlu.H. Schlichting: Das Widerstandsgesetz rauher Platten. Werft Reed. Hafen, XV. Jahrg,, 


(1934), S.1. 
[3] H. Schlichting: Grenzschichttheorie, Karlsruhe 1951. 


Thermodynamische Vorgänge in turbulenten Reibungsschichten 
| Von F. Schultz-Grunow in Aachen 


Es gibt zwei Beobachtungsresultate, die sich mit den üblichen Anschauungen über den 
Wärmeübergang nicht erklären lassen. Es sind dies einmal die eigenartige Wirkungsweise des 
Wirbelrohres, Luft in einen erhitzten und unterkühlten Teil zu trennent), 2), und dann die 

$: überraschend hohe Untertemperatur im Windschatten eines bei höherer Machscher Zahl quer 
y angeblasenen Kreiszylinders®). Man hat hier Wärmeübertragung von niederer zu höherer Tem- 
peratur, was dem 2. Hauptsatz zu widersprechen scheint. 

Es wurde gezeigt*), daß sich diese Effekte durch einen thermodynamischen Vorgang 
erklären lassen, der bei turbulenter Bewegung in einem durch Feldkräfte, hier durch Zentrifugal- 

 kräfte, geschichteten Gas auftritt. Er beruht darauf, daß die Turbulenz jede polytrope Schichtung 
in.eine adiabate verwandelt als Resultat einer vollständigen Durchmischung. Insbesondere wird 

im Bereich des Polytropenexponenten 1< x <k bei dieser Mischung Wärme von niedrigerem 

‚auf höheres Temperaturniveau transportiert, was bereits W. Schmidt?) an Hand von Messungen 

- inder Atmosphäre, die durch die Schwerkraft geschichtet ist, bemerkt hat. Nach Ertel®) arbeitet 
hierbei die Turbulenz als Kältemaschine, die Luftteilchen durchlaufen bei ihrer turbulenten Be- 
wegung einen Jouleschen Kreisprozeß. Die Wärmemengen werden also nicht bei den Tempe- 
'raturen der Herkunftsorte der Luftteilchen, wie etwa an einer geheizten parallel ängeströmten 
Platte, sondern bei einer durch adiabate Kompression im Kraftfeld erhöhten Temperrturdifferenz 
ausgetauscht, was auf eine Intensivierung des Wärmedurchgangs durch Schichtung schließen läßt. 


x 02) @. Rangue: J. Physique Radium (7) 4 (1933), 8. 112. 
0000.82) R.Hilsch: Z. Naturforschg. 1 (1946), 8. 208. 
000.8) E. Eckert und W. Weise: Forsch.-Arb. Ing.-Wes. 13 (1942), 8. 246. 


1 . *) F.Schultz-Grunow: Kältetechnik 2 (1950), 8. 273. r 
2 Fr Naar nadt: ‘Der Massenaustausch in freier Luft. Hamburg 1925, S. 18. 
*) H. Ertel: Meteorol. Z. 56 (1939), S. 109. 
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Die Schichtung durch Zentrifugalkräfte ist auffallend intensiv, im Wirbelrohr von 5 mm 
Durchmesser entspricht sie 6000 m und i im Windschatten entspricht sie 5000 m atmosphärischer 
Höhendifferenz. Allgemein bestehen an zwei 
Vergleichsorten im Zentrifugal- und Schwere- 
feld gleiche Temperaturdifferenzen, wenn die 
Zentrifugalbeschleunigung gleich der mit 
dem Maßstabsverhältnis vervielfachten Erd- 
beschleunigung ist. 

Die Wärmeleitung in geschichtetem Gas 
läßt sich nach W. Schmidt in die gewohnte 
Anschauung, daß Wärme stets von höherer zu 
niedrigerer Temperatur strömt, durch Ein- 
führen der potentiellen Temperatur bringen. 
Die Betrachtung des Kreisprozesses ermög- 
‚licht eine Verschärfung dieses Begriffes, wo- 
durch die thermodynamisch erzeugte Wärme- 
menge erfaßt wird. Der Wärmeübergang aus 
einer geschichteten Strömung an eine Wand 
ist aber mit allen Einflüssen der Reynolds- 
schen und Prandtlschen Zahl der gleiche wie 
ohne Schichtung, weil an der Wand wegen 
der geringen Strömungsgeschwindigkeit keine 
Schichtung besteht. Allerdings ist der paradoxe 
Fall möglich, daß in einem kreisförmigen Kanal 
Wärme von der kälteren Innen- zur wärmeren 

Bild 1. Vergleich der theoretisch und experimentell am Außenwand transportiert wird. Än Stelle 

Wirbelrohr erhaltenen Temperaturen ? ß } = 

einer Kanalwand kann auch eine freie Strö- 

mung treten, so daß durch einen wärmeren turbulenten Luftstrom eine kältere Wand gekühlt 
werden kann und ebenso das Umgekehrte möglich ist. 


Aus diesen Vorstellungen heraus läßt sich mit Aussagen über die Strömung eine Theorie 
des Wirbelrohres entwickeln, die auch quantitativ mit den bisher vorliegenden Temperatur- 
messungen!) übereinstimmt, was Bild 1 veranschaulicht. Die im Wirbelrohr erzeugte Kalt- und 
Warmlufttemperatur ist dort über dem Verhältnis der Kalt- zur Gesamtluftmenge aufgetragen 
so, wie sie sich theoretisch und experimentell ergibt. Die Diskrepanzen rühren von einer un- 
richtigen Extrapolation der experimentell ermittelten Kaltlufttemperatur auf y=0 her und auch 
daher, daß im Experiment nicht zwischen Gas- und Wandtemperatur unterschieden wurde). 


Zur Berechnung der Drucktendenz bei Wärmezufuhr innerhalb einer 
isothermen Atmosphäre von konstanter Grundgeschwindiskeit 


Von H. Stümke in Denzlingen (Baden) “ 


Eine reibungsfreie Atmosphäre mit der konstanten Temperatur 7 und der konstanten 
Grundgeschwindigkeit U parallel zur x-Achse einer &,y-Ebene, die mit der Winkelgeschwindig- 
keit A/2 um die 2-Achse rotiert und als Ersatz der Erdoberfläche dient, besitzt die ungestörte 
Druckverteilung ” 


ParBsNe Staa 
RT | 


wobei b=  -; R=c,—.c, ist. Eine Wärmezufuhr in dieser Atmosphäre macht sich als 
Störfunktion in der Zustandsgleichung®), *) bzw. in der Kontinuitätsgleichung®) bemerkbar. 


Zur Berechnung der daraufhin einsetzenden Drucktendenz (= lokale Druckänderung 2,) eignet 
sich die Methode der ‚mittleren‘ Zustandsänderungen®), bei der die wirklichen Beschleunigungen 


% R. Hilsch: Z. Naturforschg. 1 (1946) S. 208. 

2) Ausführliche Veröffentlichung demnächst in Forsch.-Arb. Ing. -Wes. 

®) H. Stümke: Zustandsänderungen durch großräumige Wärmezufuhr im Innern einer anfänglich 
ruhenden Atmosphäre. Erscheint in der Met. Rundschau. » 

*) H. Stümke: Theorie der Weiterentwicklung von stabilen atmosphärischen Druckfeldern. Met. ” 
Rundschau 3 (1950), S. 241—24 3 und 4 (1951), 8. 7—13. I 


5) H. Stümke: Rationelle Behandlung der Konvektion bei der angenäherten Berechnung det Wetter. 
entwicklung. Erscheint in der Met. Rundschau. 


er, angew. M. 
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relativ zur rotierenden Erde durch die Variation des „‚geostrophischen“ Windes (Gleichgewichts- 
windes) nach der Zeit ersetzt werden (quasigeostrophische Rechnung). Bei kleinen Störungen 
lauten dann die dynamischen Gleichungen 


APv=— Up, En; APu= — ar ; gG=— Me: 
Die Zustandsgleichung wird 
IR 

Pı Er ® P, -| WW PR; ER P (o: ; = up, a w P,) az S, 

und die Kontinuitätsgleichung 
0, + Pu, -+ (Pv),-+ (Pu), = S,. 

Große Buchstaben bezeichnen hier Größen des Grundzustandes, kleine die additiv hinzutretenden 
Störungen. S, und 5, sind Störfunktionen. Elimination von %, ®, w,o, und geeignete Trans- 
formationen führen zu der elliptischen Differentialgleichung 


1 
Dix + Dyy-+ Du er D®=O(X,Y,Z), 


wobei ® die transformierte Drucktendenz, Q eine Störfunktion bedeutet. Bei ebenem Erdboden 


‘und verschwindender Wärmezufuhr am Boden lautet die Randbedingung 


im[®B4A+0®,]=0 mit a Be Ye —1 u 
— ARE EEE 3 =; ı = YV a LE, 
Z0 i 2xy1-+ u? a Yxgb 


Die Forderung, daß die Drucktendenz im Unendlichen absolut bzw. relativ zum ungestörten 
Druck verschwindet, führt zu den beiden Regularitätsbedingungen 


——(Z+unN —— (Z+uP) 
nr einlpre —0. bzw. lim |B.e’' +" —(, 
VX+ +20 VxX®+Y°:+Z2 +0 
Die Lösung dieser Randwertaufgabe läßt sich mit Hilfe der Greenschen Funktion herstellen, 
die Vf. in seiner Dissertation angegeben hat®). Der Eindeutigkeitsbeweis für die Lösung ist 
nicht trivial. Man muß dabei wesentlich die Tatsache benutzen, daß x größer als 1 ist. 

Im Spezialfall des ebenen Problems (Ebene senkrecht zur Grundströmung) stimmt die 
behandelte Aufgabe formal mit einem Ansatz zur Berechnung kleiner Gleichgewichtsstörungen 
überein, der von L. Prandtl für eine anfänglich ruhende Atmosphäre gegeben wurde”), sich 
aber auch sinngemäß auf eine strömende Atmosphäre übertragen läßt. 

In der Energiebilanz tritt gegenüber dem Fall der anfangs ruhenden Atmosphäre?) als 
neues Glied ein meterologisch wichtiger Ausdruck für die Änderung der kinetischen Ener- 


‚gie auf. 


Zur Theorie der Strömungslager 
Von @. Heinrich in Wien 


Als Strömungslager soll ein Lager bezeichnet werden, bei dem das Schwimmen des Zapfens 
durch Einpressen von Flüssigkeiten oder Gasen in den Lagerspalt erzwungen wird. 

Am einfachsten liegen die Verhältnisse beim Strömungsspurlager für einseitige Abstützung 
(Bild 1). Die Flüssigkeit oder das Gas strömt durch Kapillaren, die in einem Kranz ängeordnet 
sind, dem Lagerspalt zu und fließt dann durch den Spalt radial nach außen. Für die Druck- 
verteilung im Tragfilm gilt bei Vernachlässigung der Trägheitskräfte bei inkompressibler Strö- 
mung: Ap=0. Beim Austritt aus dem Lager muß der Außendruck p, herrschen (erste Rand- 
bedingung); die durch die Kapillaren zuströmende Flüssigkeitsmenge ist gleich der Quellung 
der Filmströmung an der Einmündung (zweite Randbedingung). Dadurch ist bei gegebener 
Spaltdicke ö Druckverteilung und abgestützter Axjalschub S bestimmt. Als Sicherheit gegen 


" Berühren wird & = En eingeführt. (Sa, — Axialschub für ö=0.) Esergibt sich: © — F—P« 


Pi — Pa 
(p, = Speisedruck, 9, = Einmündungsdruck.) ö wählt man etwa 0,1 mm, der Düsenkranz- 


_ radius kann nach dem Minimum der Pumpenleistung bei gegebener Tragkraft und Sicherheit 


berechnet werden. Die Kapillarabmessungen folgen aus der gewünschten Sicherheit. Der zum 
L 9% H. Stümke: Rotationssymmetrische Gleichgewichtsstörungen in einer isothermen Atmosphäre. 


Z. Geoph. 16 (1940), S. 127—149. yet 8 
: = L. pr ande Beiträge zur Mechanik der Atmosphäre (Mömoire prösent6 3 l’Association Möt6orologique 
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Minimum gehörige optimale Pumpendruck ergibt sich zu: (p, — —= 1,757-6: 9%. (Pm = 
mittlerer Lagerdruck). Bei turbulentem Zustrom durch die Drosselröhrchen ändert sich nur die 
Zustrombedingung, minimale Pumpenleistung und optimaler Pumpendruck bleiben gleich. Das 
TEN IE OuG 


55 (n = dynamische EN a re 


Reibungsmoment ergibt sich zu: Mk = — 


radius, ec, „— Umfangsgeschwindigkeit. 
Beim doppelseitigen Spurlager (Bild 2) ergibt sich der Stützdruck aus der Differenz- 
wirkung beider Lager. Wir setzen die beiden Spaltdicken: 6, — 6, (1-8) und, ö, = d, (1—e). 
(2 ö,, = gesamtes Lagerspiel.) 1/e ist einMaß für die Sicherheit gegen Berühren, wenn man eZ1 


voraussetzt. Der Stützdruck Sist dann proportional zu e (erste Näherung). 9 = un 


ist eine 
€ 


Lagerkennzahl. Bei gegebenem p,„ und e existiert ein Minimum der Pompänldeinee sowohl 
nach der Stärke der Drosselung im Zustrom, als auch nach dem Düsenkranzdurchmesser; dazu 
gehört ein optimale Pumpendruck. Bei turbulentem Zustrom liegt das Leistungsminimum etwa 
um 1, tiefer als bei laminarem. 

Das zylindrische Traglager (Bild 3) besitzt ein oder zwei Düsenkränze. Die Exzentrizität 
des Zapfenmittels gegen das Schalenmittel @, soll schon bei ruhendem Zapfen klein sein gegenüber 


dem mittleren Lagerspalt ö,,. Dann ist n = = ein Maß für die Sicherheit gegen Berühren. Die 


a 0 0 
Druckverteilung bei ruhendem Zapfen ist durch die Gleichung: p - (ö?7p) = (0 mit den zwei 
erwähnten Randbedingungen zu ermitteln. Bei Vernachlässigung von Gliedern höherer Ordnung 
in a ergibt sich wieder: Ap=0; dies läßt sich für den vorliegenden zweidimensionalen Fall q 
leicht lösen. Es existiert wjeder ein Leistungsminimum nach der Stärke der Zuflußdrosselung 
und nach dem Verhältnis von Lagerlänge zu Zapfendurchmesser, das bei turbulentem Zustrom 
etwa um !/, kleiner ausfällt als bei lJaminarem; dazu gehören wieder optimale Speisedrücke. 


Din ne stellt eine Lagerkennzahl dar. Bei rotierendem Zapfen ist der geometrische Ort 


des Zapfenmittels für variable Drehzahl ein Kreis mit dem Durchmesser a). 
Bei kompressibler Tragflüssigkeit kann man wegen der Konstanz der Enthalpie mit dem 
Boyle-Mariotteschen Gesetz rechnen. Man erhält dann für die Druckverteilung im kompres- 


2 
sibeln Tragfilm die Differentialgleichung: p R ol5« y— 125 yr)|= =. («= Geschwindig- 


Für ruhenden, zylindrischen Tragzapfen erhält man daraus, bis auf Glieder höherer Ordnung 
in ag: “ 


} | Hr = 


Ri; Durch die Transformation pP geht sie in AP=0 über. Ebenso lassen sich die Randbedin- 
gungen durch P ausdrücken. Das Minimum der isothermen Kompressorleistung bei vorgegebener 


; 
- B 
BR keitsvektor der an den Film angrenzenden beweglichen Wand.) 


ve Kennzahl p, = —" liefert alle optimalen Größen in Abhängigkeit von p}. Bild 4 zeigt das 
e i a’ %0 

E -Minimum der Kompressorleistung in Abhängigkeit von p,,, für inkompressible (strichliert) und 
a; kompressible Tragflüssigkeit bei laminarem und turbulentem Zustrom in einem einheitlichen 


ir dimensionslosem Maßstab. $ 
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Ein Referat des ESTOrteaBeR - 


; 

3 , » 
Ze ‚Durch ein Versehen wurden die Abbildungen zu dem Referat G. Heinrich: ZUR 2 sr 
‚Theorie der Strömungslager“ im Heft 8/9 dieses Jahrganges, Seite 295/296 nicht Ä 
mit abgedruckt. Wir bitten, dieses Blatt mit den zu der genannten Arbeit gehörenden R- 
‚Figuren hinter Seite 296 des vorigen Heftes einzufügen. Er 
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KLEINE MITTEILUNGEN 


"Eine einfache Vorrichtung zur Addition von 
Produkten. 


Schlüsselgleiehung: 2 = BI 


Das Folgende stellt einen Versuch dar, die bekannte 
Zeichenmethode für diese Schlüsselgleichung weit- 
gehend zu ‚‚mechanisieren“. Statt der Zeichnung 
werden an Kurvenscharen einstellbare Lineale. ver- 
wendet, die ein Spaltstück vom Betrage des Summan- 
den begrenzen, um welches dann eine Zunge verschoben 


4440444424422. N 
a Fe Fe Be TS 


Ya Y 
_——_— . —_—————_ g- 
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Di 

$ 
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werden kann, (Keine Zwischenablesung, nur End- 


resultat.) 
Die Vorrichtung zum Multiplizieren besteht aus den 
vier Linealen a, b, c, d. Lineale b und dinC bzw. E 
‚, rechtwinklig zu a, auf @ verschiebbar. Lineal ce um O 
auf a drehbar. Mt OO =y, AC=9,0E=x wird 


>  EF=o =". Die Worte von pi (zu, 9),% (Zu Un) 


%% (23, Y3), WO &1,Yı 29, Ya %z, Y, unabhängige Vari- 
able sind, könnten etwa aus Tabellen entnommen und 
auf Linealskalen eingestellt-werden, am besten aber 
benutzt man direkt Kurvenscharen für die @-, x-, Y- 
Funktionen zur Einstellung. (Anlegen der betr. 
Lineale an die Punkte x, y der Kurven, hierzu auch 
verschiebbare Einstellmarken Z und M auf b bzw. d.) 
Lineal « ist längs der Geraden g (Führung) verschieb- 
bar, um die Einstell-Lineale an die Kurvenpunkte zu 
er. bringen, die dafür oberhalb und unterhalb von g 
Br gezeichnet sind. Nach Möglichkeit wird man auch 
Br Kurvenscharen auf den Linealen selbst anbringen. 
“ (Gelegentlich kann auch < (c,a) als Unabhängige 
zweckmäßig sein.) 


Ess 3 "Basis an g, auf der die Kurven tragenden Deckplatte 
’; - angebracht, unter jedem Spalt eine Zunge Z, die durch 
 Einstechen eines spitzen Metallgriffel in Spaltrichtung 
verschoben werden kann. 

Nach Einstellung der Lineale wird E an die Basis 

des Spaltes gebracht, in dem addiert werden soll. Bei 
_ eingestochen und diese durch Herabziehen des Griffels 
’ eh Bf 2 i 


N AT nee % 


Senkrecht zu g ist eine Anzahl enger Spalten $, 


positiven Produkt wird dann am c-Lineal in die Zunge: 


bis Eum EF=o verschoben, bei negativem Produkt 
ist bei Z einzustechen und durch Hinaufführen des 
Griffels die Zunge in umgekehrter Richtung zu ver- 
schieben. Der Betrag der Verschiebung ist unten an 
der Zunge in einem Fenster @ abzulesen. Bei Addition 
in mehreren Spalten lassen sich beliebige Teilsummen 
bilden. Kann man bestimmten Spalten bestimmte 
%-Kurven zuordnen, so wird alles wesentlich einfacher 
(z. B. harmonische Analyse). 

Die erreichbare Genauigkeit hängt vom Maßstab 
und vom Material ab. Für rasche Überschlagsrech- 


nungen wird schon ein roh aus Pappstücken gefertigtes 
Gerät gute Dienste leisten. ; 

Zum Schluß darf noch der Sonderfall angeführt 
werden, daß Werte von Funktionen nur einer Ver- 
änderlichen addiert werden sollen: 


o=p(H)+xX(&)+v(2) +&la,)‘*- 


%1, %g, %3, %4, . . ., 2m unäbhängige Veränderliche. Hier 
wird man die Funktionsleitern direkt am Spalt an- 
bringen, für jeden Spalt ein oder zwei Leitern (beider- 
seits oder nur auf einer Seite des Spaltes). Die Zunge 
ist in diesem Fall für alle Spalten gemeinsam. 


Schwäb. Hall/Hessental. H. Schminke. 


Uber druckpunktfeste Profile. 

In der Strömungslehre sind Profile von Interesse, 
deren Auftriebsmittelpunkt bei Anstellwinkelände- 
rungen nicht wandert (druckpunktfeste Profile). Es 
sind auch viele Profile dieser Art bekannt (vgl. etwa 
Fuchs-Hopf, Aerodynamik, Bd. II, Berlin 1934, oder 
Durand, Aerodynamic Theory, Bd. II, Berlin 1935). 
Ein Beispiel für ein solches Profil ist jedes symme- 
trische Profil (Abb. 1a), es gibt aber auch unsymme- 
trische Profile mit dieser Eigenschaft (Abb. 1b). Der 
Vergleich mit einem nichtdruckpunktfesten Profil 
(Abb. le) zeigt, daß das Hinterende des Profils bei 1b) 
nach oben gebogen ist, wodurch es eine doppelte 
Krümmung erhält; man spricht deshalb von ‚S-Pro- 
filen‘‘ oder ‚Profilen mit S-Schlag‘“. 


Beim Betrachten des Bildes 1 erhebt sich die Frage,. 


ob es außer den symmetrischen Flügelschnitten a) und 


Kin 
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den S-Profilen b) auch Profile der Form ce) gibt, die 
keine Druckpunktwanderung zeigen; ob also die 
S-Form des Beispiels b) in gewissem Sinne!) zur Er- 
zeugung der Druckpunktfestigkeit notwendig ist. 


ee es 


Bild 1 


Diese Frage soll im folgenden untersucht werden, da- 
bei wird der Begriff der „S-Form‘‘ noch genauer 
präzisiert werden müssen?), 

Als Methode für diese Untersuchung wird die von 
Birnbaum?) eingeführte Theorie der dünnen, 
schwachgewölbten Profile in zweidimensionaler, in- 
kompressibler, stationärer und reibungsfreier Strö- 
mung verwendet. Man ersetzt dabei das Profil durch 


Bild 2 


seine Mittellinie (Skelettlinie). Wie in Abb. 2 skizziert 
ist?), möge sich die Profilmittellinie von #—=0 bis 
x—=1 erstrecken, sie werde analytisch durch y=f(z) 
gegeben. Da es in der Birnbaumschen Theorie nicht 
auf eine Drehung und eine Parallelverschiebung an- 
kommt, kann man f(0) = f(1) =0 wählen. Mit diesen 
Bezeichnungen ist nach Birnbaum?) das durch f(x) 
gegebene Profilskelett genau dann druckpunktfest, 
wenn das Den) um den Punkt x = 0,25: 


une |z, 


verschwindet. 

a) Bezeichnet man als ‚„S-Form des Profilskeletts‘“ 
die Eigenschaft, daß f(x) im Definitionsintervall ver- 
schiedenes Vorzeichen besitzt, wodurch also f’(0) 
und f’(1) gleiches Vorzeichen haben können (Bild 3a), 
so gibt es offenbar Profilmittellinien ohne S-Form, 
die nicht gerade sind und trotzdem keine Druckpunkt. 


\f(a)de . . (1) 


1) Wenn man nämlich von dem symmetrischen Profila) absieht. 

2) Inseiner (unveröffentlichten) Dissertation hat Herr R. Horten 
u. a. auch diese Frage behandelt. Das hier abgeleitete Ergebnis 
wurde unabhängig davon gewonnen und ist dasselbe, wie es schon 
Herr Horten erhielt. Da dort jedoch wesentlich höhere "Anforde- 
rungen an das Profil gestellt werden (u.a, sind Profile mit ge- 
kniekter Mittellinie ausgeschlossen) und da die dortigen Ergebnisse 
mit einem erheblich größeren Rechenaufwand erhalten werden, 
un die nachfolgenden Überlegungen trotzdem mitgeteilt 
werden. 

») W. Birnbaum, Die tragende Wirbelfläche als Hilfsmittel 
zur Behandlung des ebenen Problems der Tragflügeltheorie. Diese 
Zeitschrift, Bd.3 (1923), 8.2901. 

*) Zur besseren Veranschaulichung ist in Bild 2 noch gestrichelt 
eine Profilkontur eingezeichnet, 
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wanderung besitzen. Ein einfaches derartiges Beispiel 
ist (Bild 3b) 


F(x)=2(1— .) (39 — 1112 + 109 2? — 36 »°)?). 
Diese Mittellinie ist druckpunktfest, wie man sich 


. durch Einsetzen in (1) überzeugt, sie hat (vgl. Bild 3b) 
im Innern der Strecke (0, 1) positive Ordinaten, ferner‘ 


ist #°(0),=39>0 und F(I)=—1<0! 


Bild 3 


F(x) besitzt jedoch auf der Strecke (0, 1) einen 
Wendepunkt (bei 2 =5/12). Nachdem nun nach- 
gewiesen ist, daß es auch druckpunktfeste Profile gibt, 
deren Hinterende nicht „aufgebogen‘“ ist, erhebt sich 
sofort die weitere Frage, ob etwa der bei F(x) vor- 
handene Wendepunkt für die Druckpunktfestigkeit 
dieses Profils verantwortlich ist. Stellt man ne 
Frage etwas allgemeiner, so lautet sie: 


bh) Kann eine Mittellinie druckpunktfest sein, die 
auf der Strecke (0, I) keinen Wendepunkt besitzt? 


Die für die analytische Behandlung des Problems 
zu treffenden mathematischen Voraussetzungen sind 
gering und in den Anwendungen immer erfüllt, da 
dort außer Knicken in der Profilmittellinie keine 
weiteren Singularitäten an dieser vorkommen. Aus 
Symmetriegründen genügt es offenbar, sich auf Profil- 
skelette zu beschränken, die nach oben gewölbt sind. 
Man trifft folgende 


Vo pe EL f(xz) seiim Intervall < 0,1 > 
mit.f(0) = f(1) = 0 definiert und dort außer in den 
endlich, vielen Punkten &(n =]1, N) zweimal 
differenzierbar. In den Punkten a, gelte 


fin -)-f(an+)=An>20 (n=1,..,N). 


Weiter sei f’’(2)>0 und über <0,1> integrierbar. 


Behauptung: Es ist M(f)>0, dabei tritt das 
Gleichheitszeichen nur für f(x2)=0 ein. 

Beweis: Man formt (1) durch teilweise Integration 
um und erhält 


un=2[eV/sü af (a)da 
Er ER \ RT) 
+2 3 @n An Van (1— @n). | j 


Nach Voraussetzung sind beide Summanden der 
rechten Seite von (2) positiv und verschwinden nur 
für f(2J)=0 und „=0(n=]1,... N), was mit 
f(x) =0 gleichbedeutend ist. Damit ist die Behaup- 
tung bewiesen. 

Es zeigt sich also, daß, außer der geraden Strecke, 


eine Profilmittellinie ohne Wendepunkt nicht druck. 
punktfest sein kann. Es erscheint daher sinnvoll, 
schon dann von einer „Mittellinie mit S- Form“ zu 


sprechen, sobald sie (mindestens) einen Wendepuncii 
besitzt. Mit dieser Sprechweise gilt dann: 

Das einzige druckpunktfeste Profilskelett oäne S-Form 
ist die gerade Strecke. 


Stuttgart. 
®) Dieses Beispiel wurde 1 o. E Dissertation 


Horten entnommen, in der eine große Anzahl 29 
mittellinien angegeben und näher untersucht wird 
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BUCHBESPRECHUNGEN 


Dr. Walter Weizel (o. Prof. der Physik a. d. Univ. 
Bonn), Lehrbuch der Theoretischen Physik. 
Zweiter Band: Struktur der Materie, XII -+ 767 8. 
mit 194 Abb. Berlin/Göttingen/Heidelberg 1950. 
Springer-Verlag. Preis: brosch. 66,— DM. 

Inhalt: Elementare Atomtheorie, Quantentheorie, 
Moleküle, Chemische Bindung, Statistik, Struktur und 
Eigenschaften der. Gase, Struktur und Eigenschaften 
der zusammenhängenden Materie, Kernphysik, Ele- 
mentarteilchen. 

Diese Zusammenstellung der. Kapitelüberschriften 
gibt zwar den äußeren Rahmen, der den Inhalt des 
Buches umspannt, gewährt aber keinesfalls einen voll- 
ständigen Eindruck von der Reichhaltigkeit des In- 
halts. Man ist erstaunt über die Manigfaltigkeit der 
aufgeworfenen Probleme, wenn man das Buch im 
einzelnen durcharbeitet. Tatsächlich wird die gesamte 
moderne Atomtheorie umfaßt, und die Hauptprobleme 
sind lehrbuchmäßig in verständlicher Form abge- 
handelt. Dadurch nimmt das Buch unter den Lehr- 
büchern der theoretischen Physik eine Sonderstellung 
ein, denn es gibt in der deutschen Literatur kein 
Werk, das bei ähnlich geringem Umfang gleich reich 
an Inhalt ist. 

In den wesentlichsten Teilen ist die Darstellung so 
ausführlich, daß fortgeschrittene Studenten leicht 
folgen können. Schwierigere Kapitel bzw. solche, in 
denen die Ableitungen nur angedeutet sind, tragen 
eine besondere Kennzeichnung. Wie im ersten Band 
stehen kurze Inhaltsangaben und eine Zusammen- 
stellung der neu eingeführten Größen am Kopf der 
einzelnen Abschnitte. Auch viele Einzelheiten ver- 
dienen wegen ihrer Vorzüge besonders erwähnt zu 
werden. So wird das praktische Maßsystem mit vier 
Grundeinheiten in erfreulicher Konsequenz auch in 
der Atomtheorie angewandt; man findet z. B. nichts 
über eine angebliche Verletzung des Kausalitäts- 
prinzips in der Quantenmechanik. 

Nach Ansicht des Referenten hat das nun voll- 
ständig vorliegende Werk alle Eigenschaften, die es 
zu einem bevorzugten Standardlehrbuch für die 
Studenten der Theoretischen Physik machen können, 


Dresden. Recknagel. 


Adalbert Duschek, Vorlesungen über 
höhere Mathematik. Il.Bd. VI-+ 386 S. 
mit 125 Abb. Wien 1950. Springer-Verlag. Preis 
26,— DM, geb. 29,— DM. 

Erfreulich schnell ist dem ersten Band, der in dieser 
Zeitschrift Bd. 31, S. 58 besprochen wurde, der zweite 
gefolgt. Er bringt in seinen sechs Abschnitten zu- 
nächst einen Abriß der Theorie der unendlichen 
Reihen, insbesondere Potenz- und Fourierreihen. Es 
folgen die Differentiation und Integration bei Funk- 
tionen mehrerer Veränderlicher. Hier haben auch die 
wichtigsten Sätze über Kegelschnitte, Koordinaten- 
transformationen, Abbildung von Ebenen und Räu- 
men, einiges aus der Theorie der ebenen Kurven und 


. die Grundbegriffe der Vektorrechnung ihren Platz 


gefunden. Besonders zu begrüßen ist die nun folgende 
Fortsetzung der Wahrscheinlichkeitsrechnung, die bis 
zu den Theoremen von Bernoulli, Poisson 
und Bayes führt. Man kann dem Verfasser nur 
zustimmen, daß dieses Gebiet heute in die Grund- 
vorlesung gehört. Es schließen sich ein kurzer Abriß 


. der Ausg ESTER ER UUDE und der Annäherung em- 


pirischer Funktionen an, Der Rest des Buches ist der 
linearen Algebra (Determinanten, Matrizen, Tensoren) 
und der Tensoranalysis gewidmet. Dabei bevorzugt 
der Verfasser eine Darstellung, wie er sie in den von 
ihm gemeinsam mit A. Hochrainer verfaßten 
„Grundzügen der Tensorrechnung in analytischer Dar- 


‚stellung‘ gegeben hat. Hier werden aber die „sym- 
7 olischer oh reikieike (Bezeichnung der Vektoren 


mit deutschen Buchstaben usw.) und die „symbol- 
freie“ (mit Indizes) nebeneinander verwendet, damit 
sich der Leser selbst ein Urteil über deren Vorzüge 
bilden möge. Den Aufgaben des Textes sind wiederum 
die Lösungen beigegeben. 

Auch dieser Band besitzt die Vorzüge des ersten: 
die Darstellung ist im besten Sinne modern. Die Stu- 
dierenden der angewandten Mathematik, Physik und 
der Ingenieurwissenschaften sowie der in der Praxis 
stehende Forscher finden in diesem Buche erstmalig 
manches, was sonst nicht in Lehrbüchern für diesen 
Benutzerkreis aufgenommen zu werden pflegte. Ihnen 
allen kann das Werk bestens empfohlen werden. 


Aachen. Günther Schulz. 


Ernst Suter, Die Methode der Fest- 
punkte, Vereinfachtes Verfahren zur Berechnung 
statisch unbestimmter Konstruktionen mit Beispielen 
aus der Praxis, insbesondere von Stahlbetontrag- 
werken. Dritte, neu bearbeitete Aufl. von 
Dipl.-Ing. Ernst Traub. XIX-+ 2168. mit 
232 Abb. und 7 Tafeln. Berlin 1951. Springer-Verlag. 
Preis geb. 21,— DM. 

Die Statik neigt in ihren Berechnungsmethoden 
immer mehr nach der analytischen Seite. ;Dadurch 
geht oft die Anschaulichkeit verloren und Fehler- 
quellen können sich leichter ausbreiten. Um so mehr 
ist es zu begrüßen, wenn bewährte graphische Metho- 
den für die gegenwärtigen Bedürfnisse der Praxis 
weiterentwickelt werden. Das ist in der dritten 
Auflage des Buches von Su ter über die Methode der 
Festpunkte geschehen, die von einem erfahrenen 
Praktiker des Stahlbetonbaues, Dipl.-Ing. Traub, 
ganz neu bearbeitet wurde. 

Das bekannte Seilpolygonverfahren — von Cul- 
mann, Mohr und Ritter in seiner klassischen 
Fassung herausgebracht — führt bei der Berechnung 
durchlaufender Träger zu den sogenannten Fest- 
punkten. Es ist seinerzeit von Suter auf die 
Rahmentragwerke erweitert worden, konnte sich 
aber in der Praxis nie ganz durchsetzen, da die Er- 
mittlung der Festpunkte, der Verteilungsgrößen und 
der Übergangszablen etwas zeitraubend ist. Durch 
eine Arbeit von Guldan über Rahmentragwerke 
und Durchlaufträger sind wesentliche Vereinfachungen 
bekannt geworden, die Tr a u b mit Vorteilauswertet. 

Gegenüber der zweiten Auflage ist nicht vom all- 
gemeinen Fall der Tragwerke mit veränderlichem 
Trägheitsmoment ausgegangen, sondern von den in 
der Praxis am meisten vorkommenden Tragwerken 
mit von Stab zu Stab veränderlichem, aber auf Stab- 
länge konstantem Trägheitsmoment. Dadurch ist 
das Einarbeiten in diese Methode sehr erleichtert. 
Im übrigen ist das Verfahren einfach und übersicht- 
lich, die Berechnung kann teils analytisch, teils 
graphisch durchgeführt werden, so daß in jedem 
Stadium der Berechnung der Verlauf der Kräfte ver- 
folgt und d&mzufolge auch leicht nachgeprüft werden 
kann. Gut gewählte Beispiele aus der Praxis des 
Stahlbetons, Hilfstafeln und Richtlinien für die 
praktische Aufstellung der Berechnung ergänzen das 
Buch. 

Dresden. Henn. 

F. Schwank, Randwertprobleme und 
andere Anwendungsgebiete der höhe- 
ren Analysis für Physiker, Mathe- 
matiker und Ingenieure, VI-+ 4068. m, 
147 Abb. Leipzig 1951. B. G. Teubner Verlagsgesell- 
schaft. Preis 22,30 DM. 

Das Buch hat sich die Aufgabe gestellt, dem in der 
Praxis stehenden Ingenieur, einen Überblick über die 
Lösungsmethoden von Problemen zu geben, die an 
ihn herantreten, und die ein mathematisches Wissen 


hi, 


300 Eingegangene Bücher — Nachrichten 


erfordern, das über seine von der Hochschule mitge- 
brachten Kenntnisse hinausgeht, insbesondere soweit 
es sich dabei um Randwertprobleme und ähnliche 
Fragen handelt. Das Buch ist also kein Lehrbuch. 
Es bringt nur das Wichtigste ohne großen Wert auf 
systematische Darstellung zu legen, Konvergenzbe- 
weise zu bringen usw. Es will nur eine Einführung in 
die Probleme und ihre Lösungsmethoden geben und 
durch ausführliche Schrifttumsnachweise den Weg zu 
eingehenderem Studium weisen. Die Kapitelüber- 
schriften mögen einen Anhalt dafür geben, welcher 
Stoff behandelt wurde: Randwertprobleme, Funk- 
tionen komplexer Veränderlicher, Partielle Ditfe- 
rentialgleichungen (das umfangreichste Kapitel, dem 
über ein Drittel des Buches gewidmet ist), Integral- 


ee 
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gleichungen, Variationsrechnung, Differenzenglei- 
chungen. In Anhängen werden einige mathematische 
Ergänzungen gebracht, die dem Ingenieur von der 
Hochschule her geläufig sein sollten. Im allgemeinen 
scheint mir die Stoffauswahl den Bedürfnissen des 
Ingenieurs Rechnung zu tragen, wenn auch der eine 
oder andere Wunsch offen bleibt. So scheint es mir 
z.B. wünschenswert, daß auf die Charakteristiken- 
theorie eingegangen wird. Die Darstellung ist gut 
lesbar und die vorausgesetzten Kenntnisse im allge- 
meinen denen des Kreises, für den das Buch geschrie- 
ben ist, angepaßt. Es wird manchem eine Hilfe sein 
und ihn zu eingehenderem Studium anregen. 
Willers. 


Dresden. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr.-Ing. Karl Klotter (o. Prof. a. d. Techn. Hoch- 
schule Karlsruhe), Technische Schwingungs- 
lehre. Bd. I: Einfache Schwinger und Schwin- 
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- Es muß in Formel 
(3) cos? # durch sin? ® 
(8b) 0,198 2/® durch 0,198 d-213 
0,043 5 durch 0,043 b-2 
ersetzt werden. Herbeck. 
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Es muß heißen: 
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8.163. 2. Zeile von unten: 2 = — —, 
Y 

8.164. 1. Zeile von oben: y= + 4 

8.164. 12. Zeile von unten: 


" we g = 
p’+2x9 + (2 000 24,57) 008 a y=0. 
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S. 164. Letzte Zeile: 


S.165. 17. Zeile von oben. Die erste Summe lautet: 
-Loo 
‘ PH a, eu. 
“=—o 


S. 165. 18. Zeile von oben. Beim ersten Summanden 
in der [ ] fehlt der Index x. 
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Math. Mech. 31° (1951), S. 131—142. 
Es muß heißen: \ 


S.131, I. Bezeichnungen: „=1— n 
v 
S.139, Formel (30): ee 
u een oe 
Ux % v % “x Re A m. 
Be 


u ih |. 
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